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Abstract. In this paper, we study the propagation of singularities (in the 
sense of C°° wave front set) of the solution of a model case initial-boundary 
value problem with glancing rays for a concave domain on an asymptotically 
anti-de Sitter manifold. The main result addresses the diffractive problem and 
establishes that there is no propagation of singularities into the shadow for 
the solution, i.e. the diffractive result for codimension-1 smooth boundary 
holds in this setting. The approach adopted is motivated by the work done 
for a conformally related diffractive model problem by Friedlander, in which 
an explicit solution was constructed using the Airy function. This work was 
later generalized by Melrose and by Taylor, via the method of parametrix con- 
struction. Our setting is a simple case of asympotically anti-de Sitter spaces, 
which are Lorentzian manifolds modeled on anti-de Sitter space at infinity 
but whose boundary are not totally geodesic (unlike the exact anti-de Sitter 
space). Most technical difficulties of the problem reduce to studying and con- 
structing a global resolvent for a semiclassical ODE on IR+, which at one end 
is a b-operator (in the sense of Melrose) while having a scattering behavior 
at infinity. We use different techniques near and infinity to analyze the lo- 
cal problems: near infinity we use local resolvent bounds and near zero we 
build a local semiclassical parametrix. After this step, the 'gluing' method by 
Datchev-Vasy serves to combine those local estimates to get the norm bound 
for the global resolvent. 
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Figure 1 . The rays containing wavefront set which hit the bound- 
ary transversaUy are reflected according to SneU's law i.e. with 
energy and tangential momentum conserved. 

Single incident ray hitting an edge (corner) generates a whole cone 
of reflected rays with angle given by the angle between the incident 
ray and the edge. 

1. Introduction 

In this paper, we study the propagation of singularities (in the sense of C°° wave 
front set) for the solution of an initial-boundary value problem in the presence of 
glancing rays on an asymptotically anti-de Sitter (AdS) manifold whose boundary 
is not totally geodesicQ Asympotically AdS spaces are Lorentzian manifolds mod- 
eled on anti-de Sitter space at infinity. 

Singularities of solutions of linear partial differential equations are described in 
terms of wave front sets of distributions. In the boundaryless setting, a special case 
of a result by Hormander |3] gives that for P a linear differential operator with 
real principal symbol p, the solution u of Pu — has its wavefront set contained 
in the characteristic set p~^{fi) and invariant under the flow generated by Hamil- 
tonian vector field Hp. In the presence of boundary, the rays containing wavefront 
set which hit the boundary transversaUy are reflected according to SneU's law, i.e. 
with energy and tangential momentum conserved (see figure [T]). See [1] for result 
on codimension-1 boundaries and e.g |16j for higher codimension boundaries. The 
situation becomes more complicated in the presence of boundary and glancing rays 
(i.e. tangency of the bicharacteristics). For codimension-1 smooth boundary, one 
has theorems of propagation of singularities along generalized broken bicharacter- 
istics (GBB)[^that include the diffractive result, by Melrose [7], Melrose-Sjostrand 
( |11) and ^12j) and Taylor [TS]. Their results rule out the propagation of singulari- 
ties into the 'shadow region', the bicharacteristically inaccessible region; in fact at 
a diffractive point, the ray carrying wavefront set does not stick to the boundary 
(see figure [2]) . This is not the case if we consider analytic singularitie^ a result 
given by Lebeau [S]. For settings with higher codimension boundaries, in the pres- 
ence of tangency, we have general theorems (i.e without diffractive improvement) 
of propagation of singularities along GBB for corners by Vasy fljB] and for edges 



^Totally geodesies: all of the geodesies of the manifold through points on the boundary are 
also geodesies in the boundary 

^ This definition of GBB distinguishes tangential rays according to their second order of tan- 
gency the boundary. If p is the principal symbol of the operator and (f) the defining function of 
the boundary, tangential rays are categorized(i.e where p = 0, </> = and Hp<l> = 0) according to 
the sign of Hp(j>, see for e.g. definition 24.3.7 in [2] 

^Gevrey-3 is the borderline case between the behavior associated with C°° singularities and 
that with analytic singularities, which is Gevrey-1. 
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Figure 2. DifFractive phenomenon in the case of tangency to a 
convex obstacle 

by Melrose- Wunsch-Vasy [14] using a relatively permissiv(^ notion of GBB. The 
diffractive problem, one aspect of which studies the phenomenon of propagation of 
singularities into the 'shadow region', remains open for higher codimension smooth 
settings. 

The interest in solving the diffractive problem on edges provides a motivation to 
study the diffractive problem on asymptotically AdS space, since asymptotically 
AdS space can be considered analogous to a 'reduced problem' on an edge but with 
'negative dimension'. In order to see this, we consider a simple example of an edge 
which is the product of a sector with an interval: 




Introduce polar coordinates (r, 6) on the sector, and let y be the variable parametriz- 
ing the edge. In (r, 6, y) coordinates the volume form on this space is : 

_ ^codimension of edge- 1 ^^^gi^^ 

When the codimension of the edge is 1, i.e. there is no 6*, we get back the smooth 
boundary setting. However, when the codimension is higher than 1, even if there 
is no dependence on 6 (e.g when the coefhcients of the metric are independent of 6 
and we solve the problem on fixed Fourier modes), the diffractive problem is quite 
different from the smooth boundary setting On the other hand, the volume form 

*This definition of GBB does not consider higher order of tangency of tangential rays at higher 
codimension boundary. The same definition was used for the analogous work done for analytic 
singularities by Lebeau [6]. See section 5 in '16' for more details on the definition and discussion 
of GBB 

^This is because of the term f codimension of edge— i ^jjg ^Qj^j-Qg form, which in turn affects 
the first order term in the Laplace-Beltrami operator. This is seen as follows: 
dVolg = ^|det g\dz, with z = (r,9,y), and 

A/dcto 2detg 
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on asymptotically anti-de Sitter space is 



^ ^— (total dimension of the ^P^^^) f^j-fly 

Viewed in this way, asymptotically anti-de Sitter space is analogous to a 'reduced]^ 
problem on edges with 'negative dimension'. In particular, the volume form for our 
model operator on M+ x equipped with an asymptotically anti-de Sitter metric 
(jlj is x~'''"'^^\l + x)~^/'^ dxdy, where n -I- 1 is the dimension of the space. 
Since we already have the general theorem (without diffractive improvements) for 
propagation of singularities by Vasy [T7] on asymptotically AdS spaces for both 
hyperbolic and glancing cases, we can study the diffractive problem for this setting. 
The metric for our model case is: 



(1) on M+ X M" : g = 



y.2 



For convenience, we work with a modification of the Klein-Gordon operator by a 
first order derivative 

This modification does not change the problem in an essential wajQ The diffractive 
condition for P is satisfied with HpX > where p is the principal symbol for 

conformal operator P ie p = — + [(1 + x)9^ — 

By Lemma 5.2 we have existence and uniquenes^ for the solution to: 



(2) 

(Pu = 

SUppW C {Un > 0} 

u = x'-Siy) + x'-+^gix,y) + h{x,y); (/(a;,?/) e C°°(M,,I?;);supp5 C {y,, = 0} 

h e H^'bioc' some a € M 



Here s± (A) = | ± y ^| — A are the indicial roots of the normal operator for P. 
The main goal of the paper is to study the singularity structure of the solution to 
([2]) and draw conclusions about the presence of singularities in the 'shadow region'. 
The approach used to study problem ^ is motivated by the work done for a con- 
formally related diffractive model problem by Friedlander [2], in which an explicit 
solution is constructed using the Airy function. The result is later greatly general- 
ized by Melrose using a parametrix construction in and [5] and by Taylor |15j . 
After taking the Fourier transform in y and denoting by 9 the dual variable of y, we 
need to construct a polyhomogeous conormal solution (modulo a smooth function 



^See remark 



1.5 



for a direction to generalize the problem which removes the 'reduced' condition 
^This is because the added term does not affect the normal operator and the 0-principal 
symbol. This also means the bicharacteristics remain unchanged. Hence the results (propagation 
of singularities and well-posedness) obtained for the Klein-Gordon operator -|- A by Vasy in 
[T7| still hold true for P 

^See also the proof of Lemma 1 



5.2 



(*). The term g is given in terms of S{y) by expression (|13 
definition (5.7) in ,17, 



for the motivation of the form of the 'boundary condition' 
For definition of HZh?" , a £ R see 

0,0, loc' 



in C°°(M""^^)) which satisfies the following boundary condition at a; = 0: 

(3) 

Lu e C°°(]R"+i) 



X "-u l^^o = 1 s±(A) = I ± Y - A 



u 



Here CphiC) is the set of polyhomogeneous conormal functions on some blown- 
up space C of x [0, l)a;. The chosen oscillatory behavior is modeled on that 
possessed by the specific solution Q Friedlander worked with in [51. 
The difRculties for our problem start with the fact that ^ under the Fourier 
transform does not have a well-understood solution as [5] does in terms of Airy 
functions. However, following the same change of variable as in [2|, the problem is 
reduced to studying the following semiclassical ODE on IR+: 

This family of semiclassical ODE is of regular singular type, hence a 'b'-operatoi|^ 
at the zero end and scatterin behavior at the infinity one. 

We use different techniques near zero and infinity to analyze the local problems: 
near infinity we use local resolvent bounds given by Vasy-Zworski '18', and near 
zero we build a local semiclassical parametrix. The next technical difficulties arise 
from the fact that the operator Q is not semiclassical elliptic, which prevents us 
from simply patching the two above local resolvents by a partition of unity, which 
otherwise results in an error that is not trivial semiclassically. To overcome this 
problem, we adopt the method of Datchev-Vasy [1] to get a global approximate 
resolvent with an 0{h°°) error. 

The remaining technical details of the problem involve constructing a local resolvent 
near zero for Q. For this we first need to construct a certain blown- up space of 
[0, 1)^ X [0, l)z' X [0, l)h- The usual method (see e.g. [13]) involves resolving the 
6-singularity at the corner of {z = z' ^ 0} and the semiclassical singularity. In 
our case, we incur non-uniformity in the behavior of the normal operator on the 
semiclassical front face as we approach the b front face (i.e. as z' — >■ 0). This 
gives motivation to do a blow-up at Z/j — )■ oo and z' 0, where Z = ^ is the 
variable associated with the b-blowup and Z^ — the variable associated with 
the blow-up along the lifted diagonal in the semiclassical face. The blown-up space 
we will work on results from the same idea of singularity resolution but with a 
different order of blowing up; namely, after the b-blowup we blow up the intersection 
of the b-front face with semiclassical face, then blow up the intersection of the 
lifted diagonal and the semiclassical face. There will be two additional blow-ups to 
desingularize the flow and to create a transition region from the scattering behavior 
to the boundary behavior prescribed by the indicial roots. We then construct a 



^Totally characteristic operators or 'b-operators' in the sense of Melrose |10j are differential 
operators generated by b vect or fie lds V;,, which is the set of all smooth vector fields tangent to 
all boundaries. See subsection 



3.4 



"in sense of Melrose scattering calculus [9] 



polyhomogeous cornomal function U on this blown-up space so that 



q,U - SKid e 4 ^C°°; SKid schwartz kernel of Id 

After removing singularities at the lifted diagonal, the remaining difficulty involves 
choosing the order in which one will remove singularities of the error to achieve the 
above effect. 

Once constructed, a parametrix provides among other things information on the 
properties of the fundamental solution, which enables one to make very precise 
statements on the propagation of singularities, in particular in the 'shadow regions' 
(see Corollary 1.4). We state below the main results of the paper, whose proofs are 
contained in subsection 113.31 

Theorem 1.1. [Wavefront set statement] 

Let U be the exact solution to to ^ on M+ x M.y then 

WF{U) C S 

where 

^ = {ix,y,^,v) : ^ ^ d^{y e - <P) ,r, ^ dy [ye-c^),de (zj • - 0) = 0, x > 0, |0„| > |^^'| > o] 

with (j) — Zg^^) sgn^n, 9 — {9' ,0n) being the dual variable of y under the 

Fourier transform. 

Theorem 1.2 (Conormality result 1). 

On yn > 0, if we are close enough to the boundary {x = 0}, then U G x^+C^y. 



As an immediate consequence of Theorem |1.2| we have 
Corollary 1.3 (Conormality result 2). 

On yn > 0, if we are close enough to the bou ndar y {x ~ 0}, then U is conormal 
relative to Hq, i.e. U G H^'j^ . (See subsection 3.4 for definition). 

As a result of Theorem |1.1| and Corollary |1.3t we have: 

Corollary 1.4. There are no singularities^ along the boundary or in the shadow 
region, i.e. figure^ holds for U . 

The paper is broadly organized into five parts. 

(1) Part 1 states the model problem together with the history, background and 
motivations of the problem. 

(2) The second part of the paper (sectionlH]) shows the reduction of the original 
problem to a semiclassical ODE on M^. To motivate the technical details 
that are required in solving the problem, section [7| gives an outline of the 
construction of an approximate solution to ([3]), i.e. the associated problem 
obtained after Fourier transform. The technical tools we will need are a 
local semiclassical parametrix construction and a global resolvent. 



The type of regularity measured at the boundary is . See subsection ( 3.4 1 for definitions 

and for further discussion 



projection ot 



boundary {t = 0} 




Figure 3. Singulartities are neither propagated into the shadow 
nor along the boundary 



(3) The third part is the technical heart of the paper where we carry out the 
construction of a local semiclassical parametrix near zero (in section ^) 



and the construction of a global resolvent (in section 10 1. 
(4) With all the provided technical ingredients, in the fourth part, in section 11 
we fill out the details of section [7] and construct the approximate solution 
to (H 



(5) In the last part (section 12 ) we study the wavefront set and conormality 
behavior at the boundary of the solution to the original problem ([2]) and 



prove Theorem 1 1.1[ Corollary 1 1 . 2 1 and Corollary 1 1.3 
Remark 1.5. [Open problems] 

Since asymptotically AdS can be considered a 'reduced' problem on edge (as explain 
in the section), one can generalize the current model problem ^ to one that is 
conformal to the problem on edges by adding the angular variable w. Specifically 
on K+ X M.y X SSuj where y = {yi, . . . , y„) we consider the following metric, which 
we call 'fibred'- AdS: 

-dx^ + {l + x)-^dyl-j:^lldy] , ^ 



n-l 



-Dg = (1 + x)x^dl^ " ^ ~ ^ ^" ~ ^ 2(a; + 1) ^^" " 

We can study the diffractive problem for 

x^dl -{n- l)xd, + A„ + x^[{l + x)9l ~ |0f ] + A 

Another direction to obtain new problems is to remove the translation invariance 
in the y variable of ([2]) . In the new problem, the coefficients will be dependent both 
on X and y so we cannot use the Fourier transform to reduce to a family of ODEs. 
This direction is analogous to the work done by Melrose ([7], [S]) to generalize the 
model diffractive problem of Friedlander in [5] . 
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(4) 



2. The model diffractive problem by Friedlander in 'layered medium' 

In this section, we briefly visit the model problem studied by Friedlander which 
were the motivation for the techniques used in our case. In |2], Friedlander studied 
the wavefront set of the solution to the following initial-boundary value problem: 

Theorem 2.1 ( Friedlander). On x Ry consider u the solution to 

.v^u d'^u d^u d^u 
dyl dyf dy^^^ dx^ 

U\x=0 = f 

^ u = ; y„ < 

then WF(C/) is contained in the union of forward bicharacteristic strips that em- 
anate from the conormal bundle of the origin over the domain {x > 0}. 

Note that he does not explicitly make a statement about singularities at the bound- 
ary. However from other results (e.g Melrose-Sjostrand [11], Taylor [15]) we know 
that there are no singularities at the boundary. Hence the following figure |3] holds. 
Friedlander's method was to construct the fundamental forward solution Kx to 

(5) 

yn < 

and study its wavefront set. To construct K^, he first takes the Fourier transform 
in y and looks for solution satisfying: 




(6) 



K L=o = 1 



9 is the dual variable to y. He then uses the following change of variable, which is 
valid as a smooth change of variables for 6'„ 7^ 

(7) x^c^e-^'^\e'\^~{i + x)e^J^ 

to get the Airy to equation 

(8) {dl^QF^Q-, kx^F{C) 

The Airy equation has as solution the Airy function with the properties: Ai((') is 
an entire function, whose zeros are on the negative real axis 

Ai(C)=exp(-|C'/')<f(C) 

_1_3, 

*J'(C) ^ 2^ «fcC ] uniformly as C -> 00 in {C : |arg C| < vr - e, < e < tt} 

k 

The unique solution to ^ is the inverse Fourier-Laplace transform of 

(9) kx{e) = Co = CU=o = e-'''\er - 

Ai(Co) 
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From here, after undoing the Fourier transform, it is standard to draw conclusion 
on the behavior of singularities of the actual solution in terms of its wavefront set. 
Melrose generalizes this work tremendously by using a parametrix construction and 
generalized to a class of Fourier- Airy integral operators in [7] , [5] . 



3. Asymptotically anti-de Sitter space 



3.1. Description of actual anti-de Sitter space. As a justification of termi- 
nology we include a brief discussion (borrowed from Vasy's paper |17| ) of actual 
anti-de Sitter space to show the form of the AdS metric in a collar neighborhood 
of the boundary. The form of the AdS metric is generalized in the definition of 
asymptotically AdS space in the following subsection. 
In E"+^, z = (zi, . . . , z„+i) consider the hyperboloid : 

2 I I 2 2 2 1 

^1 + • • • + ^ri-l ^ ~ ^n+1 — ^ 

living in the ambient space which is equipped with the pseudo-Riemanian 

metric of signature (2,n — 1) given by 

-dzl - ... - dzl^^ + dzl + dzl^^ 

Since -I- z^j_|_i > 0, we can introduce polar coordinates (i?, 9) in the two variables 
z„ and Zn+i. Then the hyperboloid is of the form: 

zl + ... + zl_,-R^ = -l 

inside W'-^ x M+ x S] and thus can be identified with x Sj. 

Next we compatify M"^^ to a ball 8"^^ using the inverse polar coordinates {x^u) 
with X = ^ in order to study the form of the metric close to the boundary. 
Denote by 

X =B^x §1 
A collar neigborhood of dX is identified with 

[0,1), xSr'xSe, 
where the Lorentzian metric takes the form 

_ -(1 + + [ - + (1 + x^)de^] 

^ x"^ 
The boundary is time-like and —dijj^ + {l+x'^)d9'^ is lorentzian of signature (1, n— 2). 




universal cover for 




Cok neighborhood : [0, l)i x SJ"^ x Rt 



Note AdS space is totally geodesic i.e. all of the geodesies of the manifold through 
points on the boundary are also geodesies in the boundary, hence we cannot talk 
about 'shadow region'. However, this property does not hold in general for arbitrary 
asymptotically AdS spaces, whose definition is given in the next section. 
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3.2. Definition of asymptotically anti-de Sitter space. These are Lorentzian 
manifolds modeled on anti-de Sitter space at infinity. In more details, they are 
manifolds with boundary X",n > 2, such that the interior X° is equipped with a 
pseudo-Riemannian metric g of signature (1, n — 1) which near the boundary X has 
the form 

-dx'^ + h 
9 = 2 ' 

X'^ 

here /i is a smooth symmetric 2-cotensor on X with respect to some product de- 
composition of X near Y , X = Y [Q,e)x] in addition, is a Lorentzian with 
signature (1, n — 2). 



3.3. Well-posedness results. A special case of the well-posedness results for so- 
lution of the Klein-Gordon equation on asymptotically AdS space by Vasy in [T7] , 

Theorem 3.1. [Vasy] Suppose f £ H^^^'l^^^ {X),m' < to e M and supp / C {j/„ > 

0}. A < Then there exists a unique U G H^'i^c which in fact lies in H^'^i^^{X) 
with 

{Dg + \)U ^ f,auppU Cl{yn>0} 

For all compact K <Z X there exists a compact K' <Z X and a constant C > such 
that 

Notation : For definition of ffo"b 'ioc' G M see page 2 in |17) . 

Theorem 3.2 (Vasy). X asymptotically anti-de Sitter space of dimension n. {Og + 
X)u = f where f € C°°(X) then 



u^x'+<-^^v,veC°^iX),s±{X) = ^^±y^^^^-A 
Note s±(A) are the indicial roots ofOg + A . 

Remark 3.3. In Friedlander's diffractive problem, the fundamental solution can 
be considered as a smooth function in x having value in distributions of y i.e. 
C°^(M.+ ^T>'y). However as indicated by the well-posedness results by Vasy quoted 
above, the solutions we will work with are not smooth up to the boundary and 
will have boundary behavior prescribed by the indicial roots s±{X). This will affect 
how one formulates the initial boundary value problem for our model operator (see 
section [5] in particular Lemma 5.1 and Lemma 5.2 1. 



3.4. Propagation of singularities. 

Theorem 3.4 (Vasy). On asymptotically AdS space, we have the propagation of 
singularities along 'generalized broken bicharacteristics ' (GBB) of g where g = x^g. 



See [T7] for more description of the GBB'''. 

In contrast with the boundaryless settings, the type of regularities one measures 
here is not C°° but 'b'-regularity conormal to the boundary (as explained next). 
This is as expected from the lack of smoothness for the solution of the Klein-Gordon 
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at the boundary (as noted in remark (3.3)). Instead oiC°°, the 'nice' functions are 



-^0 6° with: 

u e hI'^ ^ue Hl{X) and Qu G HI{X) Vg G Difff,(X) of any order 

Here, 

^tei^o^if Lu^ L^{X,x-^''+^Ug) VL e Diff^j^^) 

Diffb is the set of differential operators generated by Vb{X), which in coordinates 
have the form 



dx + ^jdyj '■ smooth vector fields tangent at the boundary 



Diffp is the set of differential operators generated by Vq{X), which in coordinates 
have the form 



dx + bj {xdy- ) : smooth vector fields vanishing at the boundary 



ax I 



4. Statement of our model diffractive problem on an asymptotically 

ANTi-DE Sitter space 

On M+ X where y = (yi, . . . , ?/„) for our model problem, we consider the metric 

^dx^ + {i+x)-^dvi-Y:;zldyj 

5= 2 

X'' 

For such a metric, J = |det g\^^^ = + x)"^/^. We consider the Beltrami- 

Laplace operator associated with g 

j=i '■^ 

For simpler computation in the model case, we will modify Dg by a first-order 
tangential derivative. We will work with the operator: 

(10) n-l 

= x'^dl -{n- X)xdx - (1 + x)x'dl^ + ^» + ^ 

As noted in the introduction section, this modification does not change the problem 
in an essential way. 

Taking the Fourier transform in y we obtain the operator L. Denote hy 6 = 
{01, ... , On) to be the dual variables to y under the Fourier transform. 

L := x'dl -{n- l)xdx + x^l + x)6l - |0f ] + A 
^^^^ __2rQ2 

X 

The principal symbol of x^^L is 



= x' [dl - (n - 1)^ + [(1 + x)9l \ef] + x-'x] 



iix,^,e)^-e + [{i + x)el-\ef] 
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5.2 



gives existence 



The difFractive condition is satisfied since i/?a; > Lemma 
and uniqueness to the following boundary problem: 

'Pu = 

supp u C {vn > 0} 

u e x'-6{y)+x'-+^gix,y) + h{x,y)- g G C°°(E„2?;); supp5 C {j/„ = 0} 

h e -ffg^Moc' some a € M 
Our main goal is to study the singularity structure of this solution in terms of 
wavefront set and investigate the presence of singularities in the 'shadow region'. 
In order to do this, after taking the Fourier transform in y, we need to construct 
a polyhomogeous conormal solution modulo smooth function in C°°(MI|.^"'^) which 
satisfies certain the boundary condition at a: = as follows: 
(12) 

Lu e C°°(M"+^) 

71 / 71^ 

2:^"-mU=o = i s±(a) - ± y ^ - a 

u e exp{-i^,^)CpHiC) 0in = [{l + x- l^f - (1 - |^'|2)3/2] sgn 0„ 

Cph{C) is the set of polyhomogeneous conormal functions on some blown-up space 
C of Me X [0, 1)2;. s±(A) come from the indicial roots of the operator and prescribes 
the asymptotic behavior of solution of the Klein Gordon on asymptotically AdS 
space (Vasy Tf). In our case, the dimension is n+ 1. This construction is carried 
out in section 11 and the singuralities are studied in section 12 and 13 
The bicharacteristics : the bicharacteristics are the same as those for the con- 
formal operator x^P, which is the operator considered by Friedlander ?7 (up to 
leading order). We will use Friedlander's description of bicharacteristics (see page 
146-148 of m). 



4.1. Geometry of singularity structure. There are similarities in the geometry 
of phase space T* (M+ x Wy) between our model problem and Friedlander's operator 
Q in [2]. Our operator P in ([2]) is a 0-differential operatoip^so its geometry lives in 
the 0-cosphere bundle^. On the other hand, there is a natural identification between 
0-cosphere bundle and the natural cosphere bundle. Up to an identification of the 
cotangent bundle near the boundary, the principal symbol for our model problem ([2| 
and that for Friedlander's problem [2] are both — -I- q{x, y, 9). Both problems will 
have the same geometry in the cotangent/cosphere bundle. In particular, as shown 
in figure |4] away from glancing, singularities are carried by a Lagrangian. The 
boundaries of the flow are the same in both cases: what happens at the boundary 
of the Lagrangian comes from the glancing point at the boundary of the domain 
from which it emanates (see figure [5]). However, the analytic objects that show 
up at the boundary will be different in each case. In Friedlander's case, this is 

l^On the characteristic set, |6»„| = \0'\ so e„ > 0. H^x = d^p = -2$. H?x = 2dJ = 20^. 
{x = 0,y = 0,1; = 0,9) with \6„\ = \9'\ > satisfies the difl'ractive condition: H^x = ^ = and 

Hfx = el>o 

^■^O-differential operators are operators generated by Vo{X) - smooth vector fields vanishing at 
the boundary. Vo{X) is also the set of all smooth sections of a vector bundle called ^TX, whose 
dual bundle is ^T*X. The quotient oi^T*X under the dilation action of R"*" gives the 0-cosphere 
bundle. 
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Lagianp coiitainiiij sinplarities Lagrangian containing sinplarities 
B in tlie interior at hyperbolic points \ 




Expect somekind of singularities 

at the glancing po int, ^ , 

/ ,,A-'f'^C' the bicharacteristics in ths ambient spate 

I which get there would have to bve the domain 
where the PDE hves 



Figure 4. In figure A, the problem has no boundary; initial data has 
singularities contained in the conormal bundle of the origin (the dotted 
cylinder is for visual comparison to the case with boundary). Figure 
B shows the structure of the fiow out from the conormal bundle of the 
diagonal. Figure C shows the flow out in the situation with boundary, 
where away from glancing, singularities are still carried by a Lagrangian 
while at glancing the Lagrangian has a boundary. 




away from glancing 

the set cair}ing sinpkities is Lagrangian 
Airy belia\'ior on the boundary 



Figure 5. Sideview of the problem with boundary (the PDE lives 
outside of the cylinder) and initial data having singularities contained 
in the conormal bundle at the origin 



associated to Kr^ in (|9|), while in our case, we will have to construct a family of 



-^asym 



.Ads that plays the same role as K^. 



5. Various results on well-posedness and formulation of boundary 

CONDITION PROBLEM ON ASYMPTOTICALLY AdS SPACE 

In this section we estabUsh the existence and uniqueness for solution to ([2|. The 
proof also explains the appearance of indicial roots s-(A) in the behavior of the 
solution at the boundary as stated in ([2]). 
Lemma 5.1. 

(1) For g £ Vy there exists u± £ x*±C°°(M-|-,2?(,) satisfying for some a G R: 



l,a 



(2) For g e V'y there exists u± e x''±C°°{R+,V'y) satisfying 

; C°°(M„2?;)=x°°C°°(M+,2?;) 



Pu± eC°°{R^,V'y) _ ■ _ ^oo^oo^ro + 

u±ex-^g±+x-i+^C°-(R+,V') ' ^^-'^1/^ 



Proof. 

Remark 5.1. By a special case of definition of h'^'^™ ,ni > 0,r7i e M in [17], if 
w = E QjVj where e H^'^{X),Q, e DiS^{X),H^f = then v G 'fT™- i^or 
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general definition of H^'J^^ ,m g M see definition (5.7) in [17j). Hence ipj € h[ 



0,b 

for some aj > 0, G 

Remark 5.2. 5 G Vy the space of tempered distributions ie g is a finite sum of 
derivatives (with constant coefficients) in y of continuous functions. And derivatives 
in y commute with P. Hence the proof for arbitrary g is completely analogous to 
that for g = 5{y). 

(1) Consider (^0 G We consider ipj G Vy being b derivatives of ipo. With 
such (pj-s, consider ii of the form: 

u = x'-ipo + ^x''-+''ipk{y) 

k=l 

Write P as 

P^Qi+x^Q2 

n-l 

where Qi = {xd^)'^ - nxd^ + A; Q2 = -(1 + a;)a:^^9^„ + a;^ ^ 9^^. 



Note that Qi has indicial roots s± with s± (A) = ^ ± y ^^^^ — A then PU is 
of the form 

fco-l 

Pu = X'-+\2S^+ 1)^1 x'-+''[i2s^+k)k^k+Q2^k-2] +x'~+''° + ^Q2^ko-l 

k=2 

With fco large enough then the last term x^-^^°^^Q2'^kQ G ^o'moc ^o^' 
some a G M If we only care of the decay in x then this term is in 
2,s-+fco+ic°o(M^,X'^). 

We can solve for ipj with 7 = 1, . . . , fco — 1 so that Pu G Hq'i^\oc by requiring 
(13) ipi^O; For j = 2, . . . , fco - 1 : (/Jfc = ^ ^^j^ Q2(Pk-2 

Hence we have constructed u satisfying 

{u — x^- ip{y) + 2;'*-+-'^ Ym:'=q x'^fk', 'Pk being b-derivatives in y of Lp 

Since Hl^l^^^ C H'lf^^, Pu G H'lf^^. 
(2) Let is such that (A^^ + l)'"(p = 5{y) for some m. We have (Ay + 1)™ 
commutes with P. Also suppPu C > 0}. Define 

Itbdy = (Ay + 1)"U 

then Ubdy satisfies for some a G M: 
fPUbdy G -ffo'bjoc 

\wbdy G X'- (Ay + l)^^ + x'- + ' Y.lU' ^H^y + ^rVk 



^^fco is large enough to cancel out the weight in the volume form L'^{x '""'"^'dg). With this 
choice of fco, the term in consideration is then some b-derivatives of a function in Lq. Note that 



Q2 consists of only b derivatives (in fact only derivatives in y). By remark 
some a S M 
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5.1 



it is in ^^q'^ for 



where (pk are derivatives h m y oi and thus in V'y. Thus we obtain the 
first statement of the lemma. 



(3) Instead of stopping at fco as in ( 13 1 we can solve in series and ask that 



^1 = 0; For j = 2, . . . , : = (^^7^Q2^fe-2 
We then take an asymptotic summation to get itbdy;oo, which satisfies: 

\wbdy;oo e X'- (Ay + 1)™(^ + x'-+^C°°{K,,V'y) 

(4) The same proof goes through if we replace s_ by s+. 

□ 

As a corollary, we can formulate a boundary problem for our model operator: 
Lemma 5.2. We have the existence and uniqueness for the solution to 
' Pu = 

supp?i C {vn > 0} 

u e x'-S{y)+x'-+^gix,y) + h{x,y); g G C°°(R,, P^); suppg C {y„ - 0} 

h e H^'hAoc^ for some a G M 

Proof. Denote by / = Pitbdy where Wbdy is as constructed in the proof for Lemma 



5.1 



For some a G M, / G By Theorem 3.1 (Vasy), we have the existence 

and uniqueness of U G ^q'^'ioc some a G M with U satisfying: 

fpu = / 

[suppU c {yn > 0} 



Lets remind ourselves of the properties of Ubdy from the proof of Lemma 5.1 

feo-l 

Ubdy = (Ay + l)'"it; u = x^-ifiiy) + x"-'^'^ ^ x''(pk 

fc=0 

where tpk are derivatives in y of ip and ip is such that (A^ + 1)™</j = 5{y) for some 
m. Hence (Ay + l)™ipk are also derivatives of S{y). Thus suppubdy C {yn — 0}. 
Let u be defined by 

U := Ubdy - U 

then u as constructed solves : 
' Pu = 

u G x^-S{y) + x-+iC-(M,,2?;) + ffi;,"!,, 
^suppu C {yn > 0} 

Uniqueness follows from the fact the g is uniquely determined as shown in the proof 
for Lcmma|5.1| There, we construct as a Taylor series in x with coefficients given 



explicitly in terms of S{y) by (13). □ 
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6. Reduction to the semiclassical problem on M+ 



As described in section[2j in "5] Friedlander studies a diffractive problem conformally 
related to ours After the Fourier transform in the tangential variables, he uses 
a change of variable Q to get the Airy equation. Although we will not get the 
Airy equation, using the same change of variable, we show in this section how our 
model problem is reduced to the following semiclassical ODE on ]R+ of the form: 

This family is at one end of regular singular type (hence a b-operator in the sense 

of Melrose) while at infinity has a scattering behavior. 

Step 1: Follow the first step of Friedlander and use the change of variable 

x^Z = -\en\-^/^\e'\^ + (1 + a;)|0„|2/3 e M 

Also denote by 

(14) Zo = -\en\-^/^\e'\^ + |^?„p/3 eR^Zo = z\,=o 

Denote 

Q = x-^'/^Lx"/^ = {xd^f + x^[{l + x)9l - + A - Y 
In terms of Z, Zq , Q becomes 

(15) Q = {Z- Zofdl + {Z- Zo)dz + \--^ + Z{Z- Zof 

Geometric interpretation of Zq: Zq = |6'nP^^(l— ]w\z) ■ nonhomogeneous 
blowing up of the corner in phase space at glancing |0„| — \9'\ and infinity 
(see figure [6]) . The functions in terms x and 9 live on the blowup space B as given 
in figure |6] 

Step 2: With S ^ Z - Zq 

(16) Q = S''dl + Sds + X-'-^ + S^ + ZaS^; 

This is an family of ODE in S G [0, oo) with parameter Zq G {—oo, oo). 

Goal: Need to know the behavior of the solution at S = and S^^ — for bounded 

\Zq\, uniformly as Zq — >■ oo and uniformly as Zq — > — oo. 

Step 3: We rewrite the problem as Zq — ^ oo and Zq — > — oo as semiclassical 
problem using the following variables z,h. 



(17) z := xo{Zo){Z Zo) + x+(Zo)(^-^) + X-(^o)(^^) 

where XojX+jX- smooth function on M such that for some 62 > 0, 

Xo = 1 on \Za\ < ^2 fx+ = 1 on Zq > 2S2 . fx- = 1 on Zq < -2S2 

suppxo C {|Zo| < 2^2} ' \suppx+ C {Zo > (52} ' [suppx- C {Zq < -(^2} 
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Figure 6. Zq is obtained by blowing up the corner in phase space 
at glancing \e„\ = \e'\ and infinity \e\-^ = 0. Fi gure on the left 
show the base of the blowup space B. Figure on the right shows 
the portion of the blow-up space B corresponding to 9 > 0. Both 
portions for 9 > and 9 < are similar. 



On the two more interesting regions, 

^ ('=.x(l-|^f )-i Zo>2d2 
^ \=xi\9f-l)-' Zo<-2~52' 

Similarly define 

(18) h xo{Zo)Zo + x+{Zo)Zo''^^ + X-{Zo){-Zof^^ 
On the two more interesting regions, 

\(-Zo)-3/2 Zo<-2~62' 

We will also use 

(19) a := x.{Z,){Z Zof^' + X+{Z,){^-^f'' + X-{Zo){^ ' 



Thus, on the two more interesting regions, a — 



Zo ' 

X3/2(|0'|2 _ l)-3/2 < _2^2 



The space the solution lives on blowup space C, see figure [7| 

Remark 6.1 (On the blow-up at z and /i = : transition region). In the finite 
regime, the boundary behavior (at 2: = 0) is transcribed by the indicial roots 
depending on A; while the boundary behavior (i.e. at 2; = 0) of the solution of 
the transport equations is independent of A. Hence a transition region is needed 
between the oscillatory behavior in the infinite regime and the regular singular 
behavior (with no oscillation) in the finite one. See figure |8] 
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Figure 7. The pictures show the portion of blowup space C and 
blown- up space B for 0„ > with the fibration structure {Zq — 
const, CT = const} over front face F^- This fibration structure will 
be used to extend the solution on into the interior of the space, 
(z = a2/3) 

The other portion for |6'„| < is similar. 

The finite region is labeled by I = {|0| < 2}. For picture of the 
base (in terms of \9\) see figure [6] 



lb,j 




transition area 



tranport equation 
legulEtr singular down to u = 



oscillatorj' exponential 
■starts to take effect 



regular singular and no osr.illation 
b 'I'll 



Figure 8. The blowup is associated to z = 0, /i = 0. The projec- 
tive coordinates are t = j^,h (In the picture t — /i^/'^,cr — z^^^. 
See (171 and (18) for definition of z and h ) 



We will need to use the following variables associated to these blowup. The bound- 
ary defining function for F^-* 
(20) 

= t := x.{Z,){Z Z,) + x^(Zo){^-^) Z'.^' + (^o) (^-^) (-Zo)^/^ 

C- Zo — Zq 
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Properties of Q in terms in z, h 

Consider the region on wlriclr > 25-2. On this region 



z -.^^ {Z - Zo)Zq^; h 



-3/2 



In these coordinates Q has the form: 
(21) 



Q^ 



h'^{zd,y + {X- 



'-)h^ + z^ + z' 



t=zz;" ^-2 



4/3 



This is a family of semiclassical ODE with 'bj^ / regular singular behavior at one 
end (at z — Q) and scatterin behavior at infinity. 

Remark 6.2. In terms of /i, z the solution used by Friedlander on this region has 
the form 



Ai(C) 



exp 



( - fz/l-l [{Z + 1)3/2 _ 1]) + + /); / g C^ih, Z) 



Ai(Co) 

This is the model of phase function we will use for the semiclassical local construc- 
tion and the construction of local approximate solution (for infinity region \9\ > 1). 



As for Zo < we have another semiclassical problem however this time much 

simpler. 

(22) 

Q_ = h-^ \h^{zd,)^ + (A - '-^)h^ -z^ + z^l 



Semiclassical b- symbol : —^^ 

isolated point at z = 0. 

Local behavior at infinity : as in the case > and we can use the result of 
Vasy - Zworski to obtain semiclassical local resolvent estimate. 

Remark 6.3. The 'phase function' of — has the form 



- Ai(Co) 



exp I 



2(^3/2 



2-^ Jexp ( - i - |/i i(z- 1)3/2 sgnein) exp(-/i-i) ^ ^ 



exp 



(_|/^-l[l_(l_ ^2/3)3/2]^ 



Z < 1 



7. Construction of approximate solution on the blownup space C : 

Main idea 

As motivated in section |4] we need to construct a polyhomogeous conormal solu- 
tion modulo smooth function in C°°(M"^^) which satisfies the following boundary 



^'^Totally characteristic operators or 'b-operators' in the sense of Melrose I10| are differential 
operators generated by b vect or fie lds V;,, which is the set of all smooth vector fields tangent to 
all boundaries. See subsection 



3.4 

"'^^in sense of Melrose scattering calculus [9] 
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Figure 9. Front face ¥c of blowup space C (see figure [?]). This 
represents the resolution of the infinite regime \9\ > 1 at glancing 
behavior |0„| = \6'\. The number i indicates the boundary face 
F^^ for i = 1,...,6. 

F^°|_ denotes the neighborhood in region ^ close to the corner 
of i = 0, /i 

is the neighborhood close to z 

^"^^ and F^;" 



lO.oo 



where t is the projective blow-up with i = f ■ F^ ^ 

0. Similar notation for 



o,r 

F^°° and F^° for the side of ¥c where Zq < 0. Denote by F° the 



region of ¥c away from boundary face F, 



See ( 20 1 for the definition of , 



(4) 

a 



5,(1) and (17) for the definition of z. 



condition at x — 0: 

' Lu e C°°(M';+^) 

' x-^-ii = 1 s±(A) ^ -±d — - \ 

u e exp(-z0i„)/:p,,(C) 0i„ = le-'[{l + x- |^'|2)3/2 _ (1 _ |^'|2)3/2] sgn^^^ 

We show in section |6] that this problem has a reduction to a family of semiclassical 
ODE on M''" , living on the front face F^ ( see figure [9] ) of the blownup space C 
where the solution in terms of {x,9) lives (see also figure [6] and [?]) . The details of 
the construction are in section [TT] In this section, we list the main ideas in order 
to motivate the technical tools involved: a global resolvent and a local parametrix. 
We construct an exact solution for Lu = near infinity (16*1 > 1) and cut this off 
away from \9\ < 1, which now gives an approximate solution to ([3|. To achieve 
this, one needs an exact polyhomogeous cornomal solution on F^ (see figure [7]) i.e. 
with no error since this function will be extended into the interior of blowup space 
C using the fibration structure over front face F^. This gives a solution with the 
required asymptotic behavior at infinity where \6\ > 1. 

Partition the front face Fc into Fc,o, ^c,+ = ^c,+ U ^c,+ and ¥c- = F^ _ U Fg^_ 
as shown in figure [9] This construction is carried out in two steps: 



STEP 1: First we construct an approximate solution that is polyhomogeous conor- 
mal to all boundaries of this face, with an error vanishing to infinite order at all 
boundaries of Fp. This is simpler for the interior where the fibers are not singu- 
lar. In the region where \Zq\ — > 0, wc have the option of constructing directly a 
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polyhomogeneous conormal approximate solution by first constructing it on region 
^ , _ then solving along the global transport equations away to infinity. In 
fact, as we see very soon that since we will have to build a local parametrix for 
this region anyway (ie for cr ~ 0), we will use this option to obtain the approxi- 
mate solution from the parametrix by taking the limit as z' — > of the parametrix 
(which is a function of h, z and z') after some rescaling. After this, we will solve the 



error along the global transport equations away to infinity. Also see remark (7.1) 
for further discussion on the difference between the two chosen methods in dealing 
with region near z = and near z — >■ oo. 

Denoted by L^^ the constructed function in this manner. QL^^ — Eh vanishes to 
infinite order at all boundary faces of Fc. 

Remark 7.1 (On the difference between the methods used for construction of 
approximate solution near z = and that near z — ^ oo). We only construct a local 
parametrix for the local region z ~ and not for the infinite region (see remark 



(7.2) for further explication). Because of this, we can read off readily information 
of the local solution from the parametrix, while having to work out directly the 
asymptotics of the approximate solution for the infinite region. 



STEP 2 : In this step, we solve the error Eh exactly by using a global resolvent 
on Fc. Denoted by Uf^ the exact solution on Fp. 

We use the gluing method by Datchev-Vasy [T] to obtain the global resolvent G 
on Fc with a polynomial bounded in h weighted norm. 11^^ a-nd a-re then 
different by 

— Ufc ~ GEh 



where the mapping property of G given by Lemma 10.1 G maps C°°(Fc) into the 
space of oscillatory functions which have the same asymptotic behavior as t/p^ and 
GEh = 0{h°°). As a result of this, Uf^ has the same asymptotic properties as C/fc 
does at infinite \6\. 

Remark 7.2. [Ingredients and methods to obtain the global resolvent] Since our 
operator is not semiclassical elliptic, one cannot simply patch together by partition 
of unity results from the local region around zero and the infinite region. To 
overcome this problem, we use the gluing method by Datchev-Vasy In order 
to apply this microlocal gluing method, on each region, one either needs a local 
resolvent that is polynomially bounded in h or at least a local parametrix with a 
semiclassically trivial error. 

For the infinity region, we have the semiclassical local resolvent estimate by Vasy 
- Zworski [TB] for asymptotically euclidean scattering settings. On the other hand, 
the local problem near z = is more degenerate due to the 5-geometry of the 



bicharacteristics as shown in figure 10 and we do not have results on resolvent 
estimate for this region. Hence we will have to construct a local parametrix for the 
region close to zero. In fact, this situation gives rise to a new research direction: to 
obtain semiclassical resolvent bound directly without a parametrix construction. 
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Figure 10. The semiclassical b- symbol for Q is — + z"^ + z^. 
The semiclassical characteristic set is singular at z = 0. 



8. Some preparations needed for the semiclassical parametrix 

construction 



Recall the family of ODEs on M+ obtained from the model problem P through 
various changes of variable in section [sj which we denoted by Q where Q — h^Q, 

q+ = h'{zd,)' + h^{\-'^) + z' + z' 



Q_=h^{zd,f + h^{X-—) + z- 



3 z2 



Q-i_ corresponds to Zq > and Q_ to Zq < 0. For definition of Zq, h and z see 



(141, (18) and (17) 



8.1. Phase function. In terms of h, z the solution used by Friedlander on this 
region has the form 

Ai(C) 



exp 



( - fz/l-l [{Z + 1)3/2 _ 1]) + 1)-1/2(1 + /);/£ C-(/i, Z) 



Ai(Co) 

This is the phase function we will use 

= h-'\^iz) - (p(z')l; viz) + 1)'^' - 1] sgn0„ 

This function can also be obtained directly from solving the eikonal equation for Q 

this is in b blow-up coordinates : Z — p = z' ,h. 

■V 

Consider the conjugation of Q;, by eh ; First we compute : 
Zdz [e'hu] = ih'^e'huZdzf + e'hZdzu 

[Zdzf (e'lu^ = -h-^e'hu{Zdz^f+ih-^ehu{Zdz)^^+2ih-^e'h{Zdzu){Zdz<p)+e'h{Zdzy' 
The eikonal equation comes from coefficients of h'^ : 

u[-iZdz^f + {Zpf + {Zpr]^0 

=>dzip^±[Zp'+p'f' 



^ = ± 

=^ (y3(Z,p) = ± 



1/2 



2 [Zp3 + p 

3 p3 
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dZ + C{p) 

21 3/2 

— + C{p) 



We can check again that 



i3/2 



solves -{Zdzvf + {^pf + (Zp)^ = 

In semiclassical blow-up coordinates : Zh = , h, p then Z = /iZ/, + 1 . In 
this coordinate the phase function is 



±-[(hZhP + p+lf'^ + C{p) 



2 



ip+lf/^h-' 



±{l + hZh 



P + 1 



3/2 



C{P) 



which we need to be equal to ±pZh when restricted to ft, = 



h->-0 6 



± 1 + hZh 



P 



3/2 



+ |c7(p)(p+l)-^/^ 



±Zh{p''+p 



3-L.2X1/2 



Hence we need 



In that case then each fixed Zh and p the limit is 



lim^(p+l)3/2ft-i 



± 1 + /l^ft 



p+1 



3/2 



T^C(p)(p + l)-3/2 



= ±i(p + 1)3/^3^,^ = ±Z,p(p + = ±Z,(p3 + p2)l/2 

3 P + 1 

The phase function should be symmetric in z,z' and singular across the diagonal. 
Hence we choose to work with the following pheise function, which is shown in 

various coordinates : 

• In coordinates on the double space crossed with h : z, z', h 
2, 



'h-'\{z + lf/^-{z' + l) 



3/2 1 



The phase —ipin satisfies the cikonal equation 
-{zd,f + z^ + z^ 

for the operator 

q.{K\) = z-'i[-ziQAz^ =h^{zd,f^h^{\-^)^z^ + z^ 



In b blow-up coordinates Z = ■^,p = z',h 



In b semiclassical blow-up Z^, 



z-i 
h ' 



p = z',h 



^in=^h-'\{hZhP+l + pf/'-{p+ir/'\ 
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8.1.1. Computation to show that the phase function (p is smooth up to *8i in ^b,h- 
Coordinates associated to the various blow-ups in this region 

• Away from the semiclassical face A, on the projective coordinates are 

fi' = ^,h,Z 

Away from the b front face, the projective coordinates on £01 are : 

i>' = ^,z',Z 

• On DJl, away from the semiclassical face A, the projective coordinates are 

jl=lh,Z' 

while away from the b front face, the projective coordinates are 



Away from the left boundary £ : The coordinates in this region are 

z' 
h 

in which the phase function takes the form : 
2 



-(Z/i'/j+ 1)3/2 + + 1)3/2 



~ 3^ {Zfl'h + 1)3/2 + (^//j + 1)3/2 

_ 2 + 3{il'fh + {fi'fh^] - [3Z/i' + 3Z^{il'fh + Z3(/i')3/i2] 

~ 3 {Zjl'h + 1)3/2 + (^//j + 1)3/2 

since Z, jj.',h>0 the above expression is smooth in jl', h, Z. 

Away from the right boundary TZ : The coordinates on this region are 

M = r,h,Z' 
h 

in which the phase function takes the form: 



+ 1)3/2 1)3/2 



(/x/lZ' + 1)3/2 + 1)3/2 



{zZ' + 1)3/2 _ (2 + 1) 



3/2 



_ 2 _i (/i/iZ')3 + 3{fihZ'f + SjlhZ' - [iilhf + iijlhf + 3/i/i] 

~ 3^ (/},/lZ' + 1)3/2 + (/t/i + 1)3/2 

_ 2 {flZ'fh'' + 3{flZ' fh + /iZ' - [(/i)3fe2 + 3(^)2/, ^ ^] 

~ 3 ijlhZ' + 1)3/2 _^ (^^Jj + 1)3/2 

since Z', /i, /i > the above expression is smooth in /i, /i, Z'. 

RemEirk 8.1. The phase function is smooth up to *8i but blows up at the semi- 
classical face. 
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Figure 11. Blown-up space VJlb.h at the Zq > end 




.2. The blow-up space. 



Denote by X = [0, 1). We have two semiclassical problems when Zq — > oo and Zq — > 
— cx) ie when /i — > (defined by (18)). Start with the space X Xh X x [— cx), cxjj^o, 
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where X XhX is the b-blown up space which is obtained from X x X hy blowing up 
the corner dx x 9x ( which in coordinates is {z — z' — 0}) in X x X . In this section, 
we wiU analyse separately each end but the parametrix will be built all at once on 
the whole blow-up space. Locally each end can be considered as X Xi, X x [0,oo);i. 

Notation : Denote 



• by £ the lift to X XtX x [-00, oo]zo of W = 0}, by 7^ of {z = 0}, 

• by the semiclassical face A the lift of {Z^^ = 0, Zq > 0} { — {h — 0} on 
the Zq > end), 

• by the semiclassical face A the lift of {Z^^^ = 0, Zq < 0} { ^ {h = 0} on 
the Zq <0 end), 

• by Diagf, ^ the lift of the diagonal {z — z'}, 

• and by J- the resulting blown-up b - front face. 

The usual approach: In order to resolve the singularity of the phase function at 
the semiclassical face we have to do these following two blow-ups: the intersection 
of the b front face with the semiclassical face, and the intersection of the lifted 
diagonal with the semiclassical face. The two operations however do not commute 
with each other; different spaces result, depending on the order in which the blow- 
ups is carried out. For example in jl3j . the authors blow-up the intersection of 
the lifted diagonal and the semiclassical face in X Xi, X x [0, oo)h, whose resulting 
front face is called the semiclassical front face S. The operator in [T3] has uniform 
property on S all the way down to the intersection of S and J^. However in our 
situation, if taking this approach, we have degeneracy as we approach J- along S. 
This degeneracy requires further blow up at the intersection of S and 
In particular, the normal operator on the semiclassical front face has the form for 
Zo > 0, 

dl + izr + izT; Z, = ^; z = ^ 

We have degeneracy of the potential and the behavior is not uniform as we approach 
the b front face (i.e. as z' — >■ 0) along the semiclassical face. Denote by A = 
{z')^ + (z')^. At arbitrary level A 7^ solution has the asymptotic form: 

exp(-^VA|Z^|) C°°. 

This gives motivation to do a blow-up at Z^, — > 00 and A — )■ i.e. z' — 0. We also 
have degeneracy when Zq < with the normal operator on the semiclassical front 
face of the form: 

dl + i^r-i^T, Z. = ^; Z^f 



The blow-up space in our setting (see figure 12): The blow-up space we will 
use is obtained from these following blow-ups. For both ends Zq > and Zq < 0, 
we will first blow up the intersection of T and A, and the intersection of and A 
in X Xt, X X [—00, oo]^^, whose resulting new front face is denoted by *Bi and *8i 
respectively. Next on this new space, we will blow-up the intersection of the lifted 
diagonal and the semiclassical face A and similarly the intersection of the lifted 
diagonal and the other semiclassical face A; denote respectively by and ^82 the 
resulting blown-up front faces. Call our newly constructed blown-up space 9Jl. (See 



figure 11 1 



We will need to carry out two additional blow-ups on ^Sl. On Zq < 0, the problem 
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bessel-tj'pe 




Figure 13. We have the Bessel type ODE in t along the flow 
{v' = const}. *Bi is a blown-up version of the product structure 
shown on the right x Kj^", . The label 'bessel-type' and 'not 
bessel-type' denote the behavior of the limiting operator on the 
corresponding boundary sides of *8i 



is elliptic there (so singularities do not get propagated) ^ these following blow-ups 
are more of formal nature. However for Zq > 0, due to the hyperbolicity of problem 
{hence propagation of singlarities) , these procedures will be necessary. To desin- 
gularize the flow on *Bi, we will blow-up the corner intersection C C\ A, denoted 
the resulting front face by ^i. On the other hand, the transport equation along 
semiclassical face A does not give rise to the boundary asymptotic at TZ dependent 
on the indicial roots. Thus, we will blow up the corner intersection at TZH A and 
denote by ^2 the resulting front face. ^2 serves as the transition region from the 
scattering behavior to the boundary behavior prescribed by indicial roots by one 
solving the Bessel equation there. See the computation at the end of the subsection 
for more details. For Zq < 0, denote by the and ^2 the front faces resulting from 



the same blow-up procedures. Denote by dJtb^h the resulting space (see figure 12 1. 

Denote by Dyti,,h the space obtained from ^b,h by blowing down ^2 and ^2- 

Computation related to the blowup of and ^2- 

The equation in z variable: 

In projective blow-up coordinates on X^/i' — j^,h so z' = jlh, the operator 
has the form: 
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In coordinates /i', Z, z' on M2 (away from the left boundary £) the operator be- 
comes: 

^2 



(z9z)' + (A--) + MZA'r + (ZA'r 



This a family of ODE varying smoothly in h. In the neighborhood where these 
coordinates are valid, the boundary face *Bi can be described by {h = 0}. Hence 
when restricted to h = i.e. the front face *Bi, we get 

(Z9zf + (A-^) + (Z/ir 

The flow on which we have the ODE is along {//' — constant}. This flow becomes 
singular as one approaches the corner where *Bi meets C and semiclassical face A. 
To desingularize the flow, we are going to blow up the edge which is the intersection 
£ n A in 93t. In addition to this, we need to resolve the behavior at Z = and 
i>' = if^')^^ = 0. The projective coordinates to describe this blow-up: 

v' 

For the computation for Zo, we just need to change the sign in front of ic the 
operator on *8i is 

(ZSzf + (A-^)-(ZA')^. 

8.3. Operator space. We follow closely the definition used in [T3], however the 
resulting operator spaces are not the same since we are not on the same blow-up 
spaces^ Denote by 4'™,^ the space of operators P whose kernel 

Kp{z,z\h)\da(.z')\ 

lifts to '^h,h a cornomal distribution of order m to the lifted diagonal and vanishes 
to infinite order at all faces, except the b front face T ^ up to which it is smooth 
(with values in conormal distributions) and the semiclassical front faces *Bi and 
*82, up to which it is h~^C°° . Similarly behavior for *Bi and 032- 

9. Semiclassical left parametrix construction for Q(/i, A) for 

< z,z' < 1 



As motivated in section [t] (see in particular remark (7.2)), we need to build a local 



parametrix for Q. For the parametrix construction, we will work with Q where 
Q = /i^Q, so Q is of the form: 

Q+ - h'izd^f + h^{\ -"^j + z' + z' 

Q.=h'{zd.f + h^{X-'^) + z^~z^ 

Q+ corresponds to > and Q_ to Zq < 0. For definition of Zq, h and z 
see (14), (18) and (17). This family of semiclassical operators Q is obtained from 
our original operator ([2| by a fourier transform and several changes of variable 



See explanation in subsection (8.2i 
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(in section |6] in particular equation (22)). As motivated in section [7j we use the 
local parametrix for two purposes. First, we obtain from the parametrix an local 
(near zero) exact solution for Qu = 0. The second (in fact main) purpose of the 
parametrix is obtained a global resolvent in order to solve exactly semiclassically 
trivial errors. 

The main method of the construction is as follows: We construct U on dJlb.h 
such that 



QU - SKid e h'^z'^iz')y 4 
where SKjd denotes the Schwartz kernel of the identity operator Id. We will be 
removing error the boundary faces of OJtbji in the following order: 

(1) Step 1: We first remove diagonal singularities by extending the lifted 
diagonal past the boundary faces dJl^.h at both ends Zq > and Zq < 
and use transversal elliptic construction. 

(2) Step 2: next we remove error at the semiclassical front faces 032 and *B2. 
The normal operator in local coordinates has the form 

this is the Euclidean laplacian with spectral parameter z' + 1. Note the 
uniform behavior down to *Bi[^ For Zq < 

M^^{Cln)^{z'f [dl^+z'-l]. 

After this step we can blow-down *B2 and S2 and work on Mb^/j. 

Remark 9.1. At the Zq < 0, we have an elliptic problem. Unlike when 
Zq > 0, for Zq < 0, we will not have singularities propagated from *B2 onto 
A, thus steps 3, 4, 5 and 6b are only needed for Zq > 0. 

Remark 9.2. [Providing the motivation for the next steps from 
step 3 to 5] Step 3, 4 and 5 served to blow down *Bi and *8i to a prod- 
uct structure. In more details, on 58 1, Qb restricts to : (see step 4 for 
computation) 

which thus can be considered to live on the product space Bi — x Rf, 



(so Bi is the 2-dimensional product structure represented in figure ( 13 ) by 
the rectangle). On the other hand, the error we need to solve away (ie the 
inhomogeneous part) lives on (i.e is nice' function on) We will replace 
this error by something that behaves trivially at ACi and J" n *Bi in 
order to be work with an error that resides Bi. As for the Zq < end ie 
on 03 1, the normal operator is 

(%)2 + (A-^)-P 

As noted by the previous remark due to the ellipticity of the problem we 
automatically have infinite vanishing at An^i after step 2. 

(3) Step 3: for this step we remove error at the semiclassical front face A and 
5^2- The process is carried out in 3 smaller steps: 



This is the reason w hy w c choose to work with our current blowup space instead of one used 
in e.g |13| . See subsection 8.2 for further discussion 
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..A 

Initial condition at {Z = 1} 
for the transport equation in Z 




Projection of Lagrangian Li 



Figure 14. Lagrangian Fi associated with phase function (pi 
Note that the integral curves never hit the left boundary £ 



Step 3a: Singularities are propagated out from *B2 onto the semi- 
classical face A. We remove the error at A along the Lagrangian Fi 
associated to the phase function — yjin where 



+ 1)3/2 -(z + 1) 



3/21 



by solving the transport equations with initial condition at the diago- 
nal. Note that the integral curves never hit the left boundary C, see 
figure (14 1. In the region close to .4 n ^2 we have z < z' so (pin takes 



the form: 



(z' + 1)3/2 _(z+l) 



3/2 



Step 3b: Next, we correct the error on by solving inhomogeneous 
Bessel equation. Although this introduces error on A, we solve the in- 
homogeneous Bessel for specific solutions that does not undo what we 
have achieved on A, ie a solution that satisfies the boundary asymp- 
totic at 7?. n 5^2 {dependent on A) and certain prescribed asymptotic 
behavior at infinity. As a result, the singularities of the solution lie on 
a different Lagrangian which we denote by F2, and not the incoming 
Lagrangian Fi. The phase function associated to F2 is — (^out where 



V'out 



(z' -t- 1)3/2 -I- (z+ 1)3/2 



• Step 3c: In the third step, we correct the resulting error from the above 
step on the semiclassical face A but this time along the reflected la- 
grangian F2 in phase space over the corner Ar\'S2- The initial condition 



is now d.i'R^^A, see figure 15 



(4) Step 4: We remove error at F. The normal operator on F is 

71 

Nm) = h'[{-^dzf + {^--r)\ 



On J- ^ we use Mellin transform to construct solution with wanted asymp- 
totic behavior at 7?. n and Cr\F. 
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Figure 15. Figure showing lagrangian r2 associated with face 
function ipouu the shaded region represents {z < 1,2;' < 1} 



(5) Step 5: By remark (9.2), we are now ready to remove error on *Bi and 

N<sAQ) = {idif + (A - — ) + P : Bessel type ODE 

%i(Q) = (idif + ^)-t^ ■ Bessel type ODE 

*Bi and 5Bi are blown- up versions of the product structure x M^, shown 
on the right of figure [TSj 

(6) Step 6: At the Zq > end, we remove error at at that is introduced 
there from the process of error removal at *Bi in the previous step. 

(7) Step 7: Finally, we remove error at TZ. 

In the following subsections, we expand in details the various steps of the construc- 
tion. 



9.1. Step 1 : Removing error at the lifted diagonal. 

In the first step we construct Uq € ^'^^ such that Eq — Q(/i, X)Uo{h, A) — 6{z — z') 
has no singularity at the lifted diagonal Diag^ with Uq supported in a neighbor- 
hood of Diagb ;j in fXUb^h that only intersects the boundaries oiDJlb^h at *Bi, *82, 5Bi, SS2 
and T. 

(23) QiK A) = h^zd,)^ + h^iX - — ) + z3 + 

This is elliptic in the normal bundle of the interior of the lift Diag^ ^ (of the di- 
agonal) to dJlb.h- We would like to carry this elliptic construction uniformly up 
to the b-front face J-, the blown-up front face 58i,*B2,?Bi and *B2. The diagonal 
Diagj, ^ meets these five boundary faces transversally, and hence can be extended to 
a neighborhood of J-, *8i, S2, *Bi and *B2 in the double space of VJlb^h- We want to 
extend Qb,h to this neighborhood (in the double space) and apply the standard the- 
ory of conormal distributions to an embedded submanifold without boundary (which 
in this case is the extension of the diagonal Diag;, f^). In order to do this, we need 
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to show that Q is transversally elliptic or is at least a multiple of a transversally 
elliptic operator down to these five boundary faces. 

• b-front face T : Denote the coordinate on the left factor of X by 2; 
while the coordinate on the right factor will be z' . Define the projective 
coordinates 

p = z' ,Z = — h 
z' 

which hold away from the left face. The b - front face is given hy F = {p = 
0} and the lift of the diagonal is Diag;, f^ = {Z = 1}. Since zd^ lifts to Zdz, 
the lift of Q(/i, A) under the b blow-down map (5^ is equal to 

(24) Qb{h, A) = r Q(/i, A) = h\Z(hf + h\\ - ^) + {Zpf + {Zpf 

The operator is transversally elliptic in a neighborhood of {Z = 1} away 
from ft, = 0. The restriction of the lift of Q(/i, A) to the front face T = Q = 
{p = 0} is given by 



N:F{Qlb{h,X)) = h' 



(Z5z)^ + (A--) 



• The blow-up front face ^Bi : now we blow-up the intersection of b front 
face and semiclassical face {ft = 0} in X"^. In coordinates 



(25) 



(26) 



p! = -,Z,ft 
h 



the operator takes the form 



-) + {Zp'hf + {Zp'hf 



{Zdzf + {\-'-) + h{Zp'f + {Z~p'f 



The operator in the bracket restricts to *Bi face (given by by {ft = 0}) is 
of the form: 



{Zdzf + (A 



l\2 



(Z/i') 



Hence (^b,h is a factor of an operator which is elliptic in a neighborhood of 
{Z = 1} down to ft = 0. 
• The blown-up face '^2 '■ the front face comes from blowing up the in- 
tersection of the lifted Diag;, rg and semiclassical face {v' = 0}. Define the 
projective coordinates by 

7 Z-1 , 

Zh = —^^v ,z 
The lift of Qb.ft, to $Xti in a neighborhood of S2: 

((1 + v'Z^)dzS + iy')\^ - ^) + {T-hv' + + {^hv' + 1)^ 

This operator when restricted to 0^2 (defined by \y' = 0}) gives the normal 
operator 

which is transversally elliptic down to 0^2 near = 0. 
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(27) 



(28) 



• The blow-up front face *Bi : In the neighborhood of this face we use 
coordinates 

z' 
h 

in which, the operator takes the form 

(Z9z)^ + (A-^) + MZMT-(Z/i')' 

The operator in the bracket restricts to SBi face, which is given by by 
{h = 0}, is : 

W)' + (A-^)-(Z/r)2 

Hence Q^^^ is a factor of an operator which is elliptic in a neighborhood of 
{Z = 1} down to /i = 0. 
► The blown- up face *82: The front face comes from the blow-up of the 
intersection of the lifted Diag^^ and semiclassical face A. We use the 
projective coordinates: 7.^ = 

,2 



^ ^ , , z' , in which the lift has the form: 



((1 + i>'Zh)dz,f + - ^) + {Zhv' + Ifp - {Zuv' + 1)^ 

/3,*Q(/),, A) restricted to *B2 ( defined by \y' = 0}) then gives the normal 
operator 

which is transversally elliptic down to 5B2 near z>' = 0. 



9.1.1. Properties of newly constructed function and the resulting error. 
From the verification above, one can find distribution Uo with the following prop- 
erties 

• The support of Uq is close to Diag;, and which intersects only the boundary 
faces J", «Bi, <B2, and ^2 

• Uo is cornormal of degree —2 to Diag^ ,j 

• Qb,hUo — Id = Eq, where supp Eq is close to Diag;, ^ and intersects only the 
boundary faces J^, Q5i and *82- Moreover, Eq is smooth across Diag;, and 
down to T {with values in conormal distribution to the lifted diagonal). 

9.2. Step 2 : remove the error at the semiclassical front face ^2 and ^2- 

We will find a distribution Ui, U-^ such that Ui is supported near *B2, and U-^ 
is supported near *B2 such that 

Qh{Uo + Ui + U,^^J = E, 
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where 

7<8i = 7si = -2; 7-4 = "1 
with /3^i^Fi supported away from £ and TZ and close to 052, whereas is 
supported away from *82, i3 and TZ. 

9.2.1. Preparation ingredients for ^2- 
Fiber structure near the front face *B2 

In the current neighborhood we are away from TZ and C The coordinates in use 

PA = 'Z-h^; p<s2 = ^' 

Expansion of a smooth function around ^82: If w G C°°{^b,h) expand v in 
Taylor series around S up to order N 

(29) ^p|,«^+z;',«^eC°°(Z,„z),«'ep^+iC°°(97lb,,) 

k<N 

Corresponding change to the operator in order to act on the expression 



of the form : Wc have Q lifts 
Clh^PlQ{h,\)=p'' 



((1 + v'Zh)dz,f + i^nx - ^) + (Z.i>' + Ifp + (Z,j>' + 1)2 



= AfmM) + P' [2z>'Z,ai + (i^rzidi + 9' {I + 9%,)dz, + (i>')'(A - y ) 

+ (3Z,i^' + 3Zl{Dr + Zl{mP + (2Z.:>' + Zl{u'f)_ 
where 

The resolvent for &i + p + 1 

Since the semiclassical normal operator is the Euclidean laplacian on (m = 

in our current case) with a family of spectral parameter 1 + p (in the case there is 
no oj) uniformly down to p — Q . 

Af<B.(Q/.) = 9z, + 1 
The asymptotic expansion of the outgoing solution for N{(^h)u = f where / is 
smooth of compact support 

^^^-.(p+i)V^|z.|^^_^0(|Z,|5) 
In addition, for g smooth in Zh and z', g of compact support in Z/j : 

Q(/i, A)(:>')'5 = (^O'^s. {Qh)g + i^'+'Qhg 
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where Q/j is a differential operator in dzf, witlr polynomial coefficients of degree < 3 
in Zf^. Denote by Rq the resolvent to Nrg^{Qii). We can apply Lemma 5.2 in the 
paper by S.B - M - V, to apply Rq to an expression of the form L^Rog, where L3 is a 
differential operator with polynomial coefficients of degree < 3. In our specific case 
with the dimension of the space is 1 and the coefhcient of the differential operators 
can be up to degree 2, we have 

Lemma 9.1. [Lemma 5.2 in S.B-M-V] Suppose Mj.j = 1,...,N are different 
operators with polynomial coefficients of degree 3 on M. Then 

R0M1R0M2 . . . RqMnRo ■■ C^iR) C-°°(R) 



has an analytic extension from 3(ct) < to 5Rcr > 0, 3cr G M so with a — + 1 
we have 

R0M1R0M2 . . . RoMnRo ■■ Cr(]Rx) ^ e-'^^l'^l \xf^ 
9.2.2. Construction for ^82. 

Now we can start the process of eliminating the error on the semiclassical front 
face, 



• We expand Eq as in ( |29[ ) 

note that each i?o.fc is compactly supported in Z/j 

• Firstly we construct 

Uifl = —p~'^RoEoQ 

From the Lemma 



9.1 



we have Uifi is of the asymptotic form : 



The exponent part e^''^''+^^^^^l^''l is the restriction to ft, = of 

exp ( - ^lii^'zr'\ii^%^z' + 1 + z')'^' - {z' + 
Denote by C/1,0 = -e~^'^'"+'(''+i)'''|Z''li?o£^o,o where 

^in = lii>'zr'\{i>'ZHz' + 1 + Z'f/^ iz' + 1)3/2 

We consider the difference between the phase function ipi^ and the pulled 
back of the e"''^''"''^^^^^!^''! using the product structure near the semiclassical 
front face ©2. 

The prep computations included at the end of the subsection gives: 

- [1 - e-^v,.+^(P+l)'^'\z.\^Eo,o - 2 (az,.e-'^'"+'(''+i)'''|Z''l) dz,RoEo,o 

- (a|^e-^^'»+'(''+i)'''l^''l) RoEo - i>'Q,,„p-2g-^^.„+»(p+i)V^|z,|^^^^^^ 

which means that we have replaced the error z''(£'')^^£'o.o by an error Eq i 
of the form zOO((i^')°). 
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Continue in this way to construct 

Then define C/i,i = e-*'^-+'(''+i)'^'|Z''IC/i,i so that Ui,i is equal to the 
restriction of Ui^i on 582- And in the same manner we construct Ui^k- 
Note the procedure of applying Rq to the new resulting error is assured by 



Lemma 9.1 in (X3i)- 

• Thus we have constructed : 

with (Pin = lii^'z'r' I (i^'^hP + 1 + p)'/' - (p + 1)3/' 

oo oo 

and : Q, - ^(:>')'~'£^o,fc 

where the statement means the equality of coefficients. 

• Borel summation to obtain Ui 



to get Ul, and we can assume Ui to be supported in a small neighborhood 
of psg^ is small. 



PREP COMPUTATIONS 

-p^Ns{Qh)Ui,o 



Away from the singular set of the phase function, compute: 



lim ({>') 



,-iVi„+«(p+l)'/"|Z^| 



1 



= - lim e-^^'"+''(''+i)'''l^''la,.(pin 

/i-i-O 

by several application of I'hopital's rule 1 



{i + zT'/'iz^pf 



-iiPin+iiP+iy^^\Zh\ ^ _ 



= - sgn(Z/i)e 



-l¥'i„+l(p+l)l/='|Zh| 



i(z>'Z,p+l+p)l/2p-*(p+ 1)1/2' 



52^g-«^.„+.KP+i)^/^|z.| ^ _ sgn(Z;,)e'^--'(''+i)'''|Z''l [z(i>'Z;,p + 1 + p^/^p - i[p + Xf^ 
- sgn(Z,)e^^">-'(''+i)'''|z-l*i(;>'Z,p + 1 + p)-'/^i>'p 



Since 



lim 

5'->-0 



i{i>'ZhP + 1 + p)i/V - i{p + 1)1/2 



38 



we have 

Zh 5 Z/i V / 

9.2.3. Construction for 582- 

Fiber structure near the front face *B2 

In the current neighborhood we are away from TZ and C. The coordinates in use 

7 Z-1 /l , 

The process is much simpler since M^^{(^h) is now an eUiptic operator 
Q(/i, X){t''fg = (;>')'A/'<B,(Q/0.9 + i>'+'Q-a.'7 

By Fourier transform we can construct the solution u G 5 (rapidly decay at infinity) 
of N.^^{Q,h)u = f where / is smooth of compact support. Denote by Rq the 
resolvent that produces such rapidly decaying solutions. 

We work with the portion of Eq close to S2, denoted by Eq ^g^. We expand Eq 
at «B2: 

note that each E^ ^^./^ is compactly supported in Z/j. Let 
Thus we have 

00 00 

fc=0 fe=0 

where the statement means the equality of coefficients. Borel summation to obtain 

hoTcl sum , ^,-,-lSr^ / ~r\-j rr 
3 

we can assume to be supported in a small neighborhood of is small. 



9.3. Step 3 : removing error at A and 5^2- 

Removing error at A and ^2 will be carried out in 3 smaller steps. Firstly, on ^Mb^h 
we remove the error at the semiclassical face A along the Lagrangian Fi associated 
to the phase function: 

|/»-l|(z' + 1)3/2 -(z + l)3/2| 

by solving the transport equations. See the following figure (14). 
Secondly, we correct the error on 5^2 by solving inhomogeneous Bessels for specific 
solutions that have the required asymptotic type at the 710^2 while resulting in 
an error that does not undo what we have achieved on A. Although this introduces 
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error on A, the resulting error will have singularities on a different Lagrangian 
which we denote by and not Fi. In order for this to happen, we require that 
the solution to the inhomogeneous Bessel equation satisfy a fixed type of oscillatory 
behavior at An ^2- In short, before correcting error on ^^2, singularities lie on the 
lagrangian Fi associated with the following phase function —ipin- 



(z' + 1)3/2 -(z + 1) 



3/2 



In the region close to Ar\^2 we have z < z' so Lpi-^ takes the form: 

After we have constructed the solution that satisfies the above stated requirements, 
the error resulting will have singularities along Lagrangian F2, which is associated 
with the following phase function —(pont- {see the remark at the beginning of step 
3b for reasons as to why if out should have the following form) 

= [(z' + 1)3/2 + 1)3/2 

In the third step, we correct the resulting error from the above step on the semi- 
classical face A but this time along the reflected lagrangian r2 in phase space over 



the corner Ad ^2- See figure 15 



9.3.1. Some prep computations. For matter of convenience we write out the form 
the operator in several coordinates in the neighborhood of : 



In i, v' , z' 



In v' = i^,z',Z : (z'z>')"(Z9z)" + {z'v')\\ - ^) + (Zz')' + (Zz')' 



(z')2 (t>')'(Zaz)2 + {v'f(\ - ^) + Z3z' + Z2 



(z')2 {v')\td,f + (y'f(\ - ^) + {iv'fz' + {iv' f 



{z'v'f (fa,~)2 + (A-^)+P^>V + P 



In this coordinates, *8i can locally be described as {z' — 0} where as A locally 
is \tr^ = 0}. On *Bi, the operator in the bracket is a Bessel equation of order 
(^-A)i/2: 

(%)2 + (A-^)+P 
The projective coordinates associated to the above blow-up ie close to 



Z 
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Note that away from TZ the flow i>' — constant looks like i/' 



h h 



z' zT Z' 

\x = constant. The forms of the ODE in various coordinates in this neighborhood: 



h\z'df, - zd,Y + h\\ - -^) + z^ + z' 



In Z', z, h : 

lni>=-,Z\z: {^zvfijld^, + vdv ~ zd,f + {zDf{\ - ^) + + 2^ 



z 

2 

In A = Z', z : {zZ'fifiZ'dz, - zd,f + {zZ'tif{\ - !L) + ^3 ^ ^2 
ln~y=%,u,z: {zDf {Ddp - zd,f + (zpf {X - + z^ + z"" 



9.3.2. Step 3a. We can assume that error has error close to A ie has supported in 
{z>' = < e}. We use the product structure close to A 



PA 



X A 



PA = ^'(Z) 

On DJlb,h the coordinates close to A and closer to the corner 5^2 H ^ : 

~^' = j{=i-'),Z,z' 
We compute the form of the operator Q;, when conjugated by 

exp ( - q{^'Zz')-^\ {Zz' + 1)3/2 _ ^ i^3/2|^ 

• on Z < 1 

2 

{^'Zz'f {Zdz - ^d~f + {l,'Zz'f{\ - ^) - 2^^'Zz'(Zz' + ifl^Zz' {Zdz - 

- i^'Zz'iZz' + ifl'^Zz' - i~v'Zz'^-{Zz' + l)-^l^{Zz'f 

2 

= i/{Zz'f [y' {Zdz - i'd^f + ^'(A - y) - 2i(Zz' + if'^ [Zdz - ^d--) 

- i{Zz' + 1)1/2 _ iZz'^{Zz' + l)-i/2 

• and on Z > 1 



9'{Zz'f {Zdz - ~yd^f + ~y'{X - ^) + + 1)'^' - ^^5) 
+ i{Zz' + lY'^ + iZz'\{Zz' + 1)-''^ 

To get the transport equations we conjugate the operator by exp(— it^in). Denote 
by 



Qb = P^(^') T^ZPfe,^ 



(Z) 
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Using the above computation of the conjugated form of Q(, we have: 
(30) 



,2 , 



9 



2i{Zz' + 1)1/2 (^Q_^ ^ ^(_z_^2 _J^^_ ^(2^, ^ -^)i/2 _ iZz'i(Zz' + 1) 



-1/ 



(31) 



n 



2 n 



+ (A-^) 



2i(Zz' + 1)1/2 (^d^^ki^f 



■z(Zz' + 1)1/2 + *Zz'^(Zz' + l)-i/2 



In coordinates = * Z' semiclassical face A is defined by {y' = 0}. 

£Ji can be considered as being of the form 

= exp(-i^i„)i^i; Fi e /9;^-^+iC°°($m6,fe); 7.1 = -1 

Consider the asymptotic expansion at the boundary face i/' = of the form 



-2; 7^ 



The construction of an approximate solution U2 so that 

U2 = exp{—i(pin)U2', U2 is cornormal and polyhomogeneous on A2 

Qb{h, A)e-'^'°C/2 - e-''^'"Fi = e-^'^'-Ra; R3 e p^Cf^,, 
On the region Z < 1, the integral curves hit the boundary TZnA hence we can solve 
the regular singular ODE and solve each transport equation for U2,j uniformly down 
to the Si n A. On the other hand, in the region Z > 1, the integral curves do not 
hit CnA and the fiber becomes singular. On this region, to find polyhomogeneous 
approximate solution we follows exactly the same methods as in the step 3c. In 
fact the transport equations in the two situations have the exactly same 
form; so we refer the readers to step 3c for the motivations and details. Here we 
will still write out the transport equations. After constructing U2J, U2 will then 
be constructed using Borel summation of these U2,j ■ 

(1) FIRST TRANSPORT EQUATION: comes from coefficient on p^^'^, 
where 7^ = —1. We solve for f/2,0 on Z < 1 as the unique solution satisfying 
the Cauchy problem : 



(32) 



Z<1; -2i(Z^' + l)i/2Z^(z')2 



Z5z + (7^ + 0)(^)2-(7^ + 0) + ^ + 



Zz' 



2 4(Zy + l) 

= — -Pi,o 
lC/4,o|z=i = -0: IC 

Since Fij are of compact support around Z = 1, the ODE is regular singu- 
lar, C/2,0 will behave like C/2,0 ~ 7^ — | at Z = 0. 
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(33) 



For region Z > 1 we use method of step 3c to obtain 1/2,0 (see details of 
approximate solution of the first transport equation in step 3c) 



Z>1; 2i(Z^' + l)V2Z — (2')' 

= —Fi,o 
lC/2,o|z=i =-0: IC 



Zz' 



2 4(Z^' + 1)_ 



For this region we obtain 1/2,0 polyhomogeneous conormal at the corner of 

where ag-^ = a<B^ + 5 + 1; a^Bi = +2-2; 7951 = -2 
(2) HIGHER TRANSPORT EQUATION : comes from coefficient on 
Pa 



For region Z < 1 



(34) 

-2i{Zz' + iy/^Z^iz'f 



Zdz + {lA + k)i^f-{7A + k) + l+ 



2 4{Zz' + l) 



U2,k 



Zaz + hA + k- f -ilA + k- 1)J + (A - -)J U2,k-i = -Fi,k 



The asymptotic behavior of U2,k at Z = which is Z'^^ 2 . Hence on the 
region Z < 1 



U2,k e {zT'^^pf+''Z-^-2C'^{Z) 



For Z > 1 



(35) 

2z(Zz' + l)V2zJ^(,')2 



Zdz 



'z + (7^ + fc)( <|y)^ - (7^ + /c) + - + ^^^^ 



2,;c 



^ ^'-(7^ + fc 



+ (Jy^')'[(zaz + (7^ + fc-l)(<z> 

using method in step 3c we construct 



1)) +(A 



n 



2,fe-l 



where a^^^ = a>8i + 5 + 1; "Si = 7!Bi + 2 - 2; 703^ 



Take borel summation to construct U2 

where = a^^ + 5 + 1; Q!<Bi = 7»i + 2 - 2; 7iB, = -2; 7^ = -1; 75^ = 7_4 - i 
Define 

(36) C/2 = exp(-i(^in)i/2 

we can take f/2 to be supported near A by taking small cut off in p_4 
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The error after this stage close to A : Define 

E2 = QU2 + El 

(37) ^ E^^eM-^^M; e p'S p^^™ pZ'l'^' pf^ 

where a<B^ = 7*81 + 2 - 2; 725^ = -2; •^^^ = 7,4 - i 

9.3.3. Step 3b : remove error on '^2- 
PREP INGREDIENTS 



• Product structure: We use the foUowing product structure in a neighbor- 
hood of ^2 

M+ ^ X ^2 ; ^2 = if X [0, 1], 

• A global defining function for can be chosen as 

In the neigborhood close to away from A (t) ^ 1 thus P' acts as the 
defining function in this neigborhood. Since i ~ j^, away from TZ 

p,, = Z(l + (f-i)^)V2 

thus Z acts as the defining function. 

Form of the ODE : The ODE in coordinates t = ^^i>\z' has the following 



form: 

Qti,,H = {z'fplSiY^ 



[tdif + (A - ^) + ij^P^.z' + P 



• PHASE FUNCTION: In the current region (fi^ has the form 

2/ [/ / I -in3/2 /7~i I I i\3/2 

V'in = ^V'V) [{Z + 1) ' - (il^ Z + 1) ' 



(z';>')"' [({z' + 1)'/' - 1) - ((fc'^' + 1)'/' - 1) 



1 



1 



One can see that the second term behaves like —t by writing out its 
taylor series {note we are working in tv' z' « 1) 



(38) 



|(zV)-i[l + iti>V + ...-l]=-f[...] 



In order to isolate this behavior, we will work with the conjugated 
version of Q by 



exp 



[(Z' + 1)3/2 _ ll ) 



however, since this commutes with Qb /i, the conjugation leaves 
the operator unchanged. 

As indicated at the beginning of step 3, in order not to undo the effect we 
have achieved so far in step 3a, we will construct solution whose error has 
singularities along lagrangian C2 and thus has to behave like exp(— zi) close 
to ^ n 032) which is opposite of the behavior ~ exp(jt) as shown in (38) 
associated to langrangian Ci. Thus in order to achieve this effect {note that 
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any choice of phase Junction has to satisfy the eikonal equation), the phase 
function associated the outgoing lagrangian C2 has to be — t/Jout '■ 



^out 



\3/2 



1 



1 



{Z'v'r^ {Z' + 1)3/2 ^ ^ ^) 



\3/2 



CONSTRUCTION 

Properties of the error we are working with : Write E2 — + where E2 is 
supported closer to ^i. We are working with the portion of E2 that is close to 5^2 
ie For matter of convenience, we will restate the properties (37) of E2 in terms 
of the coordinate in region close to ^2 ie in coordinates t, v' , and z' . 



• E^ is supported away from i = 

• EI vanishes rapidly as i — >■ cxd. 

Following to observation made on the choice of the phase function, to start step 3b, 
we first rewrite the error from step 3b as 



exp(-f(z'z>')"' (^' + 1) 



1 



Expand E2 around boundary face ^2 whose defining function is given by using the 
global face function p^^ and take into account the stated properties of E2 at the 
end of step 3a : 



E2 



1,^+1+2 ^^r^^^ ^+2 y k j^^^. ^j^g^g ^^^^ g C°°(z,5+(f)) 



k 



where C°° {z' , (t)) comprises of functions g that are smooth in z and t, is sup- 
ported away from t — and rapidly decreases as t — >■ 00, together with suppi? G 



We solve the following ODEs for C/3 ^ : 



(%)^ + (A-^) + P^V + P 



"3-2 =73-2+1 

Solving in series for (by matching coefficients of pjj) ^3,fc amounts to conjugating 
the operator for (see step 5 for this approach). However the operator if conju- 
gated only powers of has a simpler form and in fact gives a simple operator of 
the type Bessel, thus we rewrite in the above expression with pg^ — i'' (t): 



(39) 



{idif + (A 



, 72-1 I I 72 
+ t V Z + t 



(j>')"S2(;z')"'Bl ^(l>')'=(iy+"52C/, 



3,fc 



(i>')"S2 (/)"'Bi ^(,;')fc(ty+"S2+2£:2^ 
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Solving in series by equating the coefficients of v' , we start with the foUowing 
inhomogeneous Bessel of order a = A + ^ 



n 



(tdiY + {x--) + t 



{t)"^^{zT'^-U3,o = {zT'^^{t)"^^{tfE2,o 



for which we construct a solution that satisfies : 

At t ~ : {z')°'^iU3fi e i"'{z')°'^^C°°{i, z'); a = ^J-X+^ 

As i ^- oo : i^3_o g exp(-if)f- V2(^')a!Bi 50. 

where comprises of smooth functions behaving as a symbol of order at infinity 

by choosing 



(i) 

where iji* ^ being the Hankcl functions (they are two linearly independent solutions 

to the homomgeneous Bessel equation of order a — \J —\ + \) and W being the 
Wronskian: 



Ai 
'wt 



c is chosen so that t^Usfllf-o ~ ^ (This is possible, since 

i~"C°°(t)). Since the second terms rapidly decreasing as f ^ 00, 11^^ satisfies the 
required oscillatory behavior . 



Define 
(40) 



then by the prep calculation placed at the end of the subsection we have f/3,o|z'=o = 
U^^Q. In addition we also have : 



{tdif + {X-'-)+Pv'z' + t' 



(i>')"52 {i)°'^-^{z'y^Uifi 



Wc next solve for C/3 for /c > in a similar manner, ie solve for t/3 ^ 
(41) 



n 



{tdiY + {\--)+t'v'z' + e 



{tf+''^-{z'r'^^us,k 

Properties of U^^k 



t = : 

t ~ 00 : 



Define 
(42) 



U3,k e exp(-ii)f-i/2(f)-('=^+"ff2)c°°(f) 
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Finally we construct t/3 by borel sum in p:^. 

borel sum in a 



(43) 



C/3 



Multiply back by the factor \{z' + l)''/^ — l] that we have implicitly fac- 

tored out at the beginning of the construction (view the remarks in the prep ingre- 
dients subsection), we obtain 

Us = exp ( - i^-{z'v')-^ [{z' + 1)3/2 _ ^[^^ ^ exp(-i^o„t)C>3 

(44) Vont = liz'i'T^ \{z' + 1)3/2 ^ ^ ^^3/2 _ 2^ 



e Pgj^ p-^ C (iUCb,/,) 

Prep computations used in the construction 

2 (tv'z' + 1)3/2 _ 1 _ 

lim ^ = t 

o'^o 3 ly'z' 

, 2 (fz>V + 1)3/2-1 , ^ ^ 
lim cxp(-^^ t) -1 = 



idt 



2 1 
3^ 
2 1 

3^ 



4,((ff>V + l)3/2-l)-t 

4,((f^V + l)3/2 



I -t 



1 ^2^;/^/ 



= i[(tW + 1)1/2-1 
= f[(tW + 1)1/2-1 

_ Pi)'z' 

~~ 1 + {ii>'z' + 1)1/2 + 2 + 1)1/2 



2 {ii>'z' + 1)1/2 



[(iaj)2 + (A - ^) + P + Pi>'z'] {tr^- {zT""' Usfl 

2 

= [(%)2 + (A - ^) + P + PiJ'^'] e-**+" (t)"s2 (^')"^i i73_o 

=e-*^+^* [(t5f)2 + (A - ^) + (/)"^i 

+ [(ta4)2, e-''^+**](i)"»2 (^')«!Bi f/3 Q + pjp/_j/^^ys2 [73 

Properties of constructed distrubtion and the resulting error Denote by 

-£'3 = QbUs + E2 



'^aep^tV^^"" p2+".v;:^C°°(9Jl,,,) 



C/3 = exp(-iv3out)C/3 



47 



9.3.4. STEP SC. 



Step 3b introduces error back on A however this time along a different Lagrangian 
associated to the (^outi which we will proceed to remove in this step. 
The product structure We use the product structure close to A 

The 'global' defining function for face A 

PA = ^'(Z) 

Reminder of the phase functions we are using: 

2, 



(z' + 1)3/2 + 1)3/2 



In the current region z' > z, ipi^ = -h ^ {z' + 1)^/2 - {z + 1)^/2 

3 



Vo..t= ^h-^ [(Z' + 1)3/2 + 1)3/2 „ 2 



Prep computations: 



Zdz^out = (Zz' + l)i/2Zz' 
(Z9z)Vont = (Zz' + lY'^Zz' + \{Zz' + l)-i/2(Zz')' 



The conjugated operator by the outgoing oscillatory part. To get the 

transport equations to solve for the 'symbols' we conjugate the operator by 



exp(-i(^out) = exp ( - iliD'Zz'r^ {{Zz' + 1)^/2 + (/ + 1)3/2 _ 3] ) 



(45) 



(^'Zz')' (Z9z - ~vdi) + (^'Zz')'(A - y) - 2z^'Zz'(Zz' + l)i/2Zz' (Z9z - 

- z^'Zz'(Zz' + l)i/2Zz' - i^'Zz'^(Zz' + l)-i/2(Zz')' 
^'(Zz')' (Zaz - ~vdi)^ + ^'(A - ^) - 2^(Zz' + 1)^/2 (Z^z - 

- ^(Zz' + 1)1/2 _ iZz'^(Zz' + l)-i/2 



Solving in series for the coefficients of pA amounts to further conjugation of the 
above operator (451 by powers of p^. 
Denote by 
(46) 



Pb,„ - ~v' {Zdz - ^9=)'+^'(A-^)-2i(Zz'+l)i/2 (Zaz - ^a=)-i(Zz'+l)i/2-zZz'i(Zz'+l)- 



1/2 



(47) 



Pa 

^P^ 



Zdz + k{t^f-k 



(A 



- 2i(Zz' + 1)1/2 ^2^2 + fc(^)2 -kj- i{Zz' + 1)1/2 _ iZz'^(Zz' + l)-i/2 
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The error from step 3b : Write = + El^ with E^^Ef- supported closer 
to A and TZ correspondingly. For this step, we will only need to work with E^. 
Write 

E^ = exp(-^^out)i^3; ^3 e p7:^'p^pt"''''%fAi^',i) 
Expand F3 at yl. 

taylor expand 5+QS2+2 a«i+2 v;^ fe p . p c r°° f A\ 
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From the above computation ( 47 1 we can write the transport equations satisfied by 
the 'symbols' [/4_fc : 

(1) FIRST TRANSPORT EQUATION comes from coefficient on p]^'^^: 

(48) 



'-2z(Zz' + 1)1/2 [z9z + (/3^ + 0)(^)2-(/3^ + 0) + J+ 



U) 



2 4(Zz' + 1) 
C/4 o|z=o = ~0 ■ IC double check the IC 

We have denoted by Pa the leading power oi pA, k Pa — \ + ^" 



•^4.0 



Denote by Pq the operator 



l\2 



Po = -2z(Zz' + 1)1/2 



Z .0 1 Zz' 



Zaz + (/3^ + 0)( — )^ _ (/3_4 + 0) + - + 



Z)' v^--- ' / ■ 2 4(Zz' + l) 



PU2o + -F3.0 = /o which is nontrivial at *Bi and C. 

Remark 9.3. | = (z2+(z^')2)i/2 ■ Hence Pq does not live on the blow down 
of A, which is the product space Rf x M J, . 

The short coming of just solving along each regular fiber: If one 

were to carry out solving just along each regular fiber, we will get a solution 
that does not live on the product space M+z x [0, l\z' ; ( the solution lives 
on ^ [0, l]z' but does not have uniform behavior as Z ^ 00 and z' — > Oj. 
When considered on A, we have the existence of the solution constructed 
this way on every compact subset of A away from A \ {^i U *Bi) and we do 
not know the asymptotic behavior of the solution at the corner *Bi H in 
A and along 



New method to construct polyhomogeneous solution to (48): Be- 
cause of the singular fiber (03 1 n ^) U {"Si r\ A), it will involve more work 
than just solving along each fiber. What we set out to do is : (the following 
procedure is completely in face A ) 

(a) Find a polyhomongeous function Uf]^ + on A so that 

/2 = Po(C^|o + Uf,o) ^ PcnAPAn<Si 

(we will further break this into two steps : first remove the error Q3i 
then at C then) 

(b) Now we can blowdown Q3i n .4 in ^ and consider f^ to live on the 
compactified product space [0, 1]^ x [0, l]z'. 



DETAILS OF THE CONSTRUCTION 
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• Step la : We will use as the 'global' defining function for *8i in A 
P'BinA = z'{Z) 

Since we are working close to 05 1 we will only need to consider Fs^o 
multiplied with a cut-off close *Bi (in A), denoted by -^30- Write 
Fa^o = ^3,0 + -^3,0- Expand at Bid A 

3 

Denote by ^«Bi the leading power of pos^ in the constructed solution 

(74,0. 

Wc want to solve for 

- PAnm, E PAnm.K^j such that PUf^^ G pg, ; /3s, = + 2 - 2 
3 

Solving in scries amounts to conjugating Pq by powers of p^^-^^nA which 
is the same as conjugating by powers of (Z) since power of z' commutes 
with the operator. By direct computation we get : 

P^'n^'^PopLtn!,^ = -2»(|yP<B.+l)^/'|^(Z)-VI,[zaz + (/33,,+j)(Jy)' 
+ {pA + Q){^f-{pA+^) + \ + \-—^ 



■(z> 

Since is a bounded smooth global function and have value less than 



1, we can expand in series = 1 — 7zYP!8in.4 + • • • Thus 

(-^P<Bin^+l) ^ ' 

U^Q.Q and U^Q.j have to satisfy 

(49) 



2i 
2i 



Zdz + (0 + J( Jy)' + m + 0)(|y)^ - (/3^ + 0) + ^ 
Zdz + {j + J(|y)^ + {f3A + 0)(|y)^ - (^^ + 0) + ^ 

4 



rrSSi _ 7-2 pi 



^i,3 



+ - ''^^CjU\i = — ^^CjZ~^Fg^Q J +; c; known constants 



;=o 1=0 

Use Borel sum to define from 

tt'Si/^ borel sum /3a5j V^/„ \7rr'Si (-7\ 

t^4,o(P^n<Bi,Z) ~ P^nsi Z^(^-4nQSirC/40y(Z) 

3 

Denote /i = Pof/fo^ + ^3^0 

Behaviors of constructed function and resulting error : The 
behavior of ^3,0;^ at Z = (ie at ^2) dictates that of U^J . Note F^fi-j 
has support in Z < e' for some e' and vanishes to infinite order at ^2- 
Thus 

/i ^ PAnmiPd2pdi 
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• Step lb : Now we remove the error at in A. Write /i = fl + f1 
where fl is supported away from ^2- Expand /i at S^i fl A. 

We would like to find 
^Ih ~ pl!'J2p?rUl^ so that PoUlh + fl e P^nj.PS'ns, 



In this case j3^^ = j3sQ^ + \ + 1 
In coordinates 2;, Z' , Pq is equal to 



Z' 



zd, - Z'dz, + + 0)(Z')-' - + 0) + - + 



2 4(^ + 1) 



Prep computation 



1 



Po(Z')™5(^) =zd,-m + i^A + 0)(Z')-' - (/3^ + 0) + - + 



2i(0+l)V2^2 



2 ■ 4(2 + 1) 

so C/^o;0 ^^"^ ^ifiipy have to satisfy correspondingly : 

1 2 



zd, - iP^, + 0) + (^^ + 0) - {Pa + 0) + 



4(z + l) 



'^4,0;0 



-/i\o(^) 



2i(2 + 1) ^g^ _ (^^^ ^ ^) ^ + 0) - (^^ + 0) + ^ + ^^^^ 

i-1 

+2i(2 + i)V2^2^dp;7|,o;p = -A\,. 

Use Borel sum to define U^ q.q 



'-^4,0;m 



j-T-^^ borel sum 
^4,0;0 ~ 



Denote ^ = PoC/|;);o + /i 

Properties of constructed function and resulting error : Since 
we are dealing with regular singular ODE, we can choose solution 
^io-m whose behavior 2 = is determined by that of which 
vanishes to infinite order in n>t at 2; = 0. Also we can assume U^]n 
are supported away from 5^2- 

U^Q e (0^^ p ^ ; supported away from ^2 

h G 1 PyJ ; supported away from ^2 

Prom the properties inherited by /2, we can consider f2 to live on the 

compactificd product space [0, l]z x [0, l]z'- In addition /2 G 
vanishes rapidly as 2' — >■ uniformly as 2 — >■ since we are away from 
Z = Oand {z' {Z))°° {Z-^)°° . 
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• Step 2 : Since wc have been able to blow down and consider the error 
on the product space, now we solve for Pg for U^q along each regular 
fiber (parametrized by z'). The ODE in z, z' is 



2i(^ + l)V2 



(z2 + (z')2)l/2 1 



2i{z + l) 



(z2 + (2')2)l/2 2 4{Z+1). 

the second coefficient is not quite 'nice'. However this is the conjugated 
version of a nicer operator, ie we will instead work with: 



1/2 _ 



(^2 + (^')2)l/2 



zdz 



1 



2 4(^ + 1) 



(l + (-)^) 2 t/|=„ 



Ha 



F|q is supported away from *Bi and has the following property at : 
-Ff Q e /2 also vanishes to infinite order at z = 0. Same reasoning 
as for /2, wc have /f supported away from and vanishes to infinite 
order at and 5^2 hence on the space [0, 1]^ x [0, l]^' it is in z°°{z')°° . 
In short, the RHS has the property : z°°{z')°° . Thus we can solve for 
a unique (1 + { — )'^)^^^'^U^o solution which satisfies 



Conclusion : We have constructed a polyhomogeneous function on A 



-'4.0 



^^4.0 



with PoUifi = -^3,0 

The behavior of C/4,0 on A is determined by the first two terms, more 
specifically: 

^ (74,0 e P^^p<s^ P^^ 
where fi^^ = l3^^ + i + 1; l3<s.i = q;<8i +2-2; a<s^ = 7<8i +2-2; = -2 
(2) HIGHER TRANSPORT EQUATIONS : come from coefHcients of 



Pa 



(50) 

2i(Z^' + l)V2i^ 



/\2 



Zz' 



2 4(Z;j' + l) 



4,fe 



Z2 r / \ 2 7i2 -, 

+ :^(^')'[ (zaz + {pA + k- -ipA + k- 1)) + (A - -)\u^,k-i = -Fs,k 

Using the exactly same method we can construct a polyhomogeneous conor- 
mal solution Ui^k-i on A with the property 

U4,k-1 e P5,P<Bi %i 

where P^g^ = P<3^ + 5 + 1; = C(<3i +2-2; agsi = 7<Bi +2-2; 713^ = -2 
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Now we can construct U4 on VJlb.h by taking the Borel sum 

fy borel sum B, \kTT ^ tt ^ Pa oo ^SBi /^Si 

k 

= 1 + aj, + 1 

Next define 

Ui = exp(i(^out)t^4; -^4 = Q,Ui + E3 
^E^ = exp(i<pout)F4; ^4 e pS'pr.plT^'pl? 

where = P<3^ + 5 + 1; /3<Bi = a<3^ +2-2; agsi = 7®i +2-2; 7,8^ = -2 

9.3.5. Summary of properties of resulting error close to A from step 3a 
- 3c. 

The error produced from the above steps : 



^3 ^ P^i P-R PAPS2^oi^'*) 

where a^Bi = 7>Si +2-2; 705^ = -2 

Remark 9.4. Note that all of the above-mentioned errors vanish to infinite order 
at A, hence can be considered to live on the blown down space, ie their portion 
close to *Bi can be consider to live on the product space 

where Bi is the blowndown version of !Bi where we blow-down ^ n Q5i in *Bi. In 

fact we will proceed to remove the error at b- front face to be able to consider 
the resulting error to live on Bi = x M.^,. 

9.4. Step 4 : remove the error at the b front face 

We will only work with the portion of the error supported close the J^, which in 
this coordinate in supported away from Z = and =0. This error comes from 
after step 1. Write the error as 

El = Ef^^ + Ef^^ + 

where the superscripts indicate the neighborhood of support. Each of the terms is 
dealt with in the same way. We write out the details in dealing with -Bi,<Bi i-^ the 
term which is supported close to Q5i. 

Product structure near T : 

where pjr is a choice of the global defining function for b blown-up front face T 
{global in the sense the function is valid at both TZ and C end). For a choice of 
global defining function, we use a function such that in neighborhood away from 
£, pjr has the form p^r = jl'{Z) and in coordinates valid away from right boundary 
n, p^ = il{Z'). 
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A choice of coordinate away from £ : 



/i'(=^),Z,/i 



while away from TZ we use 



j,{='-),Z' = Z-\h 
The form of the operator the coordinates in these coordinates 

/i^ [(Z5zf + (A - ^) + MZA'f + (Z/i')'] 

Prep Computation : Since any power of jl' commutes with Qb, conjugation Q;, 
by fP-p amounts to conjuation of by (Z) : 

= h'{z)-^ [{zdzf + (A - ^) + Hzn'f + {za'f] {zy 

= h' [{Zdzf + 2j(Z)-2z2(Z5z) + 2jZ2(Z)-2 + j{j - 2){Z)-^Z^ + (A - ^) + h{Zi,'f + [Zp' f] 



= e 



{Zchf + 2j-( Jy)'(Zaz) + 2j( - 2)( J^)^ + (A - ^) + H^^fp^ + {^^^p^ 



Note that ^ is a global function on each compactified fiber of 

Construction: We looked for U^'^ of the form U^^ P%J2j Pj^^^j^ ^^'^ ^^P" 
ported close to such that 



+ ^ e p^p%,pI " Pi " C~(9Jtb,ft) 

Expand E^^ at the b front face J^: 

Ef^ taylor^xpand ^ ^^^^^^ . ^ ^q^oo 

3 

Solving in series and using the prep computation, wc obtain the ODE-s which 
need to satisfy. Note that these U"^- are global solutions living on each fiber 
{h = constant} on T as follows : 
(51) 

(52) 



2 



(Zaz)^+2(7^+0)(|^)'(Zaz)+2(7^+0)(|^)%(7^+0)(7^+0-2)(|^)V(A-'^)Jf/5%^=-£i%^ 



H 



{Zdzf + 2(7^ + j)(|y)'(Z9z) + 2(7^ + j)(|y)' + (7^ + 3){it + J - 2)( Jy)' + (A - ^) 
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However for simpler computation, instead of solving directly for U^j we will look 
solve for {Z)'^^'^Wf^ which in turn satisfy simpler ODE in Z. 



(Z5zf + (A--) (Z 



5,0 



1,0 



(zaz)'+(A- 



5,J 



'5,3-2 



• We will revoke the facts about Mellin transform from remark (9.51 (at the 

end of the section) to obtain specific global solution for {Z)'^^^'^U^q that 
solves 



(Z5z)2 + (A - {Zy-+^U^^ = -{Zy'-+^ET^; 7^ = 



0: {zy-+"u^i 



{zr-+'u^^ = z 



--A 
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In this case we do not have poles contribution from the inhomongeous term 
since Ei^q has trivial behavior at both ends. 

• {Z)'^^+'^Ufl satisfies the same equation with (Z)t^+°,Ei,o switched to {Z)^^+^Ei^i 

• To deal with higher order equation for C/s.j, we consider the ODE that U^^2 
which has to satisfy. 



(A - {ZV-+^Ufi = -{Zr^+'Ef^i - h^z^zr^+^-^uf^^ 



{zy-+'u^^ = z 



Also use remark (9.51 but now keep in mind that the right hand side has 
nontrivial behavior at both ends. With the same notation used in the 
remark (9.5) we have /3 = sja + 2 while 7 = ^/a — 2 so using the Mellin 
transform we can get a global solution 

(Z) ■+2c/5»2i ^ Z^ + Z^+2 ^ Z^ 



5,2 



/^+2 ^ 2^-y^+2 



More generally for C/sj- 



/i2 
Z ' 



\zdzf + {\-'^)]v+^^u^] = -{zy^+^Ei^, - h^z^z)-'^+^-^uli_3 - /i2z2(z)^^+^-2c/5»/_2 



Zi+^^U^I = Z 



Z-^00: z^+^^u^i 
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With the same notation used in the remark (9.5), for this case, /3 ~ y^ + 2 
while 7 = y/a~ j so using the Mellin transform we can get a global solution 
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U^^ with 

Z-0 ; {Zy^+^U^I h-^Z^C'^iJ') 

z^oo : {zy^+w^j ~ h-''^z'^+^c°^{r) 

Now take 

r23i borel sum 



^23, bore^sum Y^^p^yuf^^ 



and thus f/^^ inherits the asymptotic properties : 



At 7^ n J" ie where Z = : Uf^ & h'^Z^ 4 •'c"" (p^, Z, /i) 



1 e /i-^Z ^C°°(p^,Z,/i) 



At £ n J" ie where Z"^ = : [/g® 



With superscript indicating the support, the error resulting from this step has the 
properties 



At 7e n J" ie where Z = Q:Ef^€ /i°z/^^+^C°°(p^, Z, h) 



At £ n J" ie where Z"^ = : i;e''i e /iV V 4 ^ 'C°°(p^, Z, ft,) 



Similarly, to take care of Ff^^ we use pj^ = /i'(Z) and in coordinates vahd away 

from right boundary TZ. pj- = /i(Z'). 

. ^ z' 

A choice of coordinate away from C is — ),Z,h while away from TZ we use 

h 

jl {= —),Z' = Z~^,h. The form of the operator the coordinates in these coordinates 

h'[izdzf+{x-'^)+h{zii'f-m'] 

As for ^inte""--, we use coordinates . (Note S = Z - Zq) 

Z = 5', Zo 

the defining function for is = S'{Z) while the operator has the form 

{Zdzf + A - ^ + (S'Z)^ + ZoiZS'f 
By the same method, we construct C/®^ and [/interior_ g^g below for their properties. 

9.4.1. Summarizing the properties of newly constructed function and the 
resulting error. With superscript indicating the support 

Ee = Ef' + Ef' = Q{U^' + Uf' + u'^'^''°') 
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(54) 



77'2?i ^ oo u V V 4 /^oo /cm \ 



f'ep^pl^pl ^ * C-(»l,^,) 



f^mtorior ^ ^T^^^ 4 (OTfc,,,); supp C/int<="°'- only intersects J", £ 



Remark 9.5. [on solving regular singular ODE using Mellin transform] 

Firstly we state the definition and some properties of the Mellin transform. For 
g £ C^(M+), the Mellin transform of g is defined as 

Mg{s)^ / ff(x)x-- 

s^Res + ilms; M : {R+ ; dx) -°-°J^'"-" ({1^ s = 6-^};dRes) 

M{xd,g) = -isMg 

The inverse Mellin transform is given by integrating along a horizontal line, provided 
this integral exists 

s = Res + iIms; M-i(g(.s)) = [ g{s)x-'' dRes 

J Im s— constant 

So in particular 

M"^ : L2({Ims = 0};dRes) ^ x^ ; dx) 

For our purpose, we will apply the Mellin transform to solve the following ODE on 
the half line [0,oo)x with / e 5[o,oo)i where 5[o.oo) has support away from and 
vanishes rapidly at infinity. 

(55) [(x^x)^ - a] u /; a > 

Take the Mellin transform on both sides 

{{~isf - a)Mu = Mf => Mu = ^5-^ 



and we obtain a specific solution for the ODE ( 55 1 



Ims=0 



Note that the RHS has poles at ±iy/a, a > 0. In order to know what the behavior 
of M is at X = and x — > oo, we deform the contour of integration. In this case since 
we assume / is 'nice', only the poles ±Y^|a| will affect the asymptotic behavior of 
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u. In greater details, to read off the asyraptotic property when x — >• 0, deform the 
contour of integration above Im s = 0, ie denote by u by 



Im s— c>0 



+ a 



then at X ^ 0, u ~ x'^ where c > since x = x ^Res^ims Also we have 

_ - , R I f x-^'i-Mf) ^ x--«(-M/) _ x-'^(-M/) 

For c > 

u = u + i-^x^{Mf){i^/a) 

M ~ x'^, c > 0. Since the above is true for all c > 0, take c arbitrarily large so we 
can make u vanishes to arbitrary order at x = 0. Hence at x ~ the asymptotic 
behavior of u is determined by the second term ie 

At x ~ : u x^ 

Similarly, to determine the behavior of u as x — )• 00, we deform the contour to 
arbitrarily negative end of the imaginary axis ie denote by u 

\S + a / Jlms=-c,c>0 s +a 

u ~ x~'=,c > . Since the above is true for all c > 0, take c arbitrarily large ie we 
can make u vanishes to arbitrary order at x~^ = 0. Hence at x ^ 00 the asymptotic 
behavior of u is determined by the second term ie 

At X ^ 00 : u x~^ 

Thus we have found a specific solution that solves 

f ((x^x)^ - a^u = /; a > 0, / e 5[o,oo) 



(56) 



: u 



This presents another v^^ay to view {xd^Y ~ which as a hyperbolic problem is 
a second order ODE. Instead we can view it as an global elliptic PDE and to obtain 
a specific solution, wo impose boundary condition, which in this case comprises of 
conditions at {Z = 0} and at {Z~^ = 0} {in fact, due to the kernel being symmetric 
we will impose the same Dirichlet boundary at both ends). 

In the case where the inhomogeneous part / has asymptotic at x ~ and 
x~^ . These exponent will contribute to the poles. In the case there are just 
simple poles then we use the same exact technique to read off the behavior of the u 
at X ~ and x~^ ~ 0, but now the difference of u and u will be a sum of residues. 
Heuristic computation to see the contribution of the inhomogeneous part to the 
poles : if / = /i + f2 where /i ^ with f3 7^ i^/oi and fi decays rapidly at 00, 
while /2 ~ x"''' with 7 7^ ±^/a and /i vanishes to infinite order at x = 0, 7 > 

Mf ~ / x^x''-'^ dx + / x-'^x'^-i dx ~ + 

Jo Jc P + ^s -7 + zs 
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For Mu = -T. , the RHS will now have three poles : iy^, 17. The poles 

s-^ + a 

contributing to the behavior of w at x ~ : iy^,if3 ie u ^ + . On the 
other hand, the poles contributing to the behavior of w at x^^ ~ : —i\Ja, — ^7 ie 

■U-X^v^ + X-T . 



9.5. Step 5 : remove the error at *Bi and Q3i. 



9.5.1. Prep computation. 



The projective coordinates away from A are 



t^-,i)',z' 



Defining function for *Bi : 

where g<Si has the property 

• Close d2 then g<s^ € C°°{mb,h)- 

• Close to S^i and C : g<s^ € ZC°°{mb,h) ; Z = ^. 

• Close to TZ and T : g<3^ e i>'{Mb,h) ^' ^ jr 

• Close to C and F : g<Si € Zi>(9Jlf,ji) ; Z = ^; £> 

In these coordinates i, v' and z' Q;, has the form : 



'^2 



(^') 



'.~/^2 



the operator in the bracket on *8i is a Bessel equation of order (- 



A) 



1/2 



9.5.2. CONSTR UCTION. The error results from step 3, step 4 and step 5, which 
we will denote by iJtg^ . 

7^ rpTZ I TT'Si 



(57) e p7:^'pI " '^^'p'XptAiz'rt) 



Support: E^^ is supported away from A, while is supported near TZ and away 
from 5i and y^. Since i^tg^ G p5'p^/9^C°°(9Jlf,^/i), we can consider iStg^ to live on 
Bi, the blown-down space from *Bi obtained by blowing down y^n *Bi and J"n Q3i 
( where Bi is the 2-dimensional product structure represented in figure (13 1 by 
the rectangle). See figure (13). 
Expand E^^ at *Bi 

taylor expand 7,1, j 77 ~ ■ / 1 o n\ n 

-Esi ^ PiBi 2^ PsSi^^Siy ; 7551 = mm(a<Bi + 2, 0) = 
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Properties of -E-Si j on *Bi from (57): they vanish to infinite order at aU boundary 
faces except at £ n *Bi and TZd^i where they have exponent — A. 
We solve for U?^ of the form 



U. 



However for simpler computation, in stead of conjugating by power of psii we con- 
jugate by power of which commutes with the operator. 



[/e.o has to satisfy: 

On the product structure Bi ~ x Mj, , properties of g^*^ ^(i>')~^E5,o at the 
more interesting corners: 

• Away from J-,TZ and 5^2, ffSi Z = ti)'. However since Es^o vanishes 
rapidly at i, in the corner of i^^ = and i>~^ — 5^*^^ (i>')~^E5.o G 



V 4 -''+f'Bi "where 5f denotes functions vanishing rapidly in t ^ = 
and f = 0. 

Away from di,A,d2,T^, g ^ Zv = D' . In the corner of i = 0, {v')^^ — 0, 



along each flow we can use the Hankel functions {solutions to the Bessel equation) 
to obtain the specific solution to the inhomogeneous equation. 

Furthermore, since the inhomogeneous part behaves in a specific way (using the first 

two listed properties of E<Sj^-j) we can uniquely solve for '^Uefi{i,h') 
with the following properties: 



. For A-^^ ^ Z : att ~ 0,gll^-'Ue,o ™*°'- X-^(^')-V X-^+^-i 



. As i-^ 00, "^"%xp-^*i-i/2(j>')-Vx-^+>^i-250(^~)^,. ^l^g^g 

S'^{t)i}i is symbol in t with smooth parameter in v' . 

Let X be a smooth compactly supported function with suppx(t) C {|t| < 1} and 
X{t) = 1 for |t| < i. Define Uq^q by 

(58) f/6,0 = (1 - X(i^>')) exp''^""+4 U2,^ + x{i^')U2.^o 

ie we modify t/2.0 (in fact the pull back of [/2.0 from Bi by the projection tt : 
Bi X IRpg^ ^ Si) on the region {tv' > 1} to have oscillatory part given by the 
phase function we have chosen ie i^i„. 
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Consider the resulting error at this stage 



[{idif + (A - ^) + ? + Pu'z']pl^^^-'Ue,o 
= [iid,r + (A - ^) + P + Pi>'z'] p^-/ ((1 - xiii>'))e'^-''^'^Ue,o + x{ii>')Ue,o) 

= (1 - x^Ke'^''^'^ - 1) + - + tlpl7~"Ue,o 

+ [^,-)^x(^^')]pS^"'c>6,o + [(ia,-)^ + (A - '}^)+P]pl-^-'u,^, + Pi>'z'pl^^^-'Ue,o 
And from the prep calculations placed at the end of the subsubsection, we have : 



>1F 0/'^^®i-2+^ 



Equations satisfied by U2J for j > : we are dealing with the same type of 
inhomongeous ODE: 



(59) 



{tdir + ix--)+t 
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„T"Bl-2+Jrr 



T-si-a+j, 



=5J 



- lim 



7SBi-2+i-l^ ^ 7iBj-2+i-l-. 



Since the the inhomogeneous part has the asymptotic we want {by the way we 

constructed U2.0 and then by induction U2,j-i inherits the same property), IJ2J 
also satisfies the following properties: 



For2J^^-A^Z: att~0,54^^ 



7!Bj-2+j asymptotic ~J ^ -X , ,.-J ^ -X+j^^-2+j 



00, ,gS^"'+'t/6j '^^^^ exp-*f-i/2(^/)-Vx-^+^-i-2+J50(^~).,. 



• As f 

where 5"° is symbol in t with smooth parameter in u'. 



Then wo define C/gj 

(60) (76 J = (1 - x(t'>')) cxp(-i^ii„, + It - tjrWej + x{ii>')U6,j 

where in the local coordinates on the region ti>' > 1 

thus the t/e.j has only oscillatory behavior at g^i n*8i. Take asymptotic summation 
to construct Ue: 



03 borel summation 7*33 j^— 2 



so C/f 1 = eM-mim)V^' + V^' 



where V,^^ G p:^pf^p^ p^^'' p^pl ^ "'p^^ 



ri2 1 
4 ^-^ 7-81-2+2 



^6 ^ PaP^^Pj^P^Bi Pc Pn PSi 
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7Q3i-2 V 4 



Remark 9.6. To determine the asymptotics of Vg^^^ and : 

(1) At the corner of 'Si and £, the local coordinates are 

PC ~ i>; Psi ^ i>;^i = z 

Away from J^.TZ and ^2, 9<Si ^ Z = u^^i'^^. Also expressing [>' and t in 
terms of the local coordinates: t — v^^, v' — v^^. Thus 

1 [r? T - , „ 1 rn? , 

(2) At the corner of C and the coordinates : Z~^,h,ji ~ ^. Away from 
i?i)-4iT?2,'7^, 5 ~ Zi/ = Z/i"^ . Expressing i>' and t in terms of : 

t = jl]v' = Zjj,~^ 



4 --^Z-7<Bi+2^7<Bi-22-'\/ X~-*'+T'2*i-2^\/x~-^~T'2*i+2 



Here the exponent in /i does not matter since we already have infinite order 
vanishing at whose boundary defining function is jl. Thus at the corner 
of £ and T, 



t^6,,- -pI' p^ 

Denote by £^66 the resulting error, since i>' = hz~^Z~^ and account for a factor of 
{u')"^ we have 

= QbUe + E<Si 

then 

F,b e p:^pf,P^P^A ' ' pV ' ;7si = 

Note that the above process introduces error with non trivial behavior back at 
so we will have eliminate it in the following subsection 6b, using the same argument 
as in step 4. However this is time is easier since we will be only dealing with one 
conjugated operator since there is one phase function ipi^ that we are dealing with. 

9.5.3. Dealing with the error at 5Bi. 

E^^ = Ef^ e p'^pI/t^ ' pI " C-(9Jl6,,.); 
The projective coordinates away from A are 

Z , 

In these coordinates Qb has the form : 



{z'v'f 



(H-r + (A-^)+P^V-P 
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The operator in the bracket on Si, Ptg^ restricts to a Bessel equation of imaginary 
argument of order — A)^/^: 

+ (A-^)-P 

On the product structure Bi = M.f x Kj^", , along each flow we can use the Macdonal 
functions to obtain the specific solution to the inhomogeneous equation with the 
following properties: 



. For 2 V - A ^ Z : at i ^ 0, fV X-^^o 



at f — > oo, ^ exp(— f) i S°{t)j:,'; ie rapidly deceasing, where <S'°(t)p' 



is symbol in t with parameter in . 
Expand in taylor series around *Bi. 

taylor expand 735 j jp ~ n 



We solve for of the form 



Note that line (60) would be replaced by: 

= (1 - X{ii>')) cxp-^*i+*-I^ Ufj + x{iv')uf.} 

Then take asymptotic summation to construct U^,: 

Ue = cxp(-(^^J[7f 1 



frSi „oo„oo oo oo ■^•3 1-2 V 4 ^ V 4 ^. ~_ _ 

^ PaP^^P^^Pt P^^ Pc Pn ' 7<Bi - " 
9.5.4. Summarizing. 

C/^ = exp(-^^ jVf ^ + Vf ^ 



QfcC/f ^ e PAPtPf.P'^Pt.p'P'^p'Ifi'^^'-^l^^ = 



exp(-z(pn,,)l/f 1 + V^^ 



•^6 ^ Pj^ P<& 1 Pc 



-2 

I 

63 



-2+7 



'Si 



■i/Si (z rP° rP° n°° f?^'^ ^ 
^% ^ PAPd2Pj^P<3i Pc 



Write = QbU^' + E<s, then 



P-R 



\ 1 
4 -'^ 7!Bi-2+2 

Pdi 



TP r- rtOO -OO -CX) -OO - V 

^66 e PAP^.zPj'P'BiPn 



PREP COMPUTATIONS used in calculating error The chosen phase func- 
tion in current coordinates and on tv' > 1 



2 \{ii>'z' + 1)3/2 _ ^ 1)3/21 

lim —-- -— 

z'^Q 3 Z'u' 

Prep calculations for f — )• oo so for u' ^0 , ii>' > 



idt 
{idtf 



2 1 

3^ 
2 1 
3z^ 



{iv'z' + 1)3/2 _ ^ i)3/2l _ _ i_ 



(fi>'-l) 



(tV/ + 1)1/2 ^ 1 
(tV^' + 1)1/2 _ 1 



f^v'z' 



2 + 1)1/2 



1 



Also on this region on which {iv' > 1} we have 

2 {iv'z' + 1)3/2 _ ^ 1)3/2 ^~^/ _ 1 



lim 

^'^■0 3 



exp - i(piii„ + {it - ^)) 
We have similar computation for ip^^ . 



1 + {}V'Z< + 1)1/2 2 {^v'z' + 1)1/2 

lim exp f-iyJiim + fii- ir)) -1 = 

1 = 0(^0 



9.6. Step 6 : removing error at '^\. The construction from step 6 introduces 
error on which will be removed in this step. We will need to work with the 
conjugated operator, conjugated by (which is (^in on this region). As a reminder 
we write out the form of the phase function near 5i 

<^U^ = |(^Z'^)-l [(^ + 1)3/2 - (^Z' + 1)3/2] 

The coordinates in this neighborhood and product structrure :Product 
structure 

[0,<5)p,^ x(5i),l. 



- 1/ , 

Coordinates : /i = — , Z , 



Pdi = 

By direct computation, the operator conjugated by exp(— (^nm) has the form : 
P- / v/n , c x2 , A* - n'^ 



2 P 



{~P) "'Hp) 



^py, i-z'dz' + zd^f + -^py, (A - -r) + 2i(^ + i)'/'(-z'az' + zd,) 



(A) 
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We are removing E^j^ = exp(— which comes from Eqi, = E^^ + E^ with 
the index indicates the support of the function. 



We will construct U^b (supported close to S^i) with 

Ueb = exp(-i(piim)i/66 



Expand F^b at S^i 



Ueb€pT,' ^^K^,p:^pP ; 7<b, =0 



„ taylor expansion 7*b 1 — 2+ ^ +1 

^66 ~ Z^PS,^2,k 



First transport equation 

2»(0 + l)i/2;j2 ^ ^^Q^ _ (~^^ _2+i+o) + ^ + i ^] Ueb;o = -Feb;o 
IC:z4U5nU=o = 



Higher transport equation 



2i{z + 1)1/2,2 _| [^g^ _ (^^^ _ 2 + i + i) + - + 2 ^]C^66;, = "^5,0 
[(^a, - (^^^ - 2 + i + i - 1))' + A - ^] t^66;i-l 

IC : ^^Us.o = 
Take Borel summation to construct U^b 

borel sum ct^, —1 \ ^ j _ otix, —1 a®, —2 oo 

Define 

J/gj = exp(-i(^iim)t^66 

note this is supported close to S^i. 

The resulting error denote by supported near S^i has the properties 



-C'S e Pgr^ pj: p^^ p-^ P£ PSi ; ISl = u 
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9.7. Step 7 : remove the error at TZ. We proceed to remove the error at TZ, 
which we denote by E-ji . 



TP Tp'R' I O TT^l r- ^OO ^OO rtOO -OO „00 ^OO -CX) ^OO ^OO ^ V 4 ' 

Ev. = Eq^ + QbC/g e Psi Psi ^s/ffi ^7?, P/: 

('Owr convention : TZ ~ {Z = 0}). Since we have infinite order vanishing at other 
the other faces except TZ and £, we can all blow them down and consider the error 
to live on [0, 1)2 x [0, l)z' x M^. The coordinates in use are then : 

h, z, z' 

In these coordinates property of E-jz : 

En e h°°z^^~^^\z')^^~^C°°{z,z') 



The operator in these coordinates 

4 



Q(/i, z) = h\zd,f + h\\ - ^) + ^3 + 0^ 



f^-\ 1 

One use power series argument to find U^ of the form U-j z^ ^ z'^Uj^k 
such that 



QUr -Ene h°°z°°{z')'/^^^C°°{z,z') 
specifically we solve in series for the Ur-j . 
Hence we can construct Uj with 



with resulting error iJg 



£^8 e h°°z°°{z')y 4 



9.8. Concluding the construction. We have constructed U such that 
U = ijo + Ui + U2 + U3 + Ui + U5 + U^' + Uf' + f/66 + Ur 



QU - SKid e h°°z°°{z')'J^^~^ 

To compute the overall asymptotic property of J7, we first relist the property of 
each of the component : 

"ffi = "osi + 1 + 1; Q!<Bi = 7!Bi + 2 - 2; 7<Bi = -2; 7^ = -1; 75^ = 7^ - i 

=752 + 1; = 5 + Q!52 + 1; ^3^1 = /^QSi + ^ + 1; ,S!Bi = q!<bi + 2 - 2 

7»i = 
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in a simplified version 

"51 = /35i =-2+^ + 1 = -i 
P'Bi = a<Bi = 7qSi = -2 
(61) 7^ = -l; -5^ = 1 

7^2 ~ "1 ~ 2 ~ ~2' "52 ~ ~2 

7asi = 

C/i = exp(-z^i„);7i, C/i e p^-^/p^ip^-^C^^CEKb,,); C/^ € p2V^V^C°°(Mf,,;,) 
C/2 = exp(~^^i„)c72; iJa G PrPs^V^'V^^C^^C^Jlft,,) 
(^in = /i-'l^(^) - ^{z% ifiz) ^l[{z + 1)3/2 _ ggj^^^ 

f/s = exp(-i(/3out)C^3 + V's 
^out = h-^ [ip{z) + ipiz')] sgn^„; ip{z) = f [(^ + if^^ - l] sgn^„ 

C/3 e p^^^p^ *^p^'^C°°($Hb,/i); supp J/a is close to ^ and away from 7^ 



c/4 = exp(i^o„t)C/4; e pl'^pr2p£v§^ 



U?'^p7P^yn^ 'pI^ V(97l,,;,); 7^ = 



ur&p'/p^yi''~\t'~'c^mb,H) 



[^interior g ^x.^^ 4 "^v 4" (OTb,;,) ; supp C/^"*<^"°'' Only intersects J^,TZ,£ 



In coordinates close to and £ : (/?^^ = | [ — (1 z)^^^ + (1 ~ 2;')"^^^] 



^6 ^ PaP^^P^^Pj'P^^ Pc Pn ' 7qsi - ^ 
C/e = C/e"^^ + Uf' 

Close to and C : i^um = Vin 

1 

•^6 ^ PAPs2pJ^P<Bi PnPc P^i 
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n . 71 . 



^6 ^ PAP^2PrP<Si Pc PtZ PSi 

J766 = exp(-iv3ii„i)C/6f,; t/ef, G Vsi' ^pJ'P?'; supp t/eb is close to 

[/7 e /i^^z/^^+'Cz')^"^"^ 



((^in and (fiout have support close to *Bi and A, ie we only have oscillatory behavior 
at the Zq > The term containing (pn^ is supposed close to ^i. 

9.9. Summary of parametrix construction for Q. Reminder Q = h^^Q 



(62) [/ = l^o + exp(-i.^i„)C/in+exp(-z.^out)C/out + exp(-i(^iii„)C/tp,i„, +C/bdy + t^<8i 
where 

C^in e pl''^'''%^\+'pY+^p"^l'^%''i'^'c°-m,^h); suppC/in is close to A 
7251 = -2; 7-4 = = -5; 752 = -§ 



(^vun. e PaPtP^^'pI ^ %sr'^'C°°(5Hb,^); suppt/^,..„ is close to Ji; = 



f/bdy e * %y 4 %-p0^^p|^p0^c°°(aHb,;,)C°°(Wl,,^); suppf/25,is away from A,^^ 



C/<s^ G exp(-^^ Jp^p^ ^ \i 4 'p|=^p|^^p|^p^C-(9Jtfc,,0 

10. The global resolvent fiberwise on Fc 

As motivated in section [Tj we want an exact solution to Qu = on as well as 
information on its asymptotic behavior. In order to obtain this, we first construct 
an approximate solution with an error that is trivial at every boundary face of Fc 



(see figure 20 1 and in particular is semiclassically trivial i.e 0{h°°). We will need a 
global resolvent on Fc to solve completely this resulting error. In this section we 
show how such a global resolvent is constructed. 

We will construct the global resolvent fiberwise. On the fibers where \Zo\ is 
bounded, the problem is elliptic and is not semiclassical, so we can patch together 
by partition of unity the resolvent at z = and that at infinity ie at = 0. 
However, since the operator is not semiclassically elliptic, this method will not give 
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uniform bound as Zo — > ±00 and will result in error that is not semiclasslically 
trivial. When Zq — !■ 00, we use the gluing technique by Datchev-Vasy in [1 to 
combine the readily available resolvent estimate for the infinite region given by 
Vasy - Zworski in 18 with an estimate we will obtain from the local parametrix 
constructed for the region close to in section [9j The situation for Zq —00 is 
much simpler since the gluing is taken over an elliptic region, i.e. there will be no 
semiclassical singularity propagation. See remark ( |7.2[ ) for further discussion on 
the difference in methods between the two regions. 

The proof for the main lemma will contain 3 ingredients: estimate for the infinite 
region, estimate for the local region near and the gluing process. 

Lemma 10.1. [Mapping property of the global resolvent on ¥c] 

Partition ¥q by Xq and Xi where Xq — {z > 1/4}, Xi ~ {z < 1}, Xi (trapping 
region) XIjXo partition of unity, ipj ^ 1 near suppxj- 

QGi = Id+Ei, pill e 0(/i°°), IIGill < ai = ci/i"™' 
QGo -Id,||Go|| <ao^coh-' 



where asmi is given by (66). Here the norms are taken in some weighte^^ L'^ 
spaces. 

Denote by Osc^^ the set of oscillatory functions on Fq (see figure with the 
behavior: 



(63) 4 ^/i-C°°([0,l),)+exp((/.„,,)/i°°C°°([l,oo),) 

i.e. having boundary asymptotic prescribed by X at z ^ and oscillatory behavior 
at z ~> 00, where in the corresponding coordinates used on each region we have 

r-i|/i-i[(z + 1)3/2 -l]sgn(?„ Zo>0 
(l^osc^ l-lh-^[t{z-lf/^ + l\ Zo<0 

[-«|Z3/2 sgn0„ \Zo\bounded 

We obtain the mapping property of the global resolvent Gq on ¥c 

G^:C°°{¥c)^ Oscw^- ||G|| < Gai 

close to z = the weighted norm is that same as that for Gi, while close to 
infinity we use that for Gq . 

10.1. Zq > 0. 



10.1.1. Ingredient 1: Semiclassical estimate at infinite. 



h^i^d^f + e'^v+ - ll ; y+ hH'^/^d, - §(A - ^)hh^'^ - 1 



Notation: In this subsection for matter of convenience we will work with variable 
(7 = z3/2 and t = (7~^ instead of z. 

Vasy and Zworski in (|18)) prove a semiclassical estimate for asymptotically eu- 
clidean scattering with long-range perturbation. 



see estimate (|64| for weighted norm for Go and see 
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for Gi 



Theorem 10.1 (Vasy - Zworski (2000)). Let X be a manifold with boundary. Let 
A be the Laplacian of a scattering metric on X. If P = h'^A + V is a semi-classical 
long-range perturbation of h^A and R{X) = (P — A)^^ its resolvent then for all 
m G M 

||i?(A + zi)/|| <Coh-^f\\ ;e>0 

With Co independent of t ^ real A € /, / C (0, +oo) a compact interval in the set 
of non-trapping energies for P. 

Here H'^'^{X) denotes Sobolev spaces adapted to the scattering calculus. The index 
m indicates smoothness and k the rate of decay at infinity. 

Remark 10.2 (Long range perturbation). Long range potential V are second order 
semiclassical operator which in local coordinates of the form V — J2\a\<2 ^a(^j h){hDz)°' 
and near the boundary dZ, in local coordinates y € dX^ 

V = x-^Vq,Vo^ Vko.{x,y,h){hx^D^f{hDyT, 

\a\ + k<2 

Vka - < e /i^"'"^°(^),< e 5"^°(X); 7 > 0. 

In Euclidean setting the condition on the coefficients corresponds to assuming that 
the coefficients are symbols in the Euclidean base variables. The constraint on 7 
is just 7 > 0. For application in our case, 7 = |, ie the pertubation goes to zero 
more slowly than Coulomb potential. 



Denote by 



In order to apply the result we need find a global operator on Mo- whose behavior 
coincides with that of Q+ on cr > 1. Q+ is extended to IR+ in a way so that: there 
is no extra boundary at cr = 0, the problem should be elliptic on cr < ie the 
negative infinity end, and most importantly the operator is globally nontrapping. 
The following extension satisfies these criteria: 

QgiobaL+ = x(ct)Q+; Q+ = 



XGC-(M), 0(a) = 
is such that V^~''|ct>i = and on h — 0, has the profile of 5+ in figure 




16 



In particular -5+ = c > | for cr < j, > | for cr < ^ and < | for cr > i 
and on cr > 1 we have = . Next we need to verify that the semiclassical 
bicharacteristic of Qgiobai.+ is globally nontrapping. 
Non-trapping condition verification: 
The semiclassical symbol of Qgiobai,+ 
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Figure 16. Extending the local problem of Q-|_ at infinity to a 
global problem on cr = z^l"^ 




Figure 17. The bicharacteristic set of the new global operator Qgiobai;+ 



Since > 0, on the characteristic set f& we have ^ = ±i / % — v. x G C°° so the 



bicharacteristic of ilf,+ is the same as that for . 



± 




When ^ > 0, a{f) increases. When ^ < 0, a{£) decreases. Also the bicharacteristic 
set is a one dimensional submanifold, so away from ^ = the bicharacteristic is 
non-trapping. It remains to verify non-trapping condition for cr < |. 
In fact, close to time s = 0, we can work out directly the profile of the bicharac- 
teristic set. a can be considered as a smooth function parametrized by ^. This is 
seen as follows: v — ^ — ^ and on < tr < | so {i is a smooth strictly monotone 
function of cr so cr = f (w) for some smooth strictly monotone function / and hence 
cr = f(| — is a smooth function of ^. The profile of the characteristic set is as 



shown in figure (17) 



Now we are ready to apply the result by Vasy-Zworski to Qgiobai.+ i denote by 
R{\ + it) = QgiobaL+ we have: 

\\R{X + it)f\\ 1 <Coh-'\\f\\ 1 ;e>0;A=| 
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Following the method in Melrose's paper using these estimates, one can show 

1 „ 1 

that for fixed /i > the limits R{\±iO)f exist in Hsc ' ^ for / e Hsc ' ^ (^o), 

||i?(A + iO)/|| 1 ^ <Coh-'\\f\\ ^ ^i^^ ; e>0; he{0,ho) 

In terms of z and m — 2 
(64) 

\\R{X + iO)f\\ 1 , <Coh-'\\f\\ 1 ; e>0; A= |; fte (0,/io) 

Q = h^^Q. Denote by 

Go-^i?(|+zO)/i\ 

where x = 1 on cr > 1, ?/;= 1 on suppx, and suppx C supp?/; C {cr > |}. Then, 
with the norms taken in corresponding weighted space, for some ao(/i) > we 
have: 

(65) IIGoll < coh^ = ao 



Mapping properties of Gq 

Go :C°°(Fc) ^/i°°0scoo 

where h°°OsCao is the set of functions that are semiclassically trivial with trivial 
behavior at finite end (at ^ 0) and oscillatory behavior at the infinite end (as 
(T — oo). The oscillatory part is given by: 

exp(-z</.)C°°(Fg;+); 0(z) = f/j-i [{z + if/^ ~ l] sgn0„. 



10.1.2. Ingredient 2: Estimate around from local parametrix o/Q. 

From the construction of the local parametrix near 2: ^ in section ^ we get: 

QU-SKm = E£ h°°z°°{z')^f^^^^'^''^ 

U ^Uo+ exp{~iipin)Uin + exp(-Z(pout)C^out + exp(-«^iim)J/ipiim + '^bdy + f^<Bi 

where C/jgi ^^e portion of the approximate solution near *Bi (where Zq < 0). Uq 
and SKid have support close to the lifted diagonal, away from C and TZ. (See (62) 
for more details). 

For the purpose of weighted estimate computation, we write U as 

U = f/ifdo + C^osc 

U<s,do comes from the term Uq which is supported close the diagonal and is conor- 
mal of order —2 to the lifted diagonal, while Uosc come from: exp(— i</3in)J7in + 
exp(— i93out)t^out + exp(— i(y9iim)L/;pi;^. Since we will only need a weak estimate and 
deal with L°° norm, we can ignore the oscillatory terms. 

(1) First we deal with the term C/*do- Due to the cornormality of negative 
order at the diagonal we have 
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We rewrite J7*do as 

f/*do = hx{z)x{z')f{^^), f e L^; suppx C [0,1] 

sup / |C/*do| dz' ; sup / |C/*do| dz < oo 
this is because 

hx{z)x{z')\fC-^) I d,J f x{z)x{z - yh)\f{y)\ dy < 



Similarly 



hx(z)x{z')\fC—^)\dz' J' f x(z)x{z' + yh)\f{y)\dy<^' 

(2) We are going to compute which power of h and z' that will make Uin, Uout 
and L°° on ^b,h- In fact, to make sure the operator whose S.K of U 



maps into zV 4 '^C°°([0, l)z), we find c and a such that: 



(z')° lifts to 
h" lifts to 

(p® 1 P<B2 PS2 PSi Pa) {^b,h) 

z lifts to 

z = pnPTPd2P'SiC°°{^b,h) 



From (62), Uin has the following property: 

TT ^ 7'Bi+2 —1+2 74+2 asi+2 752 +2/ioo /rm \ rr • l i /I 

i/in e P<B2 Pji Psi /^52 suppC/in IS close to ^ 

7251 = -2; 7^ = -l;a5i = -^;752 = -§ 
To take care of the index at 

• *8i : we require a + c> — A 



Similarly for C/out 

and for U,^^:^^^ (note that the hidex of pc is aheady positive) 



C^vu„ e PaPtP^^'pI ^ " P5^'^2c°o(9j|^^). suppt/^,.^ is close to Ji; = 
Choose 



Denote by C°°(Fc) the set of functions which vanish to infinite order at 
both the finite end (where z = 0) and the infinite one ( where z~^ =0) and 
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is 0{h°°) at both end when Zq — >■ ±c». Such a function has the form for 
any N and N' 

e^h^'{z'feN',N{h,z'); e(z')C°°(Mj x [0, 1),,) 

{U,e{z')) = {Uh^{zr,e,A^y)) 

with a — \J ^ — A, c = \J — A we have 



The blown down of h'^U e zV 4 ^L°°{[Q, l)u x [0, l)^ x [0, 1)^0 thus L\^^ 
and hence satisfies the requirement of Schur's lemma. For e G C°°{¥c) we 
have: 

||(cxp(iv3i„);7in, e)|| r-2 — <Ch^''\\e\\^aL2 



4 i2([0,l],) 

The same computation applies for C/out and C/,^,i,„. Thus the operator is 
bounded between the following weighted spaces : 

(C-(Fc),|l.|U.i.)^z/?-^L2([0,l],) 



•II < /i"--; a = V ^ - A, = c = V - A 



(66) 



Similarly we have the bound for the error term E € h°°z°° (z') V 4 with the same 



weighted space but with 

\\E\\^o{hn 



10.1.3. Global resolvent estimate. Now we are ready to put together the two above 
estimates using the method of Datchev-Vasy in pi . 

Denote by = {cr > 3/4}, Xi — {z < 1}, Xi (trapping region). Let Xi,Xo be a 
partition of unity, with ijjj ~ 1 near supp Xj ■ 

Ingredient 1 (651 gives the estimate in the weighted bound z~2'^2~ 
z2^2+^'>L'^ and 

QGo-Id, WGoz^^^W < coh-' ^ ao 
Ingredient 2 (661 gives the weighted estimate from z^L^ x°'L^ with 



QGi^ld+E; \\Gi\\ <cih" 



ai; \\E\\^Oih°^) 



Define 



F := iPqGqXo + ipiGiXi 
QF = Id +^iExi + [Q, ^o]GoXo + [Q, ^i]GiXi 



Id +^3 + Aq + Ai 



We have 



Al = Al=Q 
Q(F - FAn) = Id +^3 - A3A0 + Ai- A^Ao 
Q{F- FAo - FAi) = Id +A3 - AsiAo + Ai) - AiAo - AqAi 
Q{F-FAo-FAi+FAiAo) ^ld+A3-A3{Ao+Ai-AiAo)-AoAi+AoAiAo = ld+E 
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potential v 








(7 = 1 - 




1 1 1 1 

(T = ^ 



Figure 18. Extending the local problem of Q_ at infinity to a 
global problem on a 



A3 - ^3(^0 + Ai- AiAq) is 0(/i°°) and by the argument in the paper ([!]) A^Ai 
is 0(/i°°), so is AqAiAq. Hence E g 0{h°"). For small enough h the inverse for the 
term exists. The global resolvent is then 

R^{F-FAq- FAi + FAiAo) (Id +E)-i 

Consider the bound for 

F - FAq - FAi + FAiAo = F - ViGiXi^o + VoGoXo(-^i + ^1^0) 

to get 

||i?|| < C(ao + ai + 2haoai + h^aQUi) < Ch^a^ai = Cai 
with the norm near infinity is taken in the same weighted spaces as for those for 
Go, while that near zero we use the weighted space for Gi. 

Mapping properties : The mapping property of R is given by that of F = 
V'oGoXo + '/'iGiXi- The oscillatory behavior at infinity of the image is given by Go 
while the boundary behavior prescribed by A at z ~ is given by Gi. Hence R 



maps to the space Oscf^ as defined in ( 63 1 
10.2. Zo < 0. 



The semiclassical scattering principal symbol is 



4+2/3 



10.2.1. Ingredient 1. Denote by 



- _ 4.2/3 _ 4 -2/3 

7 gL gU . 



We find a global operator on Mo- whose behavior on cr > 1 coincides with those of 
Q_. Define 



Qgiobai;- = x(o-)Q-; Q- = 

XGC°°(M), 0(a) 




a > 1 



a > 1 



V" is such that V^\cr>i = and on h = 0, y~ has the profile of in figure 
(18 1. In particular the problem is semiclassical elliptic for a < 1. Qgiobai;- satisfies 

75 



Figure 19. The bicharacteristic set of the global operator Qgiobai;- 

global non-trapping with the profile of the characteristic set given by figure [19] 
After this step we follow as in ingredient 1 for Zq > 0. 

10.2.2. Ingredient 2. Similarly as in Zq > but in this case it is even simpler since 
we have exponential decay at the semiclassical face. 

10.2.3. Gluing the two above estimates to get global resolvent estimate: Since the 
gluing region is in [i, |] where the problem is elliptic, we do not have propagation 
of singularities. 

Denote by Xq — {a > 1/2}, Xi — {z < 3/4}, Xi (trapping region). Let Xi^Xo be 
a partition of unity, tpj ^ 1 near supp Xj ■ Ingredient one gives the estimate in the 
weighted bound 

z-i(5+^)L2^Xo) ^ zi^^+'^L^{Xo) 

Define 

F := VoGoXo + "01^1X1 

Q^^ = Id +2P1EX1 + [Q, 4'o]GqXo + [Q, ^i]GiXi - Id +A3 + A0 + A1 

Since the problem is semiclassical elliptic, with [Q, ipo] and xo having disjoint sup- 
ports, Aq = 0{h°°). Similarly for Ai. From here we procceed as in Zq > 0. 

10.3. Zq bounded. Similar to Zq < 0. 



11. Construction of approximate solution on the blown-up space C: 

DETAILS 

In the following subsections, we fill out the details of section [7] and construct the 
approximate solution to ^ on the blown-up space C. We refer the readers to 
section [7] for motivation and the main idea of the construction. See figure [7] for 
blown-up space C. 

11.1. Step 1 - Approximate solution on Fq- 
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11.1.1. Step la : Getting the local approximate solution from the local parametrix. 
For convenience, we list the properties of interest for the local semiclassical parametrix 
constructed in ( 62 ) 



QU - SKid = Ee h°°z°°{z')^ 4 ^ 

U =Uq+ CXp{-i(pin)Uin + exp(-i(^out)C^out + exp(-i(piim)[/i^ii„ + C/bdy + 

where Uq, Ui^ and C/out are supported away from C, hence only these following terms 
are of interest for the purpose of ( 67 1 



Uvu„. G PaPtP^VpI " Psr 'C°°(9Jl,,h.); supp[/^„.„ is close to 5i; = 

where close to and C : ipuni — ^in with ipia = h^^l'fiiz) ~ ipiz')\, ip{z) — 
I [(z + 1)3/2-1] sgn0„. 



— —A 1/ — — A 

f/bdy e " Pc pf^p'^^p^^p°^C'^{Mb,h); suppUrQ^s away bom A,di,A,Si;i 



— — A 

U^^eeM-V^,)PAPl" pI' p'f^pf/^/:FC°^{m,.n) 

^^^ = h-^mz)-^{z')\- ^ = |[i-(i-z)3/2] 

We want to get an approximate solution up^ on [0, 1)^ x M^o such that 
(67) 

'Clu^.^^e^o^(.h°°z°°C°°{¥l) 



= pint ^t" + pJi) V ' ^ith u^o |p = 1; G C°°(F^)(*) 

See (20) for the definition of ^^(i) and ( |T7| for the definition of z. 

Remark 11.1. If we modify IJ by factor dependent only on h and z' then the 
resulting term still satisfies Qw = up to a trivial smooth term. 

In order to to this, partition the front face Fc into ¥cfi^ IFc,+ — Fg'^UFj^, ^, Fc _ = 
F0,__ U F^__ as pictured in figure ([20]). 

• Consider the portion of ?7bdy = close to £, denoted by C^bdy Close to C 



PC - Z; Z-VX-^ ^ z-Vx-^(/)^/^-^ 

/i lifts to p^iPSiPAPd2P<Si- ^' lifts to pjrp<SiPcPdiP^i- z lifts to PTP<BiPnPd2P<Si- 
So we have: 



V— \ — A/^i~^'^ ^^i^A/^;~^A A/— p— A A/^i~~A „ 



/l"(z') V 4 -[/j^^y ^p^V 4 ^^^^ 

here we list only the exponents of faces that do not have infinite order 
vanishing. 



If we choose m = y ^ — '^j then Ubdy is smooth down to C and TZ with 
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4 



Rtgion F^^ 




Region 1F^_ 






Region ¥c,o 




5 


Region ¥°n + 


Region 




u ■ 




.1. 





Figure 20. F^°j_ denotes the neighborhood in region ^ close 
to the corner of t = 0, /i = where t is the projective blow-up 

0,i-i=0. 

Similar notation for F^°^ and F^°_ for the other side where Zq < 



witht = f , while Fp°;; 



of Fr 



For j = 1, 



is the neighborhood close to z 
'°c^,°!andF^°_ 

, 6 the numbers indicate the label of the 



corresponding boundary face F^-* . Denote by F^ the region of ¥c 
away from boundary face F^ ■* . 

See (20) for the definition of p^ii) and ( |T7| for the definition of z. 



"c 



See (19) for definition of a. 



value in distribution on DJlb^h- Note Ubdy|.Fn£ is smooth in h. Denote by 
Cbdy ^ smooth function on DJlb h such that Cbdy is constant for h constant 
and its restriction to J-" n £ equal Wbdy|j^n£- Constructed this way, Cbdy 
only depends on h. The following term 



lim>'"(z') V 4 



satisfies the boundary condition (*) in (671 ie at F^ , smooth down to Fj- 
and F^?'' , and is supported away from fJ^'' . 

With this choice of to, consider the asymptotic properties of 



Jim/.™(z')~V 4-\^i^C/^^^, 



We have: 



With the choice ofm = \/ A we have 



lhn^/i"(z') V 4- 



z'-)-0 
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is supported in region ^, away from f[^' and F^^-* and at boundary faces 

Fp-* and F^' it has exponent 7>8j + | , 7231 — respectively. 
Similarly, 



lim^/.'»(z')"V 4-\^i^{/^^ 



is supported in region F^ _, away from F^-* and ¥q\ vanishes rapidly 



at 



(5) 

and has exponent 7(3^ at F, 



Summary 



z'->-0 

Ufo satisfies (67) and 



2 / 2 



with 

• Mfo,jj lives on [0, l)z x 

• Uyo^ ^ is supported in F^ ^ , in more details 

_ 1 / . / . 

wfo^_^ e exp(-z0)p^(,,Vp(3) (F(^^^) + Pj^(,,pJ(,) C (F(.^^)+p^(,)p^(/) C (Fp_^) 

where the phase function comes from Ihiiz'^o ^Pnm and in coordinates cr, h 
and I', cr has the form: 

MfO^ is supported in F^ _ , in more details: 



%o^__ e exp(-^)p^(e5jC (Fc^_) + Pj,(6)P;(i) C (F<.__)+Pj^,e,Pj^,;, C (Fp^_) 
We have a local approximate solution u^o on this region with error e^o : 



— — A —A/— — A 

"FO = Pjdf "fo + V ; with U^o Ip ri)=0 = 1; Mpo G C°°(F*J,)(*) 
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11.1.2. Step lb : solve away e^o^ along the global transport equation in a away to 
infinity. 

Remark 11.2. The computation in this section is carried out in a, h where z = a^l^ 
but could have been done in z, h. 



For the region Zq > : We deal with the portion of eyo^ that is supported close 
the corner ofF^^ and . To get the eikonal equation and the transport equations, 



we conjugate our operator by the oscillatory exp(— z^): 



= (j'^dlu + {2h^^i{adcr(f ) + b) ad^u 
+ 



e "-hu; 6=1 



We chose 



2/3 



1 



,3/2 



1 



sgn 0„ 



These computation will be useful in determining the behavior of the regular singular 
ODE at = 

d,^ = {a^^ + 1)1/2 |<,-V3 ^„ (-2/3 ^ 1)1/22 --l/3sgj^^^ 

= |(<7V3 + l)-V2<,-2/3 3gn^^ ^ 2(^2/3 + 1)1/2^-4/3 ^g^ 

20cr(p 6 

We will look for a correction term upoo with 

c,+ 

wf~^ = exp(-i(/))/i"3 n3 = 7<Bi + 5, ^<Bi = and 

Pwf-^+Cfo ed°°(Fc,+) 

Expand e^o at F^' 

e^o, ~ exp(-i(j!))/i"^ ^/i-'cfo 

Apply the conjugated operator to ft,"^ ^ hPup'^^.j and gather the terms according 
to power of h (starting from lowest) and set them to zero : 

• Eikonal equation is satisfied by our choice of 1^ - coefficient of /i"^ , ie the 
term in the bracket vanishes 



-2/3 



= 



since 



^ ,)4^-2/3 ^ 4 ^ 4^_,/3 _4 _ 4^_,/3 ^ 4 ^ 4^_,/3 ^ , 



• First transport equation - coefficient of /i"3+i 

4 4 

^'MF~_^;0|i7=0 = 
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where epo^.g vanishes rapidly at both ct = and a — > oo. Since ip satisfies 
the eikonal equation, we obtain the first transport equation for u^oo .q: 



D Z 



the solution has the asymptotic property upoo .g cr 2 with 6 = 1 as 
a — )• +00. 

• Second transport equation - coeflacient of 

^WFg>_^;l 1(7=0 = 

where epo vanishes rapidly at both tr = and cr —> 00. Substituting in 



the expression for ij9<j(^ we obtain 



ad, - ha^l^ + 1)-^ + \h 1 WF^^;i - I 



(aa<,)2 + (6-l)aa. + |(A-^^) 



2adaf 



;1 



[(<t9<,)^ + (6-1)ct9, + §(A-^)1 _i 

Since ad^'^ ~ cr and -^^ 2ad~^ ~ a 2 the asymptotic 

behavior for itp?? -i at cr"^ = is 



Similarly, we obtain higher transport equations : 



2id„ 



(a5<,)2 + (6-l)a5<, + |(A-^^) 



2a 



= 



L'WF~_^;j 1(7=0 = eFO,;j 



where ej-u.y vanishes rapidly at both a = and a — )• 00. Similar computa- 
tion to get 

_ 1. . 

u^x^-j a 2 •';6=1 ascr^oo 

Uvg^ is now obtained by borel summation and has the following asymptotic prop- 
erty: 

Ur<x^ = exp(— i(/>)/i^®i~^ h-'u^'^ ^-j; u^^.+u ~ a~^^~-' as cr — >■ 00. 



For the region < : We deal with the portion of Cjo that is supported 
close the corner of f[^' and With <fi=l [iia"^/^ - + l] we compute the 
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conjugated operator 



ehQ_e heh 



+ 



•2c,2 , ,.2'9<T 4 4 4, 



-2/3 



(e 6 = 1 



9 9 9~ 



We will look for a correction term wpoo with 



exp {—<j))h'' /f-' Uf^ _ -J ; ^ = ^ z(cr 



2/3 , 



1) 



3/2 



'C- — -^^ ^ - -"c,--" 3' 

and Pup^_ + epo^ € C°°(Fc,-) 
note here since the previous error term is h°°, s is arbitrarily large. 

We obtain the eikonal equation and the transport equations: 



Eikonal equation : Our choice of = | [(1 — cr^/^)^/^ — l] satisfies the 
equation 



(69) 



-4/3 



+ 



^0 = 



since 



-d ? 



^(^2/3 _ 1)3/2 ^ ^ - 1)1/2^-1/3 



- 1)1/2,-1/3)' _ 4^-2/3 + 4 ^ _4^^,/3 _ 



2/3 



-2/3 



= 



First order transport equation 



WFg'_;0 



= 



,'"Fc,_;o|ct=o 



the solution as the asymptotic property mf~_;0 ^ o" 2 as cr ^ +oo. 
Second order transport equation 



a9. + J(a2/3-l)-i + l6 
6 A 

MF^__;i|<t=0 = 



WF~_;1 + 



(a^ag-6aa^ + |A)MFg_;o 



= 



we have 



(<73^-b3„)ao ^ ((<T8^)^(l-b)<79„)ao 
29,, 2adc'p 



a 2 



6-1 



^b-1 



so 7iF^ a^cr 2" "ascr— >-oo 
• Higher order the transport equation for s + 2 + j 

WF~_;j-M + 



= 



WF~_;j 1(7=0 
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Thus Ufgj :j ^ cr 2 as cr — ?> cxi. 
Similarly for when Zq is bounded, this we work with 

(Sdsf + {b- l)Sds + A + 5=^ + ZoS^;b = 1 

and conjugate the operator by exp(— + Zq)^^'^ sgnOn) to obtain a regular sin- 
gular ODE and whose solution has the asymptotic property as — > +oo 



. 2 



11.1.3. Summary of step 1. We have constructed an approximate solution sup- 
ported in Fc.+ that is polyhomongeous conormal to all boundaries of this face, 
with an error vanishing to infinite order at all boundaries of Fc. The polyhomoge- 
neous behavior of mf^ + in more details: 



(70) 
(71) 



;7rc,.gexp(-z^)p;r.)V;rr^C-(F^-) + p^,.,pV,r C-(F°^°+)+p^<.,p^7/" C-(F^°^) 



+ exp(-z</.)p3;i,^ + V,'(.,C°°(FS+) 



1 1 

F 



c 



J7^-.- eexp(-f)p2jC-(F°-) + P°,e,p;ar C°"(lF°c.-) + P°(.)P;j, 



+ exp(-0)p;;?,,p2„C-(FS;_) 



f/^--° e ,^ C-(F°,,o) + P„,y, ' C-(F°,,o) + e-'3 



p2 ,c~(Fg;o) 



where 



= i [1 - (1 - z)3/2] ; 4> = Ih'^ [iz + 1)3/2 _ 1] . ^ ^ 2/^-1 ^^(^ _ y 



3/2 



1] 



Remark 11.3. The oscillatory behaviors agree on overlapping regions. The index 



at Fp-* corresponds to the choice p (4) = (t ^. Sec (19 1 for definition of a. 



11.2. Step 2 - Justification of the asymptotic properties of the exact solu- 
tion on Fc. With [/p^ being the approximate solutio n on F g that we constructed 
in step 1, the exact solution on F^, from Lemma 



10.1 



we have 



Fc 



GE^ 



where GEh lies in the space of oscillatory functions which have the same asymptotic 
behavior as Up^ and GE^ = 0{h°°). 
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11.3. Step 3 - Approximate solution to ([s]) on blown-up space C. 

(1) We extend the exact solution Ufc into the interior of the space for \9\ > 1, 
using the fibration structure over front face (see figure [?]). 

(2) We need to modify the extended function so that satisfy the boundary 
condition in ([s]) in terms of x. As defined in ([l7]), 



z = XoiZoKZ Zo) + X+{Zo){^^) + X-{Zo){^^) 
= Xo{Zo)iZ Zo) + X^{Z,){^^)Z'I^ + ^_[Zo){^^){-Zof' 

Zo = \en\"^{i - \e'f); z-Zo = x\e^\"^-e' = f 



We have p^(i){x = 0,6) =0 hence %^|a;=o = 1- See (68) for definition of 
Define 



(72) 



•/-correction := CqXO (^o) I | ^ V 4 ^ + c+X+{Zo) |e„|(l " \e'\Y^^ 



— X 
4 



+ c^x-iZo)[\en\i\9'\'-i)'/^ 

Co, c_ and c+ are smooth bounded function of Zq (used to make up for the 
'piling up' of the values o/x+,X-,Xo)- 

Remark 11.4. For singularity computation purpose, we will be concerned 
with the growth in terms of \9\ of the various terms in the definition of 
^correction- On the the Overlapping region between xo ^-nd x+ both terms 

are both growing by j^nj 3 '-2 V 4 \ This is because on overlapping 
support of Xo a-nd x+ ^0 > is bounded ie 



1 - \ef < \0n\-'^' ^ {\e„\{i - \ef)'/')y 4 < i^jsVx 

Similarly for the overlapping region of X- and xo- 
Define 

n ^ 

(73) = Xi\0\) ^correction X 2 Uf^ 

x{\0\) e with suppx C {\6\ > 1}. Hence x is to cut off the exact 
solution away from finite region (\6\ < !)■ So in fact Ux is an exact local 
solution near infinity |0| > 1 (however an approximate solution on the whole 
blown up space B We have Ux satisfies: 




s.iX) = ^^±^-^X 



€ eM^<l^in)Cph{C) 4>n. = -le-^ [{l + x- W?)""^ - (1 - l^f )=^/2]sgn0„ 



^'^Since we only care about the singularities behavior at infinity 8, an approximate solution in 
this sense (cut off away from finite l^l) suffices 
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Remark 11.5. To get an approximate solution to the original problem ([2]) we will 
take the inverse fourier transform of U. In the singularity computation section we 
will work with J-^^{U). 



12. Singularities computation 



In section 11 we constructed iix that lives on the blownup space C (see figure [t]) 
and satisfies 



x'^-Ux \x=a = 1 s±(A) = - ± Y — - A 

Ux e exp(*0i„)/:p„(C) 0i„ = -|0-1[(1 +:e - |0'|2)3/2 _ (1 _ |0~'|2)3/2],gj^^^ 

The function Ux lives on blownup space C but due the reduction process in section 
|6] and the construction in section |11[ it has a complete description as a polyho- 
mogeneous conormal function on the front face of blowup space C (see figure 
[rjand 20). In this section, we will study the regularity / singularity properties in 
terms of conormality and wavefront set of the approximate solution U = J^~^{ux) 
on each of the region which after Fourier transform corresponds to the partition of 



Fc (see figure 20). Using these results, in section 13 we will draw conclusion on 
the singularity structure of the approximate solution U and then the exact solution 
[/ to ([2]) in section 13 which is the original goal of the paper. 

The main observation from the singularity computation: The portion of of 
the solution in the interior of Fc under the fourier transform is 'trivial' i.e. does 
not yield singularity. The structure of the solution away from a; = is due to the 
asymptotic behavior at scattering end for different value of h (near boundary face 

f(^)). In our case, due to the fact that we have a reduction to Fq, singularities at 



the glancing point propagate into the interior (i.e away from x 
to singularities at the hyperbolic points. 



0) and contributes 







flegion J^'_ 




Region Fc,o 




Region 


Region 










h h 





For convenience we list the polyhomogeneity behavior of the constructed function 



tlx defined in ( 73 ) at each boundary face of Fc (see figure 20 ) 

n ^ 

= Xd^*!) ^correction X 2 



where 



'/'correction — CoXO (-^o) l^*™ P 



-A 



c+x+{Zo){K\{l~\W^T " 



Q'|2 



1)1/2U 4 
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We relist the asymptotic properties of Ufc at each boundary faces of ¥c (see also 



ir^ir^ ^C^C ^c^'^ 

+ exp(-z0)p;r3,^ + ^pJ„C-(F^+) 



— —A —A/— — A 



+ exp(-(/.)p;^„pi,C-(F-_) 



where 

^ i [1 - (1 - z)3/2] ; <^ = Ih-' [{z + tf' -1];4>^ Ih-' [^{z - if^' + l] 

" c 

We will be consider the properties of the approximate solution U on each of the 
following region: for > 82 (see (17) for where 61 appears): 

(1) |Zo| < 2^2, Ifi"! > 1 : denoted by region Fc,o- 

(2) ^0 > 1^2, l^*! > 1 : denoted by region ¥c,+ which will be subpartitioned into 

• F^ ^ = Fp°^ U F^°^ localizing the two corners of the blown-up front 
face. 

• F? 



Zo < -(52 

\Qn\ > {i-25,)\e'\ 

denoted by Fc,_^i which is subpartitioned as above into 

00 " 



(3) Zq < —62, \d\ > I '■ here we will need to consider 
yhich is subpartitionec 



\ z <^ s 

and < denoted by ¥c - 2 which is subpartitioned into 

\\en\ < ii-8iW\ y ,2 



andF^°_ 2. 



To localize to each region we will use smooth compact function ipfo^ ^, V'Fg'Q, "fpfO-o j 
V^po,oo , V'F-^, ^w°^°_y V'f^-^^ , V'F~_ ^w°^l_^ and tpv°^^y 

Notation : 



Denote bys_ = | — y^ — A, a = l — n, a=|a+|=— + 
We will use the projective blowup variable s = \ z near the corner of F^^^ and 
F^?-* and t = j^,h near the corner of F^^ and F^?-* . Similarly for the other end 
where Zq <0. 



12.1. On region ¥cfl- In this region we have {\Zo\ < 262', \d\ > 1} which lies in 
the conic set {\9\ > 1,1 — ^0(^0) ^ l^'l < 1 + '^1(^2)}- Hence |6'„| is equivalent to 
\6'\ and thus to 16*1. Thus also implies that |6'„|, |6''| > 

For the subregion Fji, q ie close to cr = 0, we will show 

(1) U"^'^^" e x'^-C^y + x'+C^y on yn ^ 0. 

(2) U^c,o is smooth on x > S for any S > 

For the subregion Fg^Q ie close to ct^^ = 00, we will show 

(1) U'^'^.'' is smooth for yn < and vanishes to infinite order at a; = 

(2) smoothness and and vanishes to infinite order at a; = for ?/„ > if we are 
close enough to the boundary. 

(3) Wavefront set statement away from a; = for ifc.o ■ follows bicharacteristic 
strips 



12.1.1. On region F°,o. 



Remark 12.1. As defined in (17), and since S2 > 62, z ^ smooth function(Zo)S' 
where this smooth function is = 1 on < ^2, S = Z — Zq = xlOn]"^^^ ■ On the 
current support Zq < 2S2 so the properties of the solution in S are the exactly the 
same in terms of z. 



For x{x,y) e C^{x,y) 

(74) (U^c,o),^(^,y)) = iim(27r)-" / / U^o u,x{x,y)w{ee)e'y' dxdydO 



' C,0 

where we recall the form of Ux on the region currently considered 



2 nf_ . ^ . 



where U^"^ '' is conormal in S* . In more details: 



2 fnE i 00 4 /r? ^ 



Using this form of Ux we rewrite ( |74[ ) as 

(U^c.o,x(a;,y)) = lim(27r)-"^x^±y J J arc,,^.^^xix,y)w{ee)e''-Uxdyd9 

Lemma 12.1 (Properties of ^.j-). 

1 d 

(1) For L such that *L — -j — p- ^ yj -qq > ™2 have for Ck is independent of x: 
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(2) Onxy^O, 



Proof. 



(1) First we consider the bound on afo^^.j.. Note on the current region is 
equivalent to \6\. 



2 /n£ , „ 



2 /rr X 4 /rr , 



^ laFO^^^il < |0„|3V 4 

To compute the bound for higher derivatives, we use : 

de„ - \en\-^Sds + {ZoK\-^ - 2|(?„ri/3|0'|2)az„ 
(75) j>n; dg^ ^ -mr'^%\dz. 

Using the fact that afo^ ^^.j. is smooth in Zq and S, on the current region 

Zo,S and are bounded on the region currently considered and the 
computations in (75), we have: 



1, , 4 nf_ . , 2, 



(76) |ai5,"aF0^^^;±| < \0\-3'+-3y 4 

(2) On a; 7^ wc have smoothness in x for %», = X)± (,;±- Together 

with the bound in (76 1, in fact since = x^^Sds we have: 

„ S 5 1 



The order of the symbol can be lowered : S = Z - Zq = x|6l„|2/3. On a; ^ 0, 
for all e N on 

Since |0„| ~ |0„|, |0rG(x|0„|2/3,Zo) G L°°. 

□ 

Lemma 12.2. U^^.o e x''-C;^y + x'+C^y on y„ 7^ 

Proof. Let x(a;, 2/) G with either suppx C {y„ > 0} or suppx C {?/„ < 0} 
(U^o,o,;^) = lim(2^)-" /■ /■ U^o u,^{x,y)w{e6)ey'dxdyde 



lim(27r) "I I I Z^^' aFO,o;±^Xe*^ dxdydO 



± 



\de{y ■ 9)^^ is bounded away from zero on suppx x M" since the n-th component 
of the derivative of the phase function, which is y„, is never zero on this region 
currently considered. Consider a first order differential operator L whose formal 
adjoint is : 

so * Le^y'^ = e^y'^ . fc-fold integration by parts gives 

(U^^,o,;^)^lim^(27r)-"x^±y j J L''{ar<>.^.^w)xe'y'dxdyde 

(1) For k > 3(n+ - A) , 9 ^ L^{aro^^.±w) G L^{W^). We can let e ^ 

under the integral sign : 

4- J ' 4- 



J^c.o;±| < (2^)-" / \L\a^o^ .±)\d9 < Cfc(2^)-" / \e\ ^+iv " ^ d0 

J|e|>i 



The last estimate follows from the fact by Lemma 12.1 



_1 4 nf__ 

iL'^apo .j^\<Ck\0\ 3*^+3 V 4 ^- Cfe independent of x 



'Co 

Thus we have 



(2) Fix k with k > 3{n+\/3\ + - A) so 



Consider D?,\fc,o^^: 



±) 



fl is-o . f I. a Lemma [T2TI /" _lj,,4 /nf 

|^^U^c,o;±|<(2^)-"/|ife(„ < (271)-"/ |0| 3^+3V 4 

ice of k we have 9 i-^ L' 



/|e|>i 

With the current choice of k we have 9 H> -L'^(%o ■±)9^ e L^. Thus 



(3) Fix k with fc > 3(n + | W - A + |0 so that: 



± 

Consider i:'iU^^,o;±: 



From Lemma 12.1 
choice of k we have 



_1, 4 , 2, 

D^L'^apo < l^l 3''+3V 4 ^+3'. With the cmTcnt 



(4) Same argument is used to deal with mixed derivatives DyD\.. In this case 
one needs to choose k satisfying fc > 3(n + \ — A + | / + |a|) so that 

(5) We have shown that U^^.o can be rewritten as J2± 2^"* U^c.o where D^DlU^c "'- 
L'^y for any a and I. Thus: 



□ 



Lemma 12.3. \J c.o {g smooth on x > S for any S > 0. 



Proof. Denote U^'^'" the restrictions of U^^.o to T)'{Xs). For x > S we have : 

(Uf'°,x) = hm(2^)-" / / U^o u,x{^,y)uj{e6)e'yHxdyde 
(79) ' ' ' ' 

= lim(2^)-"/ / I aro^^x{x,y)w{ee)e'y'dxdyde 



where Gbo € S^ for x > 6 hy Lemma 

C.O in 



decreasing in |0| 



12.1 



In fact since on Xs we have rapidly 



uf-'" = (27r)-" / / I a,o,yy <^de 



Note U^'^'° is smooth in y. We can similarly consider dx since dx = x^^Sds and 
SdsttfO^ ^ is still rapidly decreasing in ^ . □ 

12.1.2. On region Fg^g. 

XeCr; (U^^.°,x) = lim(27r)-" / / Uxx{x,y)w{ee)e'y'+'Uxdyde 



= lhri(27r)-" / / I av^ wieO)x{x,y)e'y'+'Uxdyd9 
Lemma 12.4 (Properties of aF°° )• 



(1) Fori such that*L=\deiy9 + 4))\-^J21=i ^''^ '^^ de ■ , we have 
where is independent of x 

90 



(2) on X "^Fc ^ "^i 2 

3'3 



4 /" ^ 



Proof. 



(1) On this region 5 > 1, \Zq\ < 2S2. Choose $2 so that Z = Zo + x\0n\'^^^ > \. 
We have : 



(80) 



where V^c.o jg smooth in Z^^l"^ and Zq. Since a = -(n - 1) < 0, 

la 4 \v?_ 

|Z-4-2V*'c,o| < C|6'„|3V 4 -\ 

We proceed to consider the bounds on higher derivatives: Denote T 
Using the fohowing computation 



j > n; de^ = T%\\9r.\-'/^TdT - m,\-^/^\e,\dzo 
xd^ = -\{T''Zo + l)TdT 

and the fact that T and Zq are bounded on the region currently considered 
we have: 



14 

\{xD^yL^ar^^\<Ck\e\ 3*+3V 4 Cfe independent of X 

(2) On X ^ Q we have smoothness in x of aF~„ and as a resuh of the compu- 
tation in (80 1 we have: 



4 /"^ ^ 

3'3 

□ 

Lemma 12.5. For ?/„ < 0, U'^"^ " is smooth and vanishes to infinite order at x = 0. 
Proof. Consider x{x, y) with suppx C {?/„ < 0} 

(81) (U^^.« , x) = 1^(27^)"" / / / i>r^^u,x{x, y)m{e0)e'y'+'^ dxdyd0 

We have \d0{y ■ 0)\'^ is bounded away from zero on suppx x since the n-th 
component of the derivative of the phase function is y„ — a positive term as follow: 

y„- |a9„Z3/2sgn0„ 

(82) = y„ - 1(1 + X - |0"f )3/2 - (1 + X - l^f )i/22|^f 
— y-n — positive term 
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If Un < 0, the above expression is never zero. On that region we can consider a first 
order differential operator L whose formal adjoint is : 



(83) 



'l = \deiye + 4>)\-'J2 



diyO + ^) 
89, 



thus t\_Qiv-^+i'l> = ^iyS+iip ^ 
fc-fold integration by parts gives 
(U^^.o,x) = lim(2^)- 



(I) ¥oi k > i{n + - \) , e ^ L*(aF-^n7) is in L^{W'). We can let e ^ 

under the integral sign : 



|U^c,o| < (2,r)-" /|L'=(aF^J|d0^Cfe(2^)-" / 

J J\e\>i 



Hence e L°°{x,y) 

(2) Consider D^\fc.o, Fix k with k > 3{r 



4. 



A+l/31) SO 



1 , 4 . / n 



|e|>i 



With the choice of k we have e ^ L'=(aF- )6"5 in in L^. Since the constant 
is independent of a; 



(3) Consider xD^U^c.a, Fix fc with fc > 8(71 + - A + I) 



'CO' "CO' ■^^ 

Consider the second term in the bracket: xdx^Z^^'^ — x\0n\il+x— \9'\'^)^^^. 
We have \9'\ < (1 + ki{S2)) by (851 and x < C, so for C is independent of 
x: 

\xdJ\<C\9,,\ 



From Lemma 
we have: 



12.4 



xD.,L''{ar^J<\9\ 



1 , , 4 71" , , - 



With the chosen k 
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For higher derivative (xl?^:) U '^■o, we simply choose k such that k > 
3{n+ 1^^^ - X + N) so 



(4) Consider a;-^U''"o 

= (z - Zo)-i|0„p/3 = t2(i - ZoT')-'\e^\'/' 

since |Zo| < 252 and < T = Z^^ < 2. So for Cat is independent of x: 



Fix k with fc > 3(n + iV + - A) so 



J ■/|e|>i 
with the choice of k then have 

The same reasoning apply if we replace U^'^-o by {xD^YiJ^'^ o . 
(5) Consider x'^ D^{xD^ fU^c,o = L)^x-^(a;L>^)^U'^c,o. Derivatives in y in- 
troduce growth in \9\, while is polynomial bounded in 9. As above, 
by the same integration by parts technique on j/„ < we show that we can 
arbitrarily reduce the growth in |0| thus 

□ 

Lemma 12.6. [On y„ > / 

For every /3 > 0, there exists a > such that [j^c.o ig smooth and vanishes to 
infinite order on yn > l3 and x < . In fact a is given by {84) explicitly. 

Proof. For a > whose relation is to be determined, consider xiv) with supp x C 
{yn > /3,x„ < a^}. 

{X < a'^ 
Zq + x|6'„p/^ > 1 
\Zo\ < 262 

we have ( can assume S2 < \) 

xK\'/' = l-Zo > 1-\Z,\ > 1-252 ^ K\'/' > ^ \0n\ > C-^X^' > (^) 
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1 



Q'|2 I 



\Zo\ 



12/3 



< 1 - \9'Y + x <2a^ (since x < a^) 



Using this inequality and 

< 1 - KoiS2) < \0'\ < 1 + Ki(52) < 2 ^ l^'l < 2 

we obtain a lower bound for: 
Now choose a such that 

^23/2^3 + 12 X 2i/2a < ^ 
Thus on supp(x(a;, y)( we have 

ie the above expression is not zero. 
Consider C with 



(84) 



3 = 1 



do. 



By Lemma [12. 4| we also have 



4 n± ,1, 



From this point, the proof uses the same argument as those of Lemma 12.5 to show 
on the region currently considered of x the boundedness of 



□ 



Lemma 12.7 (Statement about wavefront set on x > 0). The wavefront set of 
ifc.o ig contained in ^,{62) where 

S(52,<5) = {ix,y,C,ri) : ^ (y ■ ~ ^'j ,r] dy (^y ■ 6 - cf^ ,d0 (^y ■ - (f^ = 0,6 e conesuppaw^^,x > S} 
with 4> = \Z'^I'^ sgn^n and 



conesupp aw^ 



> 0, (1 - K„{d2)) \0n\ < \0'\ < (1 + Ki{S2)) \dn\}. 



1^0(^2), ^1(62) are defined in the calculation of conesupp ap^j, o,t the end of the 
subsection. 



3 V 4 ' 

Proof. flFj^Q S I by Lemma 

^ ,0 



12.4 



In addition, since we have the same type 



of symbol and same phase function as in Friedlander Theorem 6.1 [2|, we employ 
the same argument to prove our results. The main idea is that on a; > 5, $ is 
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smooth on conesupp aw-^g and is homogeneous degree 1 with respect to 9, thus is 



a phase function. 

(Uf ^x) = lini(27r)-"y J j a^^^x{x,y)^{ee)e'^ dxdydQ 

is an oscihatory integral in the sense of Hormander. The wavefront set statement 
now follows from proposition 2.5.7 in FIO I by Hormander [S]. 



□ 



Calculation of conic neighborhood 

(85) {\e\ > 1, \Zq\ < 2S2} c m > 0, 1 - ^^0(62) < \e'\ < 1 + ^1(^2)} 

for some positive constants kq and ki depending on 62- 
Proof. Consider the set {|6'| > 1, |Zo| < 2(52} 



(1) Consider the case in which Zn > ic < 2 1 — 

^m~'7r^-^<2S2 



9'\ < 1 

1 

(T+a2)i/3 



Denote A = {I + a2)i/3 go ^ ^^3 _ 



(l + a2)i/3 ^ V ^ ^2 

0<a<lsol<A< 2^/3 . On this region f'{A) < 0, f{A) is monotone 
decreasing from 1 to 0. Hence for 62 such that < 2^2 < 1 , there exists a 
unique A in this interval such that 

f{A) - 2<52 

and since a > there exists a unique kq{S2) such that with a = 1 — ko(<52) 

(l + «2)l/3-2'52 

This is the value for |^'| for |6'| = 1. 

1 _ a2 

If 101 increases from 1, want ^wtt^ to decrease, which happens if a 

(l + a2)i/3 

increases. So we have: 

{\e\ > 1,0 < l^'l <l,Zo< 252} C {\9\ > l,l + Ko('52) < l^'l < 1} 

Note that 1 — ko{S2) < 1 so ko((52) > and as (^2 1 + ko(^2) 1 ie 
KoiSo) 

■\9'\>1 

(2) Similarly we consider the case for Zq < ie 2 — 1 



Denote A = (1 + a2)i/3 ^ = ^3 _ 



(l + a2)i/3 



A^) - (1 + ^2)1/3 - ^ / - 
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o > 1 so A > 2^/^ hence f'{A) > 0. f{A) is monotone, increases from zero. 
There exists a unique A so that 

f{A) = 252 

and thus there is a unique 1 + ki(^2) so that for a = 1 + «;i(^2) we have 



(l + a2)V3 



2S' 



2 



1 - 

If 1^1 increases from 1, want ,„ to decrease, which happens if A 

' ' (1 + 0^)1/3 

decreases so 

{1^1 > 1, l^'l > 1, < -Zo < 2^2} C m > 1, 1 < l^'l < 1 + Kl(<52)} 

Note that ki{S2) > and as (52 ^ \ 2^ ie a \ 1 ie ki{62) \ 
In short we have: 

{|^| > 1, \Zo\ < 252} C {1^1 > 0, 1 - Ko{So) < \e'\ < 1 + Ki{52)} 

Remark 12.2. 

• The condition Zq < 252 and \6'\ < [1 — 25i)^^ does not imply that Zq is 

bounded below ie does not imply that \Zq\ is bounded. 

• The difference between region currently considered and coneregion cur- 
rently considered : in the region currently considered , with constraint 
< < 252 (correspondingly < —Zq < 2S2) , when |6'„| increases, 
this forces ^-^^~2°i/3 (correspondingly (i"„7-)V3 to decrease) , while in the 
coneregion currently considered, we do not have this constraint, hence \Zq\ 
will not be necessarily bounded. 

• not the same situation with Zq > 52 or —Zq > 52- 

□ 



12.2. On region ¥c,+- On the current region we have {\9\ > l,Zo > 62}, whose 
cone region currently considered is {\6\ > 0,\0n\ > On region Fp_,_ we have 

0, while we subpartition F^ _|_ into ^^0!+ (corresponds to /U ~ 0) and F^°^ 
(corresponds to ^ 0), where /i is the blow-up variable according to the blow-up 

at the corner of face ¥c at cr = and h = . 

For subregion F^°j_ we show 

(1) U + is smooth on a; > 5 for any ^ > ( and any 52 > 0). 

(2) U^c"+ e x'-C^y + x'+C^y on y„ ^ 0. 

For region and F^°^ we will show 

(1) On y„ < 0, 1)"^^,+ ^ U'^"^.+ are smooth and vanishes to infinite order at x = 0. 

(2) On Un > U , U c,+ are smooth and vanishes to infinite order if we are 
close enough to the boundary for y„ > 0. 
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(3) wavcfront set statement : contained in the forward bicharacteristics from 

X — y = 0. 



12.2.1. OnF^°^_. 

(86) (Uf'+,x) = lim(2^)-"y J j i>^o,o^u.x{x,y)vj{ee)e'yUx dy dO 

We recall the properties of the approximate solution on the region currently con- 
sidered: 

V,,o^o^^, = (|0„|(1 - |0f)V2)/?^O^c"+ 

where 0'^'^,+ is conormal in t{— f^) and smooth in h, where t — x\9n\{'^ — 
and h = |6'„|-i(l - |0'|2)-3/2. in more details: 



on n , n n . 

2a/ - A ^ Z : U^c,+ e a; 2 iV 4 -^c^^ + 3. 2 TV 4 -^cj^^ 



(l-|0f)t[|0„|(l-|^f )i/2]Vx-^U^^'.+ e [|0„|(l-|^"f )i/2]^Vx-\«+c,^^+a.-c,~^ 



Using this form of Ux we rewrite ( 86 ) as 

BfO.O 



(U,^-+,x) = limE(2^)^"^'*y 7 j a^o^o^.^^x{x,y)w{ee)e'y'dxdyd9 
Lemma 12.8. [Properties oj a^o.a ._^] 

1 d 

(1) Lef L be the differential operator whose formal adjoint is *L — -j — — ^ yj 
then for Ck is independent of x 



(2) Onx^O, 



\D'^L''{a^o.o )| < Cfcl^fs^+'V 4 ^+'" 



Proof. 

(1) First we consider a bound on 0^.0,0 ._|_. 



V^o^^. = (|0„l(i - |^f)V2) V --^u^c.'V 



(l-|0f)t[|0„|(l-|^f )i/2]V 4 -^u^c,+ e [|0„|(i_|0"f)i/2]2V 4 -V+C,^^+a;^-C,^^ 
In the region currently considered on which h,t and \6'\ are bounded 
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2\/x-^ 



To compute the bound for higher derivatives, we use: 

\0n\de^ = (l - 2|^f /^2/3|0^pl/3^ ^ (-\e.n\-' - S/l^/S 1 0„ | - 1/3 | ^' 1 2^ hdh 

xd.,^tdud^ = \6r,\{i-\e'\'Y'^dt 

Using the above computation, the fact that a^o,o , is smooth in t, h and 
on the current support h, \6'\ are bounded, we obtain: 

(87) \D^L\a^o^.^.,^)\<Cu\en' 

(2) On a; 7^ we have smoothness in x for a^o,o — ^_|_ x^^a^o.o Together 



with the bound in (87 1, in fact since dx = x ^tdt we have: 

2\ 

(XpO.o S I 

In fact, we can improve the order of the symbol in this case with we are away 
{x = 0} : in this region currently considered since Z > 62 ie |0„| > \9'\ 

so \6n\ is equivalent to \6\ thus it suffices to compare with |0„|. Since 

_ 3 

\9^\ = X 2i3/2/ii/2^ on X > S for aU N we have : 



N 



so flpO.o G S I 



0:9 



□ 



Lemma 12.9. U^'^ + is smooth on x > 5 for any S > 0. 



Proof. Denote by the restrictions of U to 'D'{Xs). For a; > 5 we have : 



(Uf' + ,x) = lim(2^) 



lim(27r)-" 



tpfO.o UxX{^Ty)™{^^)^^^^ dx dy dO 
a^o,o xix,y)vj{e6)e^^'^ dx dy d9 



From lemma (12.8) on x > 5 we have a^o.o g 5"^°° hence U^'^'+ , L'^U_5^ + e L° 

77,0 



The same situation remains when we consider dxUg'^'*. Since dx ~ x ^tdt, on 
X > 6, tdt{apo,o ) still rapidly decreases in 9. 

□ 

Lemma 12.10. U^c.+ g x^-C;^^ + a;^+C'*+ on y„ ^ 0. 
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Proof. The proof is the same that for Lemma 12.2 Consider xiXiV) G ^(xy) with 
either suppx C {un < 0} or suppx C {un > 0} 



(U^c' + ,x) = lini(27r)- 



ip-po.o Uxx{x^y)w{e9)e'^^'^ dxdydO 



= lim(27r)"" y" j j ^x'^ a^o.o^.^vD{ee)xe^y''^ dxdydO 



\d0(y-O)\'^ is bounded away from zero on region currently considered of x x Hi" since 
the n-th component of the phase function, which is 2/„ , is never zero on this region 
currently considered. Consider a first order differential operator L whose formal 
adjoint is : 



Do fc-fold integration by parts to give 



(U^c°+,x) = limV(2' 



L''{a^o,o .^vo)xe''^'^ dxdyd9 



(1) For k> ■i{n + 2'^'^ - X), 6 ^ L^^ ap-^ G Li(M") by Lemma 
can let e — > under the integral sign in ( 88 1 to obtain: 

-I- J •+ 1 



12.8 



We 



By Lemma 
of k we have 



12.8 



■|pO,0 _i_ 



With the current choice 



(2) Fix k with k > 3{n + \l3\ + 2^ - A) in (88 1. Let e ^ under the integral 
sign to obtain: 

J- ^ ' 4- 



Consider D^U^c.+''^. 



de 



/|e|>i 



With the current choice oi k, 6 i-^ L'^{a^o,o )6"3 e by lemma flZi 



we have: 



so 



(3) Fix k with fc > 3(?i + 2J^^ - A + /) in (88) 
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Consider Dl{}'^o.+ 



1 TC'J ' " 



(27r) 



From lemma (|12.8|), D^.L'' 



--A+i 



With the current 



choice of k, we have 



(4) Same argument is used to deal with mixed derivatives DyD\.. In this case 
one needs to choose k satisfying fc > 3(n + 2\J ^-j — A + / + |a|) so that 



y X- 
5.0,0 



^ CO 



(5) We have shown that U can be rewritten as ^^.x^^U 0.+' where 
D^Dl\J^c°+'^'^ e L^^y for any a and I. This means: 



2;, a 



□ 



12.2.2. For subregionF'^^^. 



e-i-0 



(U*c,+ ,;^)^lim(2^)-" / / i^^^u^x{x,y)vj{e9)e'yUxdyde 



lim(27r)-"y" ^ ^ aw^ ^vjxe'y^+"'' dx dy d9 



We list the properties of ap^^ that will be of later use: 
Lemma 12.11. 

(1) Lw-^^ whose formal adjoint is *iF~^ = \deiy ■ 9 + (j}))]"^ J2'j=i ^^MJ^'^^j' 
(j) = - Zq/^) sgn6l„, then: 



(2) On X > 6, apoo ^ 



3,0;e 



Proof. 

(1) Recall the structure of the approximate solution from the construction 



U*'c.+ = exp(i(/))/i'''^i + 5o--^vi'c,+ 
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where V*'c,+ ^ h^Aj{a-^/^) with Aj smooth, c7"2 = a;"4(l - \e'\^)l. 
On the region currently considered h, and |^'| are bounded, a < , and 
< 1 - l^'P < 1 we have 



For estimate on the higher derivatives : from direct computation we have 

Since h,a~^ are bounded (independent of x) we have: with Ck is indepen- 
dent of X, 

{xD^yi^ap^J < Ck\e\-~3''+^^\ 

(2) From the same computation, and on x ^ dx = —\x~^a~^/^d^-i/6 we 
have 



3' 

□ 



Lemma 12.12. Foryn < U «s smooth and vanishes to infinite order atx = 0. 

Proof. Consider x{x, y) with supp x C {j/n < 0} 

(U^^.+ ,x) = lim(27r)-"y y j ikv^^^Uxx{x,y)w{ee)e'yUxdyde 

= lim(27r)-" y ^ j a^^ ^ruxe'^'^^'''' dxdydO 

We have • is bounded away from zero on region currently considered of 

X X R" since the n-th component of the phase function is: — a negative term. 
This is seen as follows: 

2/„-|5,„((z3/2-Zo^/2)sgn^„) 

= yn-l[{l + x- l^f - (1 - l^'f - \e^\ [(1 + ar - )V22|^f - (1 - )V22|^'|2|0„|-i 
= Vn — positive term 



If y„ < then the above expression is never zero. Denote cf) = y-6 + (f). We consider 
a first order differential operator Lf«j whose formal adjoint is: 



do, 

thus Lv^^e'y^+''^ = e*y«+*<^. Do fc-fold integration by parts to give: 
(89) (U^~.+ , x) = 1^0(27^)-" j j J i^^^ iav^^w)x-'2^xe'y'dxdyd9 
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(1) Consider fc > 3(n + ^ - A) then 9 n- L^^ ar-^^ G L^{W^). We can let 



e — >■ under the integral sign in ( 89 ) to obtain 



|e|>i 



\e\-l^dd 



(2) Consider i:'^U*'"+. Fix k with fc > 3(7i + - A + so 



dO 



|i5/^U^c,+ | < (27r)-" f\L^^^{av^je^\de<{27r)-'' f |0r5'=+ V 4 -^+1^1 
With the choice of /c, 6* i— > L^oo (aF^^)^'^ G . Since the constant is 



independent of x 



(3) Consider a:L':,U*c,+ . Fix with k > S{n + \ \ - X + 1) 



xi?^U*c.+ = (27r)-" / xDa;L^^^aw^^ + {L^ar^^)xD^(j) 



Jy-9+i 



Since |^'| < 1 and x < C we have Ixc^^^^j < C\6n\ where C is independent 
of X. This together with the bound on xD^L'^av^^ < Ck\Of^''^^ 



given by Lemma 12.11 and the current choice k, we have 
Similarly reasoning applies for {xD,j.)^\J^c:.+ . In that case we need to choose 



k such that fc > 3(n + a/ ^ - A + TV) so 



{xD^fL^'iar^] < Cn\0\ ^ 



{xD,)^L''{ar^)eL°°{x,y) 



^ c,+ ' 
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(4) Consider x-^\J^'^.+ : We have 

Since h < 8^ and < cr^ < 1 we have for Cat is independent of x 



Fix fc with fc > 3(n + ^ - A + |7V) . 



with the choice of k then have 

The same argument goes through if we replace U^"^'+ by Dy{xDx)^yJ^'^'+ ■ 
So in fact from here we have 

x-^D'^{xD,fU^o.+ ^L^^^ 

□ 

Lemma 12.13 (Properties of \J^'^'+ on ?/„ > 0). 

For every /3 > 0, there exists a > such that U'^'^'+ is smooth and vanishes to 
infinite order at x = on {?/„ > l3, x < a^}. In fact a is given by (91) explicitly. 



Proof. 



Denote 



(U^~+,x) = lim(27r)-"y j j i,r^^u,x{x,y)vj{ee)e'y'> dxdydO 
= hm(27r)-" I I I x{x.y)vjUe)e'y''+'* dxdydO 

e— >-0 III • + 



So = {{x^y) : dg{y ■ 9 + (p) — 0, for some 9 e conesupp a^-^ ^} 



For /3 > 0, consider x{x,y) with suppx C {y„ > P,x < a^} where the relation of 
a and /? is to be determined and such that suppx is disjoint from Sg. We show 
below that if a satisfies (91 1 then we have is bounded away from zero on 

region currently considered of x x M". Once this is achieved, we follow the same 
argument as those for j/„ < in Lemma 12.12 to show that on the region currently 
considered of x^ U'^c-.+ is smooth and vanishes to infinite order at a; 0. 
In particular, on the support of x (described above), we consider the first order 
differential operator L whose formal adjoint is : 

(90) %^^^ = |d,(y.0 + 0))|-2^^i^^^±ila,,^ L^^^^e^y'+^^^e^y'^+^'^ 
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□ 



It remains to show that: for /3 > if Un > f3 and if we are close enough the 
boundary ie x < for some a, then dg{y ■ 9 + (p) ^ : 



Proof. 

2/„-|ae„((Z3/2-Zo^/2)sgn0„) 



Vn 



(l+X- l^f )3/2 _ (1 _ l^f )3/2l _ [(1 + a; _ |^/|2)1/22|^.|2 _ _ |2 ) 1/2 3 |2 



Consider the region {a; < a^}, where a is to be determined. 
Denote by G(x) the foUowing expression: 

G{x) ^l\{l + x- - (1 - l^f )-V2l + + l^f )i/22|^f _ (1 - \§'\Y/^2\e'\ 



On the region currently considered we have a — x^^^{l — \6'\^) ■^/^ > 1 so 1— < 
X. Since x < c? , < 1 — \&\^ < o? . Use this to estimate |G(a;)|: 



\G{x)\ < -a3(23/2 + 1) + 2a(2i/2 + i) 



Choose a < 1 such that 



(91) 



-a3(23/2 + l) + 2a(2i/2 + l) < 1/3 
3 2 

^t/-G(a;)>y-|G(x)|>0 



□ 



As a corollary we have : 

Lemma 12.14. Singsupp U'^c,+ contained in the union of the forward bicharac- 
teristic from x ~ y ~ 

Lemma 12.15 (Wavefront set statement). The wavefront set o/U^c,+ is contained 
in S((5) , where 

^(6) = {ix,y,^,7j) ■.^^dAyO-^),V^dy iy9-<l>),de (y ■ 9 = 0,x > 6,\0n\ > \e'\ > 

<^=|(Z3/2-4/2)sgn0„. 



Proof. By Lemma 



12.11 



on X > S, (i <^' + ^ S I . Compared to Lemma 5.2 

in Friedlander paper, we have the same type of symbol, the same phase function 
on the same coneregion currently considered, as well as the same result on the 
singular region currently considered of U'^C'+ . Hence we have the same conclusion 
onWFU'^c,+ . □ 
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12.2.3. On region ¥°(j'^. 



(U*c,+ ,;^) = lim(27r)-" / / i^^,^u,x{x,y)vo{ee)e'y'dxdyde 



= lim(27r)-"y j j a^»,^x{x,y)vD{ee)e'y^+''>' dxdyde 



Lemma 12.16. 

1 , 2 



(1) Lj.0,00 whose formal adjoint is L^o.oa = \d0{y ■ +(f))\ ^ X^J^i dg. then 



(a;£>,)'L^o,ooapo,oo <(7fc|^|-3'=-3' 



Fc,+ "'0,+ 
_2 

(2) On X > (5, a„o,oo e S", 3"*". 
Proo/. 

(1) Recall the structure of the approximate solutions: 



^ -iniO , OO - . 1 _ ^ .-,0 , OO 

\fc,+ = exp(i^)z/^'«i + 2£7'^'»i V'^c.+ ; 7<si = 



where V c.+ is conormal in z and s(= ^). In fact it is polyhomogeneous 
with index set {j}. On the region currently considered z and s are bounded. 
On the region currently being considered, 



1^? 



To compute the bound on higher derivatives, we use the following com- 

^ Tn.0,00 

putation and the conormality of V '^•+ . 



1 + s3z^\en\^\e'v j sds + 2s3z3\en\3.\e'Yzd^ 

\e^\de, =-s^zl\e^\l\e'\sd, + 2slzlej\en\^zd^ 
d^ = s^'yr.\{i-\eT)^{-sds + zd,) 



to conclude that 



(2) We use the same computation and the fact that dx = —x~^sds + x~^zdz, 
to conclude that on a; > 5: 



□ 
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Consider x smooth in x, y with supp x C {yn < 0} 

(U^o,+ ,;^) = lini(27r)-" / / [ i;j,o.^u^x{x,y)nj{ee)e'y^dxdyde 



lim(27r)-"y" J J a^o.a.x{x,y)w{ee)e'y''+''^ dxdydO 



Since we have the exact same phase function, same result on the properties of 0^.0,00 



(compared to ofj?^) and 

x-^ = s|6l„|(l - |^'p)i/2 < constant |6'„| 

we use the same method as in Lemma |12.12[ |12.13[ |12.14| and |12.15| to get the 
following results: 

Lemma 12.17. On j/„ < U c' + is smooth and vanishes to infinite order at x = 0. 

m0,OO 

Lemma 12.18. Sing supp U c,+ jg contained in the union of the forward bichar- 
acteristic from x = y = 0. 

Lemma 12.19. For every /3 > 0, there exists a > such that U is smooth and 
vanishes to infinite order on y^ > (3 and x < o? . This means that we have smooth- 
ness for y-n > if we are close enough to the boundary, ie we have conormality for 
yn > 0. 

Lemma 12.20. The wavefront set of U + is contained in where 

S(<5) ^ {ix,y,^,7j) : ^ ^ {y ■ - ^) ,71 ^ dy {y e - <j>) ,dg iy ■ e ~ (j>) ^ 0,x > 6, \0, 

0= |(Z3/2-Zo'/')sgn0„. 

12.3. On region Fc,-a. On the current region {Zq < -62, \9n\ > (1 - 2(5i)|6i'|} 
so \9n\ ^ \d\, hence it suffices to consider and compare with the growth of |6'„|. We 
will show: 

(1) On F^°_ ^: U^c°-.^ is smooth on x > 6 for any (5 > while U^c°-,i g 
x^-C-^+x^+C- ony„^0. 

(2) On F-_ 1 

(a) U^<^'-'i is smooth and vanishes to infinite order at a; = for y„ < 0. 

(b) smooth and vanishes to infinite order at a; = for y„ > if we are 
close enough to the boundary. 

(c) Wavefront set statement away from x — . 

(3) On F^°! U'^cJ.-.i is smooth on x > S for any S > 0. U'^C.-.i is smooth 
and vanishes to infinite order at a; = for y„ 7^ 0. 

12.3.1. Subregion¥^'^_ ^. 



(92) (U''°c^'.-.i,x) = lim(27r)-" / / j iPj.o.o U'''^.-^xix,y)T:u{e9)e'y'' dx dy d0 
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We recall the structure of Ux on the region currently considered: 



where (fc,-,i jg smooth in t and smooth in h. Here t — x\9n\{\6'\'^ — 
l^nT^d^'P - lY'^- In more details: 



[|^„|(|0f - 1)1/2] V --^QF^d^,, g [|0„|(|^f ^ i)i/2]^V T~ V+C,-^ + x^-C: 



t.h 



Using this form of we rewrite ( 92 ) as 



(U^c.-,i,;^) = lim^(27r)-"a;"± j j j a^o.o^ ^.^x{x,y)^{e0)e'y^ dx dy dO 

Lemma 12.21. 

(1) L whose formal adjoint is *_L — \y\^'^ 



(2) On X > 6, Opo.o = ^_|_ x^^Opo.o e 5*]^°° 



Proo/. 

(1) First we consider the bound on 0^.0.0 



[l^nKl^f - 1)'/'] V ^~^O^c°-.i e [l^nKl^f - i)i/2]'V V+c,^^ + x'-q 

In the region currently considered on which h,t and \0'\ are bounded and 
since s_ > 0, so 



To compute the bound for higher derivatives we use the following compu- 
tations: 

Klde,^ = (-1 - 3(-Zo)-i|0„p/3|^'|2^ hd,, + (1 - (-Zo)-l|0„p/3|^.|2^ 

\9r^\de^ = 3{-Zor'\e^\-^/^\9'\hdh + i-Zo)-^\e.^\-^/%\tdt 
xdx = tdt; a, = |0„|(|0f- 1)1/29* 

Using the fact that 0^0,0 ._|_ is smooth in t, h, on the current region t, h 
and l^'l are bounded, and the computations above to get: 
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(2) We have 0-2(^0), o'o(l^l) G Sl . From the prep calculation, on a: 7^ 0, we 

3,0 

have 



_ 3 

In fact we can improve the order of the symbol |0„| = x 2 t^/'^h^/'^ . Thus 
on X > (5 for all N 

From the previous remark and using the fact that |0„| ~ on x > 5 
we can lower the order of the symbol ie € 57°° . 

□ 

Lemma 12.22. U c.-.i is smooth on x > 5 for any S > 0. 



Proof. Denote U^j*^' the restrictions of U'^c,-.i to ^'{Xs) .For a; > 5 we have 
(U?'-'\x) = lim(2^)-" / / / uMx,y)w{ee)^yUxdyd9 



lini(27r) "J J j a^o,o_x{x,y)w{ee)i 



e'y^dxdydO 



where € 5*1 °° for x > 5. The proof is the same as that for Lemma 

Lemma 12.23. [fc'^-.i ^ x^-C^ + x^+C^ on y^ ^ 



12.9 



□ 



Proof. Same proof as that Lemma 12.10 □ 



12.3.2. RegionV 



c,-,i- 



On the region currently considered we have: cr > 1 (so Z > 0), < —82, 1 < 9' < 
(l-2(5i)-i 



QFc,-,! = exp{-ilZ^/^sgnen)exp{-l{-Zof^^)h^<T-^\/^J- '; for any iV 
where V'^c-.i is smooth in cr~^/'^ and /i. Denote by = ^Z'^^^ sgn6„. 

(U^c,-,i,;^)^lim(27r)-" / / l^r^_ u,,x{x,y)uj{ee)e'y'dxdyde 



lim(27r)-" / / j a^^ _ ^x{x,y)Tu{ee)e'y''+'^ dx dy d9 



Lemma 12.24. 
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(1) For Lf^ whose formal adjoint is* Lf^ = \de(y-9+(i>))\ ^ y]"_i ^''^aa^'^'' dg 
we have 



(2) exp(-|(-Zo)3/2) e 5*; aTidapo.oc € S'l 



— X 
4 



3,0;e 



Proof. On the region currently considered exp(— 1(— Zo)'^/^) is a smooth function 
of (— Zq)"^/^ whose derivatives are bounded. Using the same argument as in Fried- 
lander lemma 5.1, we have exp (— 1(— Zq)"^/^) G S\ . The proof of the statement 

follows from direct computation similar to those in Lemma |12.11[ 

□ 

Lemma 12.25. U^'^'-'i is smooth and vanishes to infinite order at x = for 
Vn < 0. 



Proof. The proof is similar to that for Lemma 12.12 The first ingredient we need 
is the bound of in terms of \9\ using x~^ = a~^h?'l^\QrS^^'^ . The second is the 
fact that the n-th component of the derivative of the phase function is j/„— positive 
term as in Lemma [T2T2) 

2/„- |9e„^^/'sgn0„ 

= y„ - 1(1 + X - \B'ffl^ - (1 + a; - \B'fYl^2\B'fd^ sgn0„ 
= Un — positive term 

□ 

Lemma 12.26. For every /? > 0, there exists a > such that U^"^ - ! is .smooth 



and vanishes to infinite order on y„ > P and x < o? . In fact a is given by (93) 
explicitly. 

Proof. 

(U^~.-.Sx) = lmi(2vr)-"y J J u,x{x,y)vo{e0)e'y' dx dy dO 

= lim(27r)-" y" J j a^^ _ ^x{x,y)vj{ee)e^y'^+'^ dx dy dB 

Denote by 

So = {{x, y) ■ dg{y ■ 9 + (j)) ~ 0, for some 9 G conesupp a^^ ^ ^} 

For /3 > 0, consider xi^^ y) with supp x ^ {yn > P,x < a^} where the relation of a 
and /3 is to be determined and such that supp x is disjoint from Eo- Wc show below 
that if a satisfies (93) then we have \dg{y ■ 9 + is bounded away from zero on 



region currently considered of x x M". Once this is achieved, we follow the same 
argument as those in Lemma [l2.12| to show that on the region currently considered 
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of X, \J^c,-,i is smooth. In particular, on the support of x, wc can consider the first 
order differential operator Z/Fg>_ ^ whose formal adjoint is: 

= \de{y ■ e + 4>)r E ^^""gl^^^ de, ^ iF~_,e^^-^+'^ = e^^-^+^^ 

thus Lf- e*2'-''+*'^ = f,iye+i4>_ 

(',-,1 

Show that for j/„ > , if we are close enough the boundary then de{y-0 + 
Proof. 

1 < |^'2| < _ 2J^)-2 ^ 1 _ (1 _ 25^)-2 < 1 _ |^/2|2 < 

^ X + 1 - (1 - 26i)-^ < a; + 1 - |^'^|^ < a; 

Since cr > 1, xd^'p - > 1 so x + 1 - \0'^ >0^0<a; + l- l^'p < x. 
Using this, we obtain a lower bound for: 

de^iye + 4>) = yn-Ul + x- l^f - 3(1 + X - l^f )i/2|^f sgn0„ 

> yn - Ix^/' - 3xV2(i - >yn- - 3a(l - 25,)-' 

Choose a so that 

(93) ^a3^3a(l-2(5i)-2 < 

□ 
□ 

Lemma 12.27. Wavefront set ofD^'^---^ is contained in 

S((5i,(5) = |(x,j/,^,7?) :^ = dx{y-e-^,ri = dy{y-e-(i)^,de{y-9-^=0,e& conesuppav^_^,x > 6^ 

where 4> = ^Z^/'^ sgn^„ and 

conesupp aF~_^ = {|^| > 0, (1 - 2Si)\e'\ < \6n\ < \e'\} 

Proof. The proof follows the same as that for Lemma 6.1 in Priedlander □ 

12.3.3. Region F^°^ ^. Recall the structure of the approximate solution on this 
region : 



V,o,.^«. = (|^„|(|^'p- 1)1/2) Vx-o^o.-.> 

U""?^,-.! =exp(.^)si/2^^'SiU"'2^.-.i; 758^ =0 
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where (j) = - '^i' ^^^cr [l - (1 - ct^/^)^/^] and 0*"c=.-.i is smooth in z and s. 
Lemma 12.28. 

g^|fl„|(|8r-i)'''^(i-[i-^]^^^) < g-|e„|x3/=s(=.) 
with g >0, x^/'^g{x) ^ when x 0. 

Proof. We look for a lower bound for 1 — (1 — 2)^/^. We have 

{\9f -!)-'> [{1-25,)-'-!]-'^' 
Plugging this into the expression for a we get 

z = x{\e'\^ -l)>x[{l- 2(5i)-2 _ 1] -3/2 

-1 3/2 

^ 1 - [1 - zf^^ > 1 - l-x[{l- 26i)-^ - 1] = g{x) > 
Note that ^(a:) when x j^O. 

In the region currently being considered, we have z < 1, so we get a lower bound 
for - 1): 

^3/2(|^/|2 _ i)-3/2 ^ 1 ^ ^3/2 < (|^/|2 _ ^^3/2 

Using these two lower bounds we obtain: 

with x3/2g(a;) ^ when x ^0 □ 
Lemma 12.29. 

1 d 

(1) For L such that *L = -p-nr y^Vi^^ we /iat;e; 



(2) On a; ^ 0, a^o.oo G 



C,-,l 

-C30 



Proof. 



On the region currently considered > (1 — 2(5i)|^'| so |^„| ^ |^|. We 
have : 



%o,-^^ < |6i„|V 4 
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Similarly to the proof of Lemma [12.21| we have: 



1 , 2 



From the Lemma 12.28 on a; ^ we have rapidly decreasing in order of \6„ 
given by e V l j / so tor all A* we have: 

\On\^a^o.^ ^ < Cat 
and the same holds for derivatives in 9j and dx- Hence on x ^ 



□ 



Lemma 12.30. U c.-.i smooth on x ^ 0. 



Proof. Denote U^'^'"'^ the restrictions of U ^ - i to 'D'{Xs). For x > 6 we have : 



(U,--\x) = lim(2^)-" 



7/;jpo,o Uxx{x,y)'uj{e9)e^^'^ dx dy d6 



lim(27r)-"y" ^ j a^o,^_ x{x,y)w{ee)e'y^ dx dy dO 



where 0^.0= € S*, °° for a; > 5 by Lemma 
as that for Lemma [12.91 



12.29 



The rest of the proof is the same 

□ 



Lemma 12,31. U ^-'-'^ is smooth and vanishes to infinite order at x = for 



Proof. Consider x{x, y) with either suppx C {?;„ < 0} or suppx C {?/„ > 0} 
.,^_l:™/o_^-" ; ; / ^^^^^ Uxx{x,y)vj{ee)e'y" dx dy dO 



(94) 



lim(27r)-" 



a„o,oo vjxe^^^dxdyd9 



We have \d0{y ■ 9)\'^ is bounded away from zero on the region currently considered 
of X X I^" since the n-th component of the phase function is which is never zero 
there. 

We use the first order differential operator L whose formal adjoint is 



(95) 



d 

we. 



thus ^Le^^'^ — e^y'^ , and is used for the /c-fold integration as in Lemma 
(U^c.-.i,;^) = lim(27r)-"y" J j L^{a^a^^_w)xe'y'^ dxdydO 

To show rapid decay in x, we use the boundedness of x~^ as follows 
X-' = s|0„|(|^f - 1)^/' < s\e^\ [(1 - 25^)-' - 1] 
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12.10 



The rest of the proof follows as in that for Lemma |12.10[ on y„ < to show the 
boundedness for 

□ 



12.4. On region Fc,-,2- On the current region {Zq < -S2, \9n\ < (1 - Si)\0'\}. 
Since \e'\ > (1 - Si)-\ z ^ x{\e'\^ - 1)-^ < x {{I - 61)^ - l)"^ ie bomrded. 

We will show 

(1) U^c"-,2 is smooth onx^O and U^c"-.2 ^ x'-C^y + x'+C^y on ?/„ ^ 0. 

(2) U '^.-.2 is smooth on x 7^ and is smooth and vanishes to infinite order at 
x = for y„ ^ 0. 

12.4.1. F^°_ 2- 

(96) (U^°^'.-.^x) = linr(2^)-"y" J ^ Vf«^o_ u,x{x,v)w(ee)e'yUx dy dO 
We recall the structure of iix on the region currently considered: 



where U <^---^ is conormal in t and smooth in /i . In more details: 



[|0„|(|^f ~ 1)1/2] /¥^U^C."-^ ^ [(|^,|2 _ 1^^1)1/2] 2V^^,^^«3^ ^ ^,_^oo 



Using this form of we rewrite ( 96 ) as 



(U^c,-,2^^) ^ lim^(27r)-"2:'*± ^ ^ j aF;;;° ^,±wxe"^'^ dx dy dO 
Lemma 12.32. 

1 d 

(1) Let L he the differential operator whose formal adjoint is *_L — - — — ^ yj - 



then for Ck is independent of x 



(2) Denote by — ^_|_ x'^^a^o.o On x > S, for all N and all a 



m'^xDj^d^a^o.o 



< C 



N,m,a 



Proof. 
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(1) First we consider the bound on a^o.o 



[K\{\0'\' ' 1)1/2] V^O^^.''-.. ^ [(|^,|2 „ 1^^1)1/2] 2V^^.+ ^oo^ ^ ^._^oo 

On the region currently considered, h,t and |0„| are bounded so 



[|0„l(l^f- 1)1/2] vx-^u^c,-,. < |0'|2Vx-^ 



For bound on higher derivatives we follow the same method for example in 
Lemma [l 2 . 2 1 1 and get 



\DTL''ia^a.o_ J\<Ck\d\ 5' 

(2) On the region currently considered, \9'\ is equivalent to \9\ so it suffices to 
compare with \9'\. 

\e'\ = x-H{x-^ht-iY^^ 

^ |6i'|^G(t, h) = x-'^t^ix-^ht ~ l)5^G(i, h) 
Note on the region currently considered we have x > d, (— ^o) > ^2 f^nd h, t 
bounded so the above expression is bounded for G bounded. Thus we have 
X > 5: for all N and all a 

□ 

As an immediate result of this we have: 

Lemma 12.33. U ^-.-.a ig smooth on x ^ Q. 

Lemma 12.34. \J^c"-.2 ^ x'-C^y + x'+C^y for ?/„ 7^ 0. 

Proof. Same proof as for Lemma |12.10[ the main idea is that the integrands in 
consideration have growth in \0'\ ^ \6\ but we can reduce arbitrarily the growth in 
1^1 when y„ ^ 0. □ 

12.4.2. Region F^fj"^ .^. 

(U^c.-,2^^) ^ lini(2^)-« /" f f^^o,^ u^xix,y)m{ee)e'y''dxdyd0 



lim(27r)-" / / / flpcoo wxe'y^dxdyde 



e-i-O J J J C,-,2 

Recall the structure of the approximate solution on this region : 
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where 



'^z^^ [l — (1 — z)^/"^] and C^c -,! jg smooth in z and s. 



Lemma 12.35. U ^.-.2 jg smooth on x ^ 



Proof. 



\0n 



so 



(1 - x)(^?n) exp(-<^)aFO^,^^ < V 4 

On the other hand on a; 7^ we have rapid decreasing in \9'\ coming from the term 
exp(0): Using 1 - (1 - x)^/^ > cx for X G [0,1] for some c > to obtain a lower 
bound for: 

l^nKl^f - 1)'/' (1 - (1 - ^)'/') > \9n\{\e'\' - If/'Z 

^\9^\{\e'\'^i)'/'x{\§'\^-i)-' 

= \9n\i\0'\'-iy^^x = i\e'\^-\9„\^)x 
> {\9'\'-il-dif\e'f)x^di\9'\x 

^ |6l'|^apo^,oo ^ < Cat; ViV e N 
Same result if we consider derivatives in 9j and dx- From here we get the smoothness 



on a; 7^ (see e.g Lemma 12.3) 



□ 

Lemma 12.36. U ^,-,2 smooth Un o-nd in fact vanishes to infinite order at 
X ^0 for yri ^ 0. 



we have: 



Proof. Same reasoning as in Lemma 12.29 to get for L such that *L = — -r- ^ yj 



\y\' 



ide. 



{xD^YlL'^a^o.^ I <Cfe|0TV 4 



The proof then is the same in ( 12.31 ). The main idea is that the integrand's growth 



is in terms \9'\ ~ \9\ which we have control over using method of stationary phase 
when yn ^ 0. For the rapid decay in x we use the boundedness of x~^ by 

x-^^s\9^{\e'\^-\Yl^<s\9'f 



□ 
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13. Main results on wavefront set and conormality 



In this section we prove the main result for the paper which is Theorem |1.1| and 
Corollary |1 . 2| and Corollary 1 1 . 3| which give a description of the singularity behavior 
of the solution to ^ 



'PU ^ 

supp U d {yn> 0} 

U = x'^-diy) + x''-+^gix,y) + h{x,y); g{x,y)€C^ 



c,I?y);supp5 c {yn = 0} 



h e Hii'b.ioc^ ^^O'' some a G 



Here s_ (A) = ^ — y '\ — A are the indicial roots of the normal operator for P. The 
existence and uniqueness for such a solution is given by Lemma |5.2[ Theorem 
Corollary |1 .2| and Corollary 1.3 describe singularities of the solution in the interior 
in terms of wavefront and at the boundary in terms of conormality. These results 
are used to arrive at Corollary |1.4| to address the original goal of the paper: whether 
or not there are singularities in the 'shadow region'. 

In this section, we put together all of the lemmas in section[T2] to get results for the 
'approximate' solution U = F~^{ux): then for the exact solution U . Note that from 
its construction and the computation in section 



satisfies (see Theorem 13.2): 
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U=.F-i(".)ex^-C"(M+,2?;) 



'PU = f; f e C°°(M+ X W^) 

ll/„<0 S X C.j. y + X '*'Cx,y 



13.1. Wavefront set statement for the approximate solution U. The proof 
follows exactly from that for Theorem 4.1 in Friedlander. Wavefront of (which 
is the restriction on U to a; > (5) is a union of the following wavefront sets: 



• From Lemma 
S(<5) = 



12.20 



From Lemma 



12.15 



the wavefront set of U '^■+ is contained in S((5), where 
iS)^{ix,y,^,Tj):^^dAy0-^),V = dy{y9-(f>),de {y ■ 9 ~ <j>) ^ 0,x > 6,\0n\ > \0'\>o} 

i=|(Z3/2-Zo^/2)sgn0„. 

the wavefront set of [J^c,+ is contained in T,{d), where 
nS) = {{x,y,^,v)-^^dxiy0'^),V = dyiy9-<f>),dg (y • - 0) = 0, x ><5, |0„| > |^?'| > o} 

1(^3/2 _Zo^/^)sgna„. 
• From Lemma 12.7 the wavefront set of U^'^ " is contained in I]((52, (5), where 

ihJ) = |(x,?/,^,77) : ^ = dx (y ■ - (f^ ,1] = dy (^y ■ - (f^ ,dg (^y ■ - (f^ = 0,6 £ conesuppaF~Q,a; > 6^ 
with 4> = ^Z^^^ sgn9n and 



conesupp ayo 



> 0, (1 - Ko{S2)) \0n\ < \0'\ < (1 + ^l{S2)) \0n\}- 
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• From Lemma 12.27 wavcfront set of U'^c - i is contained in S((5i,(5), where 
^i^i^S) ^ ^^{x,y,^,r]) : ^ ^ (y ■ e - }^ ,7] ^ dy (^y ■ e ~ (fj ,d0 (^y ■ e - ^'j ^ 0,9 € conesupp ap^^ ^ , a; 

with 4> = 1^^^^ sgn6'„ and 

conesupp of-,^ = {\e\ > 0, (1 - 2Si)\e'\ < |0„| < \e'\}. 

In short 

(97) WF(U5) C S((5)Ul]((5i,(5)Ul]((52,'5) 

I]((5) is independent of 6i and 62- Denote 

((5i,,52)eA = {(s,t):0<s,t<i} 
We have as ^2 ^ then ki{62), K2{S2) and 

fl S](<5i,5)=fl EiS2,S) 

A A 

= { (a^, y, e, r/) : C = a. • - ^) , r; - • - 0) , (y • - 0) - 0, l^'l = |0„| , |0| > 0, X ><5} 

All the first order derivatives of |^o^^ — 4' ~ 4> when |(?„| = \6'\ so 
fl E(<5i,<5)=f £(52,(5) CS(<5) 

A A 

This is in fact the boundary of 'S{5) which projects on the glancing bicharacteristics. 



Thus ( 97 1 can be replaced by 

WF{(Js) C S(5). 
Since this holds for all (5 > we have 
Theorem 13.1. With <j) = j{Z^/^ - Z^^^) sgn6'„, 

WF{[J) cJ: = [ix,y,^,r,):^^d,{y9~(f>),r,^dy {y0-^),de (y • - 0) = 0, x > 0, |0„| > \0'\>o] 

13.2. Conormality statement for the approximate solution U. 
Theorem 13.2. 

On y„ < 0, U G x'-C^y + x'+C^y. 

On yn > and if we are close enough to the boundary {x — 0}, then we have the 
same conclusion i.e. U G x^-C^y + x'^+C'^y 

Proof. The proof is a combination of all the results on 'classical conormality' (i.e 
in x'^-C^y + x^+C^y) proved in section [T2| 

On every region considered in section |l2| we have conormality on y„ < 0. 
For yn > 0, on each region in section [12[ either we have smoothness (with infinite 
order vanishing at a: = 0) or 'classical cornomality' behavior when we are close 
enough to the boundary {a; = 0}. □ 
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As a result of Theorem 



13^ U = T-\u^) e x''-C°iRt,V'y) satisfies 

PU = f; f e C°°(E+ X m;;) 



13.3. Conormality and wavefront set statement for the exact solution U. 



Remark 13.3. Conormal in this proof refers to 'classical cornormality' i.e x^^C. 



- poo 
x,y 



Proof for Theorem \1.1\ and Theorem 

Let U solve (i.e the exact solution to) the original problem ([2| (See Lemma 5.2 for 
the more information on boundary condition). 

'PU = 
supp U C {y„ > 0} 

U = x'~d{y) + x'-+^g{x,y) + h{x, y); g{x,y) e C°°(E„I?;); supp.g C = 0} 

h € Hl'^y^^\ior some a e M 



From its construction and Lemma 
solves 



13.2 



we have U = J^H"^) e a;"-C"(M+, ) 



PU = f; f e C°°{R+ X m;;) 
^"^-u|^^„ = %) 

We would like to to consider the cut-off of the approximate solution as follows: Let 
X be a smooth function such that x(t) = 1 for t < — 1 and = for t > — i. We 
have 

^(x(2/„)U) = x(2/«)PU + [P,x(y„)]u 

Note supp [P, x(yn)] C { — 1 < t < — i} i.e away from y„ = 0. On this region 
(y„ < 0) we have [P, x(yn)]U is conormal. Hence 



p (1 - x(y«)) U = PU - P (x(y„)U) G x^-CZ + ^'^CZ 



P[ (1 - xivn)) U j = f e x-C- + X^^C^y 
In fact (1 - xiVn)) U solves ^ x-'- (1 - xiVn)) U \^^^ ^ S{y) 

(1 - xiVn)) U |,„<-r= 



Consider the difference 
W solves 



= [/ - (1 - x(yn)) U 



Ply -f e X'-C^y + X'+C^y 



x-'-Wl = 

\x—0 







ly„<-i 

By the well-posedness result of Vasy [Theorem 8.11 [TT]], W is also conormal. On 
the other hand, by Theorem |13.2[ U is conormal on ?/„ > and if close enough the 
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boundary {x — 0}. Thus we have proved theorem 1.2 which says on ?/„ > and if 
we are close enough the boundary {x — 0} then 



U ex^-C^ + x^+C^ 



On the other hand, from the form of boundary condition satisfied by C/ in ([2]), on 
2/„ 7^ 0, [/ € H^'^ for some a € JR. This means that we cannot have the term with 
x'^^ in the expression of U on this region, therefore on ?/„ > and close enough to 
the boundary {x — 0} we have: 



U e x-^+C^ 



hence we have proved Theorem 1.3 



Also from this, since W is smooth on a; ^ 0, away from the boundary the singular- 
ities U are the same as those of U; hence we have proved Theorem |1.1| □ 

As an immediate consequence of Theorem |1.2| we have Corollary |1.3[ 



Proof for Corollary 1.3 For j/„ > and close enough the boundary {x — 0}, 
U e x^-C^y + x^+C^y. On the other hand, from the form of boundary condition 
satisfied by C/ in (jij), on j/„ ^ 0, U G H^'^ for some a G M. This means that on 
2/„ 7^ when we are close enough to the boundary, 

rl,00 



□ 
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